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ABSTRACT 

Analog vibration I'ontrol techniques are reviewed 
and are compare with digital techniques. The 
advantages of the digital methods over the analog 
methods are demonstrated. The following topics 
are covered: 

(1) Methods of computer -controlled random 
vibration and reverberation acoustic 
testing. 

(2) Methods of computer -controlled sinewave 
vibration testing. 

(3) Methods of computer-controlled shock 
testing. 

Basic concepts are stressed rather than specific 
techniques or equipment. General algorithms are 
described in the form of block diagrams and flow 
diagrams. Specific problems and potential prob- 
lems are discussed. The material is computer 
sciences oriented but is kept at a level that facili- 
tates an understanding of the basic concepts of 
computer-controlled induced environmental test 
systems. An introduction to computers is included 
as an appendix for a review of computer technology 
as it exists in early 1975. 

I. INTRODUCTION AND GENERAL BACKGROUND 

Before the introduction of computer-controlled test 
systems for induced vibration environments, 
industry used only closed-loop analog techniques 
for environmental vibration test control. Figure 1 
is a block diagram of a conventional closed-loop 
vibration analog test system. A power amplifier 
is used to drive an electrodynamic vibration 
exciter (shaker). Accelerometers or other trans- 
ducers mounted on a test fixture or test specimen 
send ac signals to the conditioning or charge 
amplifiers. The output signals from the charge 
amplifiers arc proportional to the acceleration at 
the accelerometer mounting points. The signals 
pass through analog servo control electronics, and 


this analCj, electronic device continuously corrects 
for the response spectrum by correcting the input 
level and frequency to the desired test 
specification. 

With the advent of the first hardwired special- 
purpose computers for signal processing, digital 
control of induced environmental testing was seri- 
ously considered. It was the development of these 
analyzers that has led to computer-controlled test- 
ing. The first computer-controlled test system 
was an attempt to control random-type vibration 
tests. 

The vibration test laboratory has been a natural 
place to apply automated digital signal processing 
and control techniques. The highly specialized 
analog instruments that have traditionally been 
used are generally limited in performance, diffi- 
cult to set and operate, expensive to maintain, and 
are prone to early obsolescence. The first 
attempts to overcome some of these deficiencies 
with digital techniques used large digital computers. 
Technical analysis performance was much 
improved, and the cost of processing large amounts 
of data^as decreased significantly, but other 
problems were introduced. The centralized data 
processing facility is inherently remote from the 
vibration test laboratory and must be manned by an 
army of specialists. To keep costs reasonable, 
data must be processed off-line, often resulting in 
turnaround times of 2 or 3 days. The remoteness 
and the off-line processing requirement prevent 
the operator from monitoring and interacting with 
the test while it is being performed. Attempts 
have been^made to overcome these deficiencies 
with use of remote time- shared terminals and high- 
rate data transmission links. The high coat of 
these techniques limits their use to the very large 
test facilities. 

The next generation of signal prc'^ssing equipment, 
first available commercially about six years ago, 
was essentially a special-purpose hardwired com- 
puter dedicated to performing specific functions 
rapidly and economically. The cost ^as low enough 


This paper presents the results of ^ne phase of research carried out at the Jet Propulsion Laboratory, 
California Institute of Technology, under Contract NAS 7-100, sponsored by the National Aeronautics 
and Space Administration. 

©Copyright 1974 by P. Chapman. 


1 


and the equipment small enough so that it could 
be put back in the laboratory under direct control 
of the test operator. Specialized computer oper- 
ators were not required because the instruments 
had a familiar type of control conoole» minimizing 
the need for the operator to understand the intri- 
cacies of digital processing. The hard-wired 
approach was necessary to achieve the required 
speed of processing economically, but had the dis- 
advantage of being highly specialized and not flex- 
ible. Cost considerations severly limited the 
functions that were designed into the machine and 
thus limited its usefulness. 

The development of the low-cost minicomputer has 
allowed the newest generation of digital vibration 
test instrumentation to combine the best aspects of 
these two approaches into one economical, versa- 
tile. high performance system (computers are 
discussed in the Appendix to this paper). This 
system can replace most, if not all, of the analog 
instrumentation in the typical vibration test 
laboratory. 

The basic concept of the computer -controlled 
closed-loop vibration test system is shown in 
Fig. 2. The purpose of this paper is to discuss 
the concepts associated with Fig. 2 with regard to 
computer-controlled random vibration and rever- 
beration acoustic testing, sinewave vibration test- 
ing. and computer-controlled transient (shock) 
testing. Initially, a brief review of the existing 
analog systems (Fig. 1) will be stated; throughout 
the paper, analog techniques will be compared to 
the digital techniques and advantages or disadvan- 
tages of the two techniques will be discussed. 

Basic concepts will be stressed rather than spe- 
cific techniques or equipment. General algorithms 
will be described in the form of block diagrams 
flow diagrams. Specific problems and poten- 
tial problems will be discussed. 

II, BASIC CONCEPTS FOR BOTH ANALOG AND 
COMPUTER-CONTROLLED TEST SYSTEMS 

A, Power Amplifiers, Shakers and Conditioning 
Amplifiers 

Figures I and 2 both illustrate a power amplifier 
driving an electrodynamic (or an electrohydraulic) 
shaker (exciter). The shaker generally has some 
type of fixture attached to it and. in turn, the fix- 
ture holds the test specimen. The forcing function 
to the shaker is either a random voltage signal, a 
swept sinewave. or in the case of shock testing, a 
transient signal. Accelerometers or other trans- 
ducers such as strain gages respond to th*^ dynamic 
environment. The response of the test item is 
indicated only at the points of attachment of the 
transducers and no assumptions shotdd be made 
about the test item dynamics except at the trans- 
ducer points. 

The signalo generated by the transducers are pro- 
portional to acceleration or strain and must be 
conditionea by special amplifiers. These ampli- 
fiers are analog-type devices which produce volt- 
ages proportional to acceleration or strain. The 
gain and sensitivity settings of these amplifiers 
can be controlled by multiturn potentiometers for 


the conventional conditioning amplifiers or by the 
computer in the case of some computer -controlled 
test systems. In either case, the signals coming 
out of the conditioning amplifiers are analog sig- 
nals, and how these signals are processed for test 
control purposes depends upon whether the test 
system is of the type illustrated in Fig. 1 or 
Fig. 2. However, the end result is to control the 
shaker in such a manner that the test specification 
is met. Both types of systems sense the response, 
compare this response to a desired reference, and 
then correct the output signal to the shaker such 
that the new response approaches the desired test 
specification. In the case of transient testing, the 
test is open loop rather than closed loop. A trans- 
fer function is computed between the input to the 
power amplifier and the response as seen at the 
output of the conditioning amplifier. The desired 
transient is mathematically and physically modified 
by the calculated transfer function in such a man- 
ner that when this modified transient is sent to the 
shaker, the response will, in fact, be the desired 
pulse in the time domain or the desired shock 
respoase spectrum in the frequency domain. 

An acoustic reverberation chamber can be substi- 
tuted for the shaker and a shaped random acoustic 
spectrum can be controlled on a closed-loop basis 
using the basic concepts of either Fig. 1 or Fig. 2. 

B. Analog Electronic Servos 

Figure 1 indicates a sine servo system or random 
automatic equalizer /analyzer. The sine servo 
system uses negative feedback to maintain the 
shaker control signal at the desired operating level. 
The transducer signal is rectified (Fig. 3) and the 
resulting DC signal is used to control the amplitude 
of an independent sinewave generator via a com- 
pressor circuit, A compressor speed control 
circuit provides a means of adjusting the .rate at 
which the servo will correct for changes within the 
control loop, either automatically or manually. 

The DC output from the compressor speed control 
is used to control the amplitude of the sinewave 
signal passing through the compressor. The Inde- 
pendent sinewave generator signal is heterodyned 
with a high-frequency carrier and single- sideband 
techniques and/or heterodyne techniques are used 
to translate the frequency higher in order to facili- 
tate the ease of operating on the sinewave signal. 
The high-frequency signal then pasr^es through the 
compressor circuit where its amplitude comes 
under control of the cot.ipressor speed control cir- 
cuit. The compressor o itput is then demodulated 
to provide the lower -frequency servo-controlled 
test signal output to the power amplifier which in 
turn drives the shaker to the test specification. 

For random excitation control the output of the 
conditioning amplifier is broken up into narrow- 
band frequency increments by a contiguous series 
of bandpass filters. The output of these filters is 
individually rectified and the resulting DC voltage 
is used to control the gain of individual amplifiers 
associated with each bandpass filter. 

For random excitation control (Fig. 4). the outputs 
of the conditioning amplifier and the noise gener- 
ator go to frequency translators, which beat the 
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random noise signals to a higher manageable 
frequency (as in the case of sinewave control). 

This higher -frequency spectrum is applied to a 
bank of contiguous oandpass filters. In the for- 
ward loop, the noise generator signal is applied 
to variable gain amplifiers associated with each 
bandpass filter. The gain of each of these' ampli- 
fiers is controlled by the detected signals from an 
identical set of contiguous filters in the feedback 
loop. The spectrum shape is initially determined 
I by a set of slide wires associated with each 

variable -gain amplifier. The output from each 
variable -gain amplifier is summed and demodu- 
lated such that the resulting signal approaches the 
test specification spectrum. 

Digital Electronic Servos 

\ It can be seen by comparing Figs. 1 and 2 that 

portions of the analog and digital test systems are 
: the same. In the analog case, knobs are set by 

the test operator that will allow the servo to 
implement the test specification. In the digital 
system, the test specification must be inserted 
f into the computer in the form of instruction codes. 

These codes arc typically transmitted to the com- 
put*-r by an operator /system interactive conversa- 
tional language and some minimal knob setting or 
switch selection. It is in this respect then that 
the two types of servos conceptually depart in 
similarity. The purpose of the conversational 
language program is to obtain, by means of a tele- 
type or other type computer interactive terminal, 

I all the necessary information needed to perform 

an induced environmental test to any given test 
specification. These software programs are gen- 
erally written in assembly language, and it is the 
means by which the systems operator communi- 
cates with the digital control system before, dur- 
ing, and after the environmental test. The 
conversational language is written as a sequence 
of subprograms. 

The pre-test portion of the program consists of a 
sequence of questions asked by the computer 
which must be answered by the systems operator. 
The questions are asked and answered in a speci- 
fied format. If a question is answered incorrectly, 
an error message is generated by the program and 
the systems operator is allowed to correct his 
mistake. The pre-test portion of the conversational 
language program is divided into subcategories: 

1. general Isgt Infgr.matifin. The computer 
might ask for the test title, number, and date, 
test engineer, wideband test level, and test dur- 
ation. The questions are typed out by the 
computer -controlled teletype or graphic terminal. 
The systems operator answers each question. The 
information obtained is tested for format errors, 
then stored in memory for later use during the 
actual test. The next question is then formulated 
by the computer. 

2. J ^requency i^rofile. Specifications for an 
induced vibration environment test are given in 
the form of a frequency profile versus amplitude. 
When displayed, the prohle can be viewed on a 
linear-linear, log-linear or log-log coordinate 


system. For random-type tests, the conversa- 
tional language program first asks the systems 
operator for all breakpoint frequencies in any 
order and then obtains a spectral or slope value 
for each range of the frequency profile. The 
acquired information is tested for format consis- 
tency and a spectral value is calculated for m 
spectral lines and n breakpoint frequencies over a 
bandwidth of p Hertz, Parameters m, n and p are 
variable. This information becomes the reference 
spectrum used by the digital system for control. 

At this time, a printout of the calculated wideband 
spectrum level is made. For sinewave testing, 
the program asks for displacement, velocity and/ 
or acceleration levels over the desired frequency 
ranges. For transient shock testing, the program 
asks for shock response data information, 

3. Time Profile. This section of the conversa- 
tional language allows the systems operator to 1 

specify the startup a* d shutdown times and how | 

long he wants to control the test at partial and full | 

levels. The systems operator can then choose to | 

obtain a paper tape output and/or a printout of all 1 

information and test Instructions to this point, I 

prior to the actual test. I 

.1 

The post -test conversational language includes the f 

capability to display, plot, and print out the test I 

results and test analyses in a predetermined for- I 

mat, which completely documents the environ- | 

mental test. In addition, the identical test can be | 

conducted again either immediately or any time | 

thereafter, since a permanent record of that par- ! 

ticular test specification exists on paper tape. | 

The control program is generally written in ? 

assembly language. The program uses the infor- I 

mation obtained by the conversational language ; 

program to perform all real-time computations i 

needed for system control. | 

i 

The digital control software determines whether 

the test specifications are presented via paper ^ 

tape or from the conversational program. Once 

the test specification data is obtained by either 

mode, the control program enters an initialize j 

subroutine. This subroutine initializes the com- | 

puter by zeroing all the input /output data blocks ^ 

and sets up all necessary counts; then the system 

is ready to start the test at the operator's 

command. 

The test specimen response signal is detected by 
an accelerometer mounted either directly on the <* 

test item or on a fixture interfacing the test item 
to the shaker. The accelerometer signal is con- 
ditioned by a transducer charge amplifier. The 
output of the charge amplifier is a voltage propor- 
tional to the test item response acceleration at the 
mounting point of the accelerometer. 

The charge amplifier voltage is the input signal to 
a computer gain controlled instrumentation differ- 
ential amplifier. This amplifier is used to isolate i 

the digital system from the shaker system and to 
allow gain adjustments for full-scale analog inputs 
to the analog- to-digital converter. 


3 



The output of the instrumentation amplifier is 
coupled to a low-pass anti ••aliasing filter whose 
cutoff frci^ucncy is act at leaat half the sampling 
frequency of the analog-to-digital (A/D) converter. 
This satisfies sampling theorem requirements. 

The filter is used in the Butterworth response 
mode with a rolloff of sufficiently steep slope. 

An ana)og>to>digital converter samples the output 
of the low-pass anti-aliasing filter at a rate deter- 
mined by the test specification. The digitized 
information is processed by the general-purpose 
computer and the time series analyzer. The input 
digitized data is compared to the reference test 
specification inserted in memory via the conver- 
sational language program prior to the actual test. 
Corrections are made by the computer, and the 
new forcing function signal in digital form is sent 
to the digital-to-analog (D/A) converter where it is 
converted back to an analog signal. 

The continuous analog signal from the D/A con- 
verter is passed through an output filter whose 
characteristics and cutoff frequency are identical 
to the input anti-aliasing filter. The purpose of 
the output filter is to assure a test frequency 
spectrum as specified by the test specification 
and to make the output of the D/A converter a 
c'^ntinuous signal. 

The filtered analog signal is conditioned by a 
computer-controlled line driving amplifier whose 
gain can be changed to produce the required drive 
voltage signal necessary to drive the vibration 
system power amplifier. 

The power amplifier drives the vibration exciter 
(shaker) with the modified test specimen response 
signal such that the new response approaches the 
desired test specification. The convergence time 
and control system response are determined by 
the mechanical dynamics of the test specimen and 
the algorithm used to determnie the error signal. 
The algorithm constants can be changed by the 
systems operator, 

III. WHY DIGITAL CONTROL OF INDUCED 
ENVIRONMENTAL TESTS? 

The rationale for developing the first random test 
system was based upon the following ideal 
advantages^ U: 

(1) The test specimen response can be 

"instantly" displayed on a lirge cathode- 
ray tube. Statistical erroi will be negli- 
gible and the confidence level will be high. 

A quick-look analysis of the response spec- 
trum is essential, particularly as the 
total test time decreases. For example, 
the total test time for the Mariner Mars 
1969 flight spacecraft was 30 sec. Future 
random noise vibration tests are expected 
to be even shorter. Most conventional con- 
trol systems can equalize the shaker and 
test specimen reasonably fast, but cannot 
adequately display the total response spec- 
trum in a short period. 


(2) For complete automation, a prepro- 
grammed test specification can be tape-fed 
to the compute 1 , Te»l parameters (levels, 
test time, etc. ) can be verified prior to 
the actual test. Human error will be 
eliminated. 

(3) The programmed tape input and the use of 
pseudorandom noise (to be discussed) 
ensures test repeatability. All computer- 
controlled environmental laboratories could 
generate exact test environments. This 
would be a step toward test standardization. 

(4) All analog systems require periodic cali- 
bration because of inherent drift. A digital 
system will minimize these drift and cali- 
bration problems. Test results will be 
more accurate and reliable. 

(5) The quick-look readout display capability 
of a digital system will allow the observa- 
tion of test specimen irregularities. The 
computer system will sense the irregularity 
and abort the test quickly, if necessary. 

The criteria for test abortion will be pre- 
programmed into the computer before 

the actual test. 

(6) Protection of the test specimen (test hard- 
ware) against transient acceleration and 
excess steady- state acceleration, velocity, 
and displacement is of prime importance 
during all types of vibration testing The 
digital system can be programmed to 
monitor these conditions. Should abnormal 
conditions arise with respect to the pro- 
grammed "normal" inputs, the computer 
will abort the test in a controlled manner to 
protect the lest specimen against damage. 

In addition, the computer could monitor 

the force /current requirements of the 
vibration exciter, which in turn will provide 
test specimen protection against the 
occurrence of abnormal conditions in the 
closed-loop control system. 

(7) Changes in control philosophy and improve- 
ments in testing methods can be easily 
implemented with appropriate changes to 
the computer program. 

(8) All the important information necessary to 
obtain a complete test report would be con- 
tained within the nondestructive memory 

of the computer. The test report would be 
generated after the test by interrogating 
the computer. Peripheral equipment 
associated with the digital system would 
print out the results of the test in a pre- 
determined format, thus eliminating costly 
computer and analysis time. 

(9) The computer control system can also be 
used for shock pulse and shock spectrum 
analysis and synthesis. Shock testing with 
electrodynamic shakers will be more prac- 
tical with an "on line" computer. 
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The first random eysteme were constructed, and 
these concepts proved to be real advantages. As 
the sine and shock capabilities were added as fea- 
tures to the computer -controlled test systems, 
still more advantages over the analog systems 
were obvious. These advantages were mainly with 
the fact that the concept of the test system could 
be changed by changing or adding software in 
order to meet special-purpose test requirements. 
These advantages will become obvious in later 
sections of the paper. 

IV. GENERAL FUNCTIONAL REQUIREMENTS 
FOR COMPUTER-CONTROLLED VIBRATION 
TEST SYSTEMS 

The functions that the test system is to perform 
can be classified into the following major 
categories(^): 

(1) Test Signal Generation. 

(2) Data Acquisition and Storage. 

(3) Data Analysis. 

(4) Servo Control of Tests. 

(5) Test Documentation. 

Each of these categories will be described below. 

A . Test Signal Generation 

This is a relatively new function that digital sys- 
tems have taken on. It Is now possible to econom- 
ically generate all the test signals for induced 
vibration tests. This includes random, sine, and 
shock signals with the required parameter 
variations. 

1. Random Signal Generation . There arc several 
methods of generating random noise signals using 
digital techniques. First thoughts on this topic as 
applied to early proposed computer -controlled test 
systems for random noise testing are described in 
Ref. 1, and are repeated here. 

To be better able to imder stand the concepts of 
random noise generation for computer-controlled 
test systems, it is necessary to describe some 
basic electronic components associated with digital 
noise generation. One such component is the flip- 
flop or bistable multivibrator . 

A flip-flop is a device with two outputs, 0 and 1, 
mutually exclusive (bistable). The device produces 
no change in output signal unless a signal (pulse) is 
applied to its input in accordance with Fig. 5. 
Assume the 0 and 1 outputs are on and off 
respectively. A pulse applied to the S or "set" 
input terminal will turn on the I output (1 high) and 
turn off the 0 output (0 low). A pulse applied to the 
R or "reset" terminal \/ill turn the 1 output off and 
the 0 output on. If the 1 output is high and a pulse 
is applied to the S terminal, or if the 0 output is 


high and a pulse is applied to the R terminal, no 
change of output state will occur. If input signals 
or clock pulses (equally spaced pulses) are applied 
to the T or "trigger" input of the flip- Hop, the 
output wiU switch to its other state as shown 
(Fig. 5). 

Another electronic device which can be used for 
noise generaiion is the voltage comparator. The 
voltage comparator is a high-gain, differential- 
input, single -ended output amplifier. The purpose 
of this device is to compare a signal voltage on 
one input with a reference voltage on the other and 
to produce a logic 1 at the output when the signal 
voltage is equal to or greater than the reference 
voltage. When the signal voltage is less than the 
reference voltage, the comparator output is low 
(logic 0). 

With the flip-flop and the comparatcr in mind, it 
is now possible to discuss binary random noise 
gene ration. 

a. Binary random noise generation . Figure 6 
demonstrates conceptually a method of generating 
what is commonly referred to as binary random 
noise. A wideband noise voltage from a noise diode 
is compared with a fixed voltage in such a manner 
that the output of the comparator triggers the first 
flip-flop somewhat randomly at a rather high fre- 
quency. Both outputs of the first flip-flop go to the 
S and R inputs of the second flip-flop. In addition, 
the T input of the second flip-flop is connected to a 
pulse generator (clock) that is running 20 to 100 
times higher in frequency than the highest desired 
frequency component to be generated. The output 
of the first flip-flop and the output of the clock are 
essentially independent events. The output of the 
second flip-flop is a telegraph wave of several 
different pulse widths with some DC voltage compo- 
nent. This telegraph wave is applied to a precision 
switch. The output of the switch is voltage -clamped 
to ±5 V. The output waveform from the clamp 
circuit theoretically has zero mean. It is impos- 
sible to predict with absolute certainty whether the 
output waveform at any time is +5 V or -5 V. There 
is a 50% probability that the binary waveform will 
change state on any given clock pulsed Clearly, 
this type of electronic circuit and its output provide 
an electronic analogy to a coin toss experiment. 


The spectrum c&nd amplitude distribution character- 
istics^ are illurtrated in Fig. 7. The binary 
random pulses have a sin^ x/x^ continuous - 
frequency spectrum. The spectrum is determined 
by the constant amplitude (±5 V) of the binary 
pulses and the clock frequency F^ (or period T), 
as can be seen in Figure 7a. If the binary pulses 
are passed through a low-pass filter whose cutoff 
frequency is approximately 0.45 F^, the frequency 
spectrum of Figure 7b can be generated(^). The 
slope of the curve above 0.45 F^ is determined 
by the order of the low-pass filter. The bandwidth 
of the white noise portion of the spectrum is 


^A good reference on random noise characteristics is W. Tustin, "Random Vibration Potent Test Tool," 
EDN Magazine, Cahners Publishing Co., April 1967. 
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entirely dependent upon the clock frequency. The 
emplitude distribution of the filtered pulses 
appears Gaussian, as seen in Fig. 7c 


is not a trivial matter, but tables can be found In 
the literature covering up to at least 33 stages of 
shift register 


The random noise produced by the method 
described above is indistinguishable from random 
noise generated by more conventional techniques, 
provided the clock frequency is much higher than 
the highest desired frequency component. 


b. Pseudorandom generation. To under • 

stan^ the next type of digital noise generation, it is 
necessary to briefly define two additional elec- 
tronic components: the shift register and the 
exclusive -OR gate. 


A shift register is formed from a series of flip- 
flop stages which have been modified so that clock 
pulses applied to the R input terminals cause a 
binary pulse (bit) to be shifted through the stages. 
For example, assume that the first flip-flop of 
Fig. 8 has been set to represent the logic 1 state 
(I output high) and all other flip- flops are in the 
logic 0 state (0 output high). When a clock pulse 
is applied to the R or reset line, it will not affect 
those stages at logic 0, but will clear (reset) the 
first flip-flop at logic 1. When the first flip-flop 
is cleared, it produces an output pulse which 
enters the delay circuit, usually part of the flip- 
flop. After the clearing pulse duration is termi- 
nated, the pulse signal leaves the delay circuit and 
enters the T or trigger input of the second stage, 
flipping the second stage to logic 1. A second 
clock pulse entered on the reset line will clear the 
second flip-flop (all other flip-flops are in the 0 
state), sending a signal into the third flip-flop, 
which in turn trips that flip-flop to logic I after the 
clock pulse duration has terminated. In short, 
successive clock puises progressively shift the 
binary bit to the right through many flip-flop 
stages as are contained in the shift register. 

An exclusive -OR gate is a two -terminal input and 
single -terminal output that gives an output if, and 
only if, one input is high (logic 1) and the other 
input is low (logic 0). The exclusive -OR symbol 
and the logic truth table arc shown in Fig. 9. 

With the concepts of the shift register and the 
exclusive -OR gate in mind, it is appropriate to 
discuss the second, and probably the most impor- 
tant, type of digital noise generation: pseudo- 
random noise generation. 

A shift register and an exclusive -OR gate, if con- 
nected as shown in Fig. 10, will generate a spe- 
cial binary random sequence, provided that 
initially not all the flip-flops are set to logic 0 (0 
o\ tput high). If this special binary sequence is 
such that the shift register generates the maximum 
number of logic 1 and logic 0 combinations possible 
for n stages of shift register, the binary sequence 
is called a maximum -length sequence (^). 

Although Fig. 10 shows a 4 -stage shift register, 
longer (more stages) shift registers are generally 
employed and more than one exclusive -OR gate 
coidd be used to generate a maximum-length 
sequence of binary pulses. Where the exclusive - 
OR gate should be connected to the shift register 


As seen in Fig, 10, the binary pulse train contains 
IS bits of all possible logic combinations (except 4 
logic Os) for a 4- stage (n = 4) shift register. In 
general, with an n-stage shift register, it is pos- 
sible to generate • 1 bits or clock periods 
before the sequence is repeated, provided the 
exclusive -OR gate is properly connected. The 
length of time the sequence lasts before it repeats 
is determined by M • (AT), where M = 2*^ - 1 and 
AT is the period of the clock pulse. For example, 
suppose a 2S-stage (n = 25) shift register and a 
200-kHz (AT = 0. 5 • 10'^ s) clock are used to gen- 
erate a maximum-length sequence. Then, 

M = 2^5 . 1 or M = 33, 554,431. The period 
M • (AT) is equal to 167.77 s or approximately 
2. 8 min. That is, the binary sequence would 
exactly repeat itself every 2, 8 min. (The sequence 
has a periodicity of 0. 006 Hz. ) 

The binary pulses from the maximum -length 
sequence generator appear to be quite similar to 
the binary pul»et» gcuerated from the binary ran- 
dom noise generator. In fact, one might say that 
the binary random pulses are a maximum-length 
binary sequence of infinite length. However, there 
is one important difference between the two types 
of pulse generation: the binary pulses from the 
maximum-length sequence generator are periodic 
and deterministic, in contrast to the binary random 
pulses, which are nonperiodic and probabilistic. 
Periodic pulses produce line spectra, whereas non- 
periodic pulses produce continuous spectra (7), 

The line spectra characteristics of maximum- 
length sequences are illustrated in Fig. 11, Low- 
pass filtering of the binary sequence produces 
pseudorandom noise (’’random " noise that is peri- 
odic). There is no statistical variance associated 
with pseudorandom noise because it is a periodic 
deterministic signj . If the spectral lines are very 
close together, this type of noise has all of the 
desirable characteristics needed for vibration and 
acoustic environmental testing, including first- 
approximation Gaussian-amplitude distribution. 

But, most important of all, anyone can generate the 
same binary random chain and therefore the same 
noise signal, provided the code is known. 

Figure 12 illustrates an actual hardware pseudo- 
random noise generator. The random noise comes 
from the effects of the low-pass filter. A binary 
maximum-length sequence exists at TP3, This 
same sequence can be generated within a computer. 
That is, a special hardware register need not be 
built. The sequence can be generated using ^ soft- 
ware subroutine program and working registers 
within the control system computer. The results 
will be number states that will be either +1 or 0, and 
these states will vary in a pseudorandom manner. 
Figure 13 is a flow diagram of such a program. 
These numbers can be used to randomize a com- 
plex frequency spectrum to produce phases that 
yield random noise. In referring to Fig, 13, it 
can be seen as an example that two 16 -bit registers 
and the overflow single flip-flop C bit register are 
program-controlled in accordance with the flow 
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diagram. The new C bit (C^) is either a 1 or 0, 
depending upon the previous state of bit 32. The 
new bit shifted into the first stage of the first lb- 
bit register is the exclusive-OR of previous bits 32 
and 3. As the program is entered, bits 32 and 3 are 
exclusive -ORed. This results In two possibilities: 
either a logic 1 or a logic 0. At this time, the com- 
puter clock shifts the results of the exclusive-OR 
output into the first stage, i, e. , bit 1. Simultane- 
ously, all the other bits have a 50% probability of 
changing state. A test is then made on the C bit, 
and the program can then branch to another sub- 
routine depending upon whether the C bit was a 1 
or 0, In other words, the C bit (or any of the other 
31 bits) varies in binary state, either 1 or 0, in a 
pseudorandom manner. The resulting sequence 
will repeat every four days when clocked at a 
sufficient rate to yield 2000-Hz bandwidth random 
noise* 


Since M has units of octaves, N/60 octaves/sec can 
be multiplied by units of seconds to yield octaves 
only as in £q. (2). That is, 

M 5 octaves (3) 

where, for logarithmic sine sweeps, At * 

£q. (2) is of the form 

f. = f. , (4) 

I 1-1 ' 

where 

X = N • ln2/60 

N = sweep rate in octaves/min 
At = 


There are many computer algorithms available to 
generate random numbers which in turn can be 
used to randomize the phase of a complex 
frequency spectrum to produce random noise. 

As can be seen, then, there is clearly no problem 
in generating random noise signals using either 
special hardwired hardware or the control system 
computer itself. 

The philosophical question exists, of course, as to 
whether central processing time should be allo- 
cated to the task of noise generation or whether 
special hardwired hardware should assume this 
task independent of the computer. 

2, Sinewave Generation. As with random noise, 
sine wave generation can be accomplished either 
by the control system computer and a software 
algorithm or by special dedicated computer- 
controlled hardware such as a programmable 
function generator. 

a. Computer-generated sinewaves . Induced 
environmental sinewave testing is almost always 
accomplished by generating a logarithmic sine 
sweep as the forcing function. The sweep is 
generally specified in octaves per minute. The 
definition of an octave is 

In fj{. 

M = — — octaves (1) 

where M is the number of octaves in the frequency 
range f. to f^ and In is the natural logarithm 
operator* 

For logarithmic sweeps, the sweep rate is speci- 
fied in N octaves/min and this is the same as N/60 
octaves/sec. 

By multiplying both sides of Eq. (1) by In 2 and 
then exponentiating both sides, Eq, (1) becomes 


. In 2 
2 = M ® 


( 2 ) 


f. is the previous full cycle sinewave 
frequency, i.e., fj Hz 

f. is the next full cycle sinewave frequency, 

^ i.e., f^ Hz 

To implement Eq. (4) is certainly no strain using 
the system computer, except that exponentiation 
takes time and it turns out that some good approx- 
imations can make logarithmic sweeps easy and 
fast to implement. It can be seen that X is a 
constant. The function f^ can be approximated by 
an exponential series; i.e.. 


f. 

i 


= f 


i-1 


+ . 


1 + XAt + 


(V.At)^ 

“ZT" 


+ 


(XAt)^ 

“7T“ 


. (VAt)“ 
n! 


(5) 


Typically, N is in the order of from 1 to 4 octaves/ 
min over the frequency range from 5 to 2000 Hz. 

An inspection of the higher-or r terms of the 
series indicates that these terms add practically 
nothing to the calculation of f^« This can be seen 
in Table 1 for these two sweep rates and frequency 
limits: 


Table 1. Exponential series terms 


Frequency 

(Hz) 

Series 

term 

1 octave /min 
f. jXAt, Hz 

4 octaves /min 
f. ,XAt, Hz 

i- 1 

5 

1 

5 

S 

5 

2 

1. 16 X ;o’^ 

4.62 X 10'^ 

5 

3 

1. 33 X 10"^ 

4.27 X 10'® 

2000 

1 

2000 

2000 

2000 

2 

1. 16 X 10*^ 

4.62 X 10'^ 

2000 

3 

3, 34 X 10‘® 

5. 34 X 10'^ 

(Insignificants of higher-order terms of the 
exponential series) 
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# 
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Therefore, Eq, (5) is approximated by 

f. 5fi 1 (1 +\At) (6) 

But At = l/f._j 
Therefore, 

f. = r X ^ ^ upsweeps (7) 

and 

r = f._j - X for downsweeps (8) 

Equations (7) and (8) then are equations that could 
be implemented by the control system computer to 
generate a logarithmically swept sinewave in real- 
time. 

The actual generation can be implemented in the 
following manner. Assume an allocation of p core 
locations to store p values of a half sine function 
(0 to TT radians). The sine v^ilues can be generated 
and stored in these allocated core locations by 
performing tne following integer arithmetic 

TABLE. = SIN ^-jQ. I (9) 

where Q. = 1, 2, 3, . , . p, and p is the number of 
points (values) of a half sinewave. 

For a 16 -bit machine, the sinewave function can run 
over the numeric integer range from +32767 to 
-32767 withO being at 0% +32767 at +90*. 0 at 180*, 
-32767 at 270* and 0 at 360*. For a constant clock 
rate, the lower frequencies will have more points 
p than the higher frequencies. The clock rate can 
He calculated by specifying the minimum number of 
points p from ^he allocated core locations to be 
outputed. That is, 

CLOCK PERIOD = ( 10) 

where p is the number of points to be used to 
define the half sinewave at frequency T. 

In order to increase the frequency in accord with 
Eqs, (7) and (8), it is necessary to output the sine 
values in the allocated core locations at a faster 
rate, and If the clock that works the digital-to- 
analog (D/A) converter operates at a constant rate, 
this means that fewer values must be put out as the 
frequency fi increases. This means that the qual- 
ity of the sinewave suffers as the frequency 
increases. However, the fact that the D/A con- 
verter output is passed through a low-pass filter 
guarantees that even a very minimum number of 
values will yield a low distortion sinewave. As an 
example, choose a maximum of 409f values to 
define a hal' sinewave and a clock period of 4 ps 
(250 kHz). From Eq, (10), the low-frequency 
limit would be about 0. 3 Hz, and at 2000 Hz only 
62 points could be used to define the half sinewave. 
Owing to the fact that the function is passed through 
a low-pass filter, the resulting function would bo a 
clean sinewave of low distortion. 


It can be seen from Eq. (10) that a scheme to 
change the clock period could result in a variable - 
frequency sinewave that was always defined by the 
same number of points (values). This concept u 
difficult to implement with s. ^ wave coTiputer 
generation techniques but has been implemented 
with a special hardware sinewave generator 

b. Hardware sinewave generators . The first 
computer -controlled sinewave te^ system used 
a minimal amount of core storage since It v s 
strictly sine control. Core storage at that time 
was costing about $3500 for 4K words (K = 1028) of 
memory, so the manufacturer chose to design a 
special computer-controlled function generator to 
generate the sinewave forcing function, although 
additional "advantages" were listed (Ref. 8) to 
provide the rationale for this special hardware. 
This function generator produced a 128-step 
synthesized waveform across its frequency range 
by varying its clock frequency. That is, it 
employed a variable clock oscillator as the master 
clock, and the clock frequency was controlled by 
the control system computer. The variable clock 
technique v/as accomplished through implemen- 
tation of a digitally controlled phase-lock-loop to 
precisely control the master clock oscillator 
frequency relative to a highly stable fixed- 
frequency reference oscillator. The phase lock- 
loop was connmanded by a modulus term while the 
octave was commanded by an exponent term. 

These terms are related by 


where 

f^^^ is the output clock frequency 

= proportionality constant 
A = 0, 5 to 1 in 1024 steps 
N = 0 to 14 in 15 steps 

There are several types of programmable frequ • cy 
synthesizers suitable for sinewave generation in 
computer-controlled test systems in addition to the 
one described in Ref, 8 above. 

One type, the direct digital frequency syn^esizer, 
is shown in simplified form in Fig. 14 ' ^ The 

digital-to-analog converter and low-pass filter are 
the only sections of the system that are not built 
entirely with digital IC's. 

The phase accumulator generates the linearly 
increasing digital phase value of the sinusoid at a 
fixed output sample rate. This phase value 
advances for each output sample by an increment 
directly proportional to the frequency setting. The 
read-only-memory (ROM) is a sii * function table 
that converts the phase information provided by 
the accumulator into digital samples of a sinusoidal 
waveform. 
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These fAmplet ere fed to the D/A, which produces 
a staircase approximation of a sinewave in analog 
form. The iow-pass filter removes the higher- 
frequency sampling components, reducing the wave- 
form to a pure sinusoid. 

With reference to Fig, 14, the phase accumulator 
contains a frequency register. The purpose of the 
frequency register is to load and store from the 
control system computer the digital value of the 
instantaneous desired frequency, as given for 
example by Eq. (4). 

The accumulator consists of an N-bit binary full 
adder followed by an N-bit data register. At every 
clock pulse, the adder sums the values in the data 
register and the frequency register and updates 
the data register with this new sum, which repre- 
sents the phase of the sinusoid being generated. 
Since the accumulator increments at a fixed rate 
determined by the reference clock, the value in 
the frequency register is specifically the rate of 
change of phase. 

The accumulator overflows at modulo 2^, and this 
cycle of 2^ represents one cycle of the sinusoid. 

The relationship between N (number of accumu- 
lator bits), Fi (lowest frequency possible), and Fc 
(CLOCK frequency), is given by 



where F] represents the resolution of the 
frequency setting. 

The phase samples provided by the accumulator 
are used to address the ROM in which the sine 
function is stored. Since the BOM is a relatively 
expensive device, only the first quadrant values 
are stored. External circuitry is then used to 
extend these values for all quadrants. 

One version of the accumulator and ROM signal 
waveforms is illustrated in Fig. 15 ' ^ 

The most significant bit (MSB) of the accumulator 
will be at logical "0" for phase values in quadrants 
1 and 2 and logical "1" for values in quadrants 3 
and 4, This behavior represents the polarity of 
the output sinusoid and hence is called the SIGN bit. 
The next lower bit is logical "0" during quadrants 
1 and 3 and logical '*1" for quadrants 2 and 4; 
therefore it is called the QUADRANT bit. These 
top two bits specify one of four quadrants and are 
used to extend the first quadrant values of the ROM. 

Since there are K address bits into the ROM, the 
next lower K accumulator bics b.low the 
QUADRANT bit are used. The digital value of 
these bits during one sinusoidal cycle is plotted in 
Fig. 15, There are four sawtooth cycles, one per 
quadrant. Using these bits directly for the first 
quadrant will address the ROM correctly, but for 
the second quadrant these bits must be comple- 
mented so that the slope of the sawtooth is inverted, 
effectively addressing the ROM from 90* back to 
0*. Likewise, for the third quadrant, the address 


is not complemented and for the fourth it is. The 
digital value of the address to the ROM is shown in 
Fig, 15. This waveform represents correct phase 
addressing for a full-wave rectified or magnitude - 
only version of one sinusoid cycle. To accomplish 
this address complementing function, a K-bit paral- 
lel complementer circuit is required with comple- 
ment control provided by the QUACHANT bit. 

The digital output of the ROM is shown in Fig. 15, 
and this corresponds to the expected magnitude- 
only waveform. The sign or polarity information 
is CO* ained in the SIGN bit. Since conventional 
DA^ require straight binary or 2's complement 
digital input codes, this sign-and-magnitude ver- 
sion must be converted. This is accomplished 
with another complementer circuit which comple- 
ments the output of the ROM during quadrants 3 
and 4. Control of this f^inction is provided by the 
SIGN bit. 

The M bits from the amplitude complementer plus 
ihe SIGN bit comprise the M+1 bit word of the 
complete sinusoid ready for D/A conversion. It is 
converted and then passed to the low-pass filter. 

The output of the low-pass filter is a clean sine- 
wave which car. be used as the for cing function to 
the power amplifier. Its amplitude is modified by 
the control system as required to meet the test 
specification. 

In summary, then, it is possible to either generate 
the sinewave entirely by the control system com- 
puter or use special hardwired hardware controlled 
by the test system computer to synth esize the sine- 
wave driving signal. 

3. Transient Waveform Generatio n. Transient 
waveforms are used for computer-controlled shock 
testing using the test system power amplifier and 
shaker. These transients are generated by the 
F. ^t Fourier transform processor and test system 
computer*. There are several classes of tran- 
sients currently being used for tliis type of shock 
testing (12)^ some of which are derived from 
transfer function techniques 03) and described 
in a later section of this paper 

B, Signal Conditioning 

Analog signals, generally from multiple sources, 
must be conditioned, sampled, digitized, and 
stored in some form of mass memory, such as 
core, magnetic tape, disc, or multiples of these 
forms of storage. 

Figui-J 16 illustrates the functions required to 
convert a single-channel analog signal to digital 
form prior to entering the control system computer 
for processing or temporary storage, 

1. Analog Signal . The analog signal comes from 
the transouc^ signal conditioning amplifier. This 
analog signal is a voltage proportional to strain or 
acceleration. Tf t vibration or reverberation 
acoustic control, this signal's frequency content 
will usually be below 10 kHz, and in most instances 
the frequency content of concern is usually below 



f 


I 


2500 Hz» since the destructive properties of a 
forcing function to a structure or component are 
proportional to the nonrigid body displacements 
resulting from the forced vibration, and the dis- 
placements decrease with increasing frequency by 
the square inverse frequency. That is, for sinusoi- 
dal forcing functions, 


Peak- to -peak displacement 
where 


S 

0.0511 f 


I 


g = normalize 1 acceleration, and 


(13) 


f = frequency of the forcing function. 

At 2000 Hz, a l-g vibration results only in a 
theoretical displacement of the test item of about 
5 fiin. One might worry about these small dis- 
placements in the manufacture of semiconductor 
devices but certainly not spacecraft or aircraft 
structural members. Therefore, in most instances, 
the interest of concern bandwidth of the analog 
signal will be about 2000 Hz, although the analog 
signal is the result of a response at the mounting 
point of the transducer, and in most cases this 
response signal will contain frequencies greater 
than 2000 Hz. 


2. Programmable Amplifier /Attenuator . The 
amplitude of the analog signal can vary over per- 
haps a 40 dB (voltage) range in the case of a 
response signal or open-loop test, in contrast to a 
control channel. The purpose of the programmable 
amplifier /attenuator is to modify the analog signal 
amplitude in such a manner as to guarantee a 
known maximum (full-scale) input into the remain- 
der of the hardware illustrated in Fig. 16. If the 
devices illustrated in Fig. 16 are rated at ±10 V 
maximum full-scale input, it is necessary to guar- 
antee that, at least in the case of random noiso 
testing, the analog signal leaving the programmable 
amplifier/attenuation be no greater than about 2 V 
rms. This constraint is necessary in order to 
guarantee that the 3- sigma noise spikes get digi- 
tized and are taken into consideration in determin- 
ing the rms control level. 

For sinewave testing, the two-thirds down from 
full-scale restriction need not be observed if the 
control amplitude algorithm is based upon a peak 
detector concept. If the software detector concept 
is based upon an average algorithm, then it is 
necessary that no clipping take place in the data 
channel. 


As can be seen, for large-amplitude signals it may 
be necessary to attenuate, or in the case of low- 
amplitude analog signals it may be necessary to 
amplify, in order to meet the above constraints. 

The amount of attenuation or amplification is 
’’remembered" by the test system computer and is 
used as a scale factor in calculating desired as 
well as actual test levels. The initial setting of the 
amplifier /attenuator is determined by the system 
computer and is baaed upon transducer sensitivity 
and full-scale test levels specified by the test 
operator during the conversational mode (inputing 


the test specification) of the control software 
program. 

3. Programmable Low- Pass (Anti-Aliasing) 

Filter . The output of the programmable arnplifier/ 
attenuator must be passed through a low-pass 
filter prior to digitizing in order to minimize 
aliasing errors. The concepts of aliasing are 
discussed below. Let it suffice for now to state 
that it is necessary to sample the analog signal 
from the transducers and conditioning amplifiers 
at a rate that is at least twice the highest fre- 
quency of the analog signals and preferably four 
times the highest frequency contained within the 
analog signals. This can he guaranteed by low- 
passing the analog signals prior to digitizing. For 
acoustic random, and sine testing, the filter should 
have Quite rworth response; and its cutoff frequency 
as specified during the conversational mode of the 
control program should be 0 dB and not *he nor- 
mally accepted half-power point of -3 dB, The 
conventional transfer function of this type of filter 
is illustrated in Fig. 17 and the distribution of the 
poles of the filter are shown in the Fig. 18 S-plane 
diagram. 

The Butte rworth- type filter is a "maximally flat" 
amplitude response filter. This type of filter is 
generally used for data acquisition systems to pre- 
vent aliasing errors in sampled- data applications. 
The attenuation rate beyond conventional (-3 dB) 
cutoff frequency is -n6 dB /octave of frequency 
where n is the order (number of poles) of the filter. 
For acoustic and vibration control, n should be at 
least 6, giving an attenuation of at least -36 dB/ 
octave outside the filter passband. 

For shock testing, however, the Bessel filter (also 
known as linear phase) is more suited as the low- 
pass filter within the data channel of Fig. 16, 

Figure 19 indicates the transfer function of this 
type of filter, and F.|i?. 18 shows the distribution 
of poles in the complex frequency plane. The 
Bessel filter is a linear phase filter, and because 
of this linear phase characteristic, the filter 
approximates a constant time delay over its pass- 
band frequency range. Bessel filters pass tran- 
sient waveforms with a minimum of distortion and 
over shoot-to- step inputs. The maximum phase 
shift is -nw/2 radians, where n is the order 
(number of poles) of the filter. The conventional 
cutoff frequency f^. is defined as the frequency at 
which the phase ^nift is one -half of this value. For 
accurate delay, the cutoff frequency should be at 
least twice the maximum analog signal frequency. 
Table 2 lists the conventional -3 dB cutoff frequency 
as a function of n, the number of poles. 


Table 2, Bessel filter cutoff frequency 
characteristics 



2-pole 

4 -pole 

6 -pole 

8-pole 

-3 dB 

frequency 

0. 77 fc 

0. 67 fc 

0. 57 fc 

0. 50 fc 


10 



An '*ideal" filter is normally described in the 
literature as being one with flat response to a cut- 
off frequency, infinite attenuation beyond that 
frequency, and linear phase or constant delay over 
the passband. Such a hlter exists only in the minds 
of mathematicians. 

Ideal filters can be approximated by networks that 
have very large numbers of elements* Alternatively, 
they can be simulated U/ lengthy digital computa- 
tion. But any real filter will have a finite rate of 
cutoff and will never attenuate at an infinite rate. 
Also, it will have phase nonlinearities. Therefore, 
a real filter cannot multiply all of its different 
frequency components by an equal amount with 
constant delay. And it can't completely eliminate 
all the undesired signals. 

A real filter, then, will inaccurately multiply the 
many spectral components of a time-varying signal. 
The filter's output will differ from the input signal 
e ;n if its nominal bandwidth is greater than that 
of the input. To determine the effect of a filter on 
dynamic accuracy, one must know the spectral 
distribution of the incoming signal and the transfer 
function of the filter. Then, by multiplying the 
original spectrum by the transfer function of the 
filter, one can derive probable signal accuracy 
over a period of time for signals having varying 
amplitude and frequency components. 

4. Encoding. The encoding of the analog signal 
from the filter consists of a sample -and-hold (S/H) 
device and an analog-to- digital (A/D) converter. 

For multiple channel control or analyses, a multi- 
plexer would be inserted after the programmable 
hlter in Fig. 16 and would be controlled by the 
system clock. Each signal channel would have its 
own programmable amplifier /attenuator and low- 
pass filter, but would time share the encoder. 

Since the encoder takes a finite time to convert the 
data from analog to digital form, rapidly changing 
data can cause errors. If the data change during 
the sampling interval is greater than tne value of 
the least significant bit in the binary word, a 
fpurious output results. The solution is a circuit 
that samples the analog signal fast, then holds this 
value as a dc level for the full conversion time. 

This is the basis of the sample-and-hold circuit. 

Figure 20 illustrates the oasic idea of a sample and 
hold circuit. The basic sample-and-hold circuit is 
a normally open switch that closes for a very short 
time interval so that the analog voltage value can 
be stored in a memory element. This value is then 
held constant, independently of the input signal, 
while it is encoded. The memory element (often a 
cc.pacitor) is buffered from the output to minimize 
decay error. 

In the case of several control or anal/ses signal 
channels, and where the change in data voltage is 
more than the least bit during the multiplexer 
period, the data channel is sampled at the begin- 
ning of the multiplexer's normal sample interval. 
When the multiplexer rise-time is too slow and 
does not allow enough encoding time, the sample- 
and-hold cl' cuit takes its analog sample at the end 
of the multiplexed interval, with the voltage closest 


to its final value. The encoder then digitizes that 
voltage while the next data channel is switching to 
the multiplexer output and is rising to its final 
accurate value. 

The critical step in any analog -to -digital data 
system is the encoding of analog data. Most 
encoders use the binary system; their output is the 
digital equivalent of the analog voltage expressed 
as a binary number, either as time serial or time 
parallel pulses or voltage levels, A majority of 
the existing analog-to-dlgital encoders are closed- 
loop analog-to-digital servo systems. 

The A/D encoders can be divided into six general 
classifications (l^). 

(1) Space encoders. 

(2) Chronometric encoders. 

(3) Successive -approximation encoders. 

(4) Successive-comparison encoders. 

(5) Time parallel encoders, 

(6) Combination or special technique encoders. 

For vibration control systems, successive- 
approximation encoders are satisfactory. This 
encoding technique is widely used today and is the 
fastest conventional time serial logic known. More 
encoders of this type are probably used in data 
acquisition Systems than any other type. 

Starting with the most significant bit, the weight of 
each binary bit is fed to a voltage summer and the 
output of this summer is compared against the 
unknown inpiit voltage. If the s\immer output is 
less than the unknown voltage, the bit is left turned 
on; if the output is greater than the unknown voltage, 
the bit is turned off. This continues until all of the 
reference bits have been "stacked" against the 
unknown voltage. At the end of the series of bit 
weight tries, a voltage s\im is attained which equals 
the unknown input analog voltage with an accuracy 
of ±1/2 the least significant bit weight. 

Most successive-approximation encoders are 
closed-loop analog-digital servos consisting of 
three main elements: a D/A converter (either 
resistive ladder, capacitive or inductive); an 
amplifier to detect the difference between input 
analog voltage and the converter output; and a 
digital servo. 

Resistive and inductive D/A converter ladders are 
capable of at least 0. 01% absolute accuracy and up 
to 0.001% relative accuracy. Capacitive D/A con- 
verters yielded 0.05% accuracies. 

The logic of the current summation successive- 
approximation encoder is id''ntical to that of the 
voltage summation successive-approximation 
encoder, but it has the advantage that different 
input voltage ranges are handled simply by changing 
the value of the range resistor. Full-scale voltages 
from 100 mV to 1000 V have been used in actual 
hardware. 
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In summary, then, the purpose of the sample -and- 
hold circuit is to sample the analog signal from the 
filter and hold the amplitude value of this signal 
until the A/D converter encodes this analog ampli- 
tude voltage to a digital value. Both units are 
usually controlled by the same clock, but the clock 
pulse may be delayed in such a manner that the 
S/H samples the analog signal just prior to digital 
conversion. Once the conversion has been accom- 
plished, another sample is taken and the process 
starts over. Naturally, after a conversion, the 
resulting data word in digital format must be tem- 
porarily stored for processing. The rate at which 
samples of the analog signal are taken Is called 
the sample rate (or sample frequency). This rate, 
as mentioned above, must be greater than twice 
the frequency of the cutoff frequency of the low- 
pass filter, and preferably four times that fre- 
quency, For example, if a random test is to be 
conducted and the test specification requires spec- 
trum control to 2000 Hz, the sampling frequency 
should be 8000 Hz. This sampling frequency would 
automatically be selected by the control system 
computer once the upper breakpoint frequency was 
specified during the conversational mode of the 
control program. In addition, the low-pass filter 
would also be selected for a 0 dB cutoff frequency 
of 2000 Hz. 

5. Sampling Theory . Dynamic continuously 
varying signals from the filter have a high proba- 
bility of changing in a manner that is largely 
unpredictable. The greater the unpredictability 
over a given period of time, the greater is the 
amount of information that may potentially be 
received in that time. 

On might think that, if we recover information that 
determines what was happening all the time, then 
we must have made observations all the time. This 
is not true. Any physical system has built-in 
behavioral limitations that determine how fast its 
output signal can change. In mechanical systems 
this is related to inertia. In electrical systems it 
is related to bandwidth. 

It is often thought that the highest sampling speed 
is needed during the period when the input infor- 
mation IS changing most rapidly. But this is not 
so. When something is moving rapidly, it is likely 
to continue to move in the same direction. Only 
when an event ceases does it become impossible to 
predict what is going to happen next under the 
influence of random external events. 

Consider Fig, 21, which shows a time-varying 
signal being measured in a time AT at time t. 
Suppose we wish to know the amplitude at time t. 
The ability to give an accurate answer Is based on 
two factors. One Is the ability to determine time t 
accurately; the other is the ability to accurately 
measure the amplitude at the precise time t. The 
true point sample usually referred to in most 
mathematical analyses cannot be achieved. Rather, 
there will always be some uncertainty about the 
exact signal amplitude at the instant of sampling 
and about the exact time of the measurement. 

There will always be some jitter in a timing circuit 
and somo finite time required to make a 
measurement. 


Quite clearly, if the signal changes by an amount 
AA during the time AT, and if AT represents the 
uncertainty of the time, then there is an uncer- 
tainty of AA in the single -point sample. The effect 
of variation of AT is not the same for all kinds of 
analog-to-dlgital converters. The effect for a 
successive- approximation converter without 
sample -and- hold differs from that with sample - 
and-hold. The effect differs yet again if one uses 
a ramp converter or a dual- slope converter. 

Now, if we make a single -point sample on a time- 
varying signal, the only information that we obtain 
is an approximation to the amplitude at that partic- 
ular time. We do not recover the entire signal. A 
little intuition reveals that, if we were trying to 
recover the entire time-varying signal over a 
period of time, the overall effect of finite measure- 
ment time in a series of measurements would be 
different from that of a single -point sample. For 
example, if the signal amplitude were continually 
changing, part of the time it would be increasing 
upward and part of the time it would be decreasing 
downward. Analysis of values resulting from a 
succession of spaced measurements of finite 
duration will yield more information than a single 
measurement made at that particular time. This 
is because the error can be smoothed out. 

When an analog signal is digitized it is "quantize^." 
There are only a finite number of possible output 
codes that any digitizer can yield. An input analog 
signal may pass through all possible contiguous 
levels with almost infinite resolution, but the out- 
put digital code must "jump” from one value to the 
next. This is illustrated in Fig. 22, in which it is 
assumed that when a variable reaches the halfway 
mark between successive quantization levels the 
output jumps to the next quantization level. For 
this condition, assuming an otherwise perfect 
digitizer, there are points for which the error will 
be zero. But, most of the time, the resultant code 
will be either higher or lower than the actual value. 


The actual error function for a llnearily changing 
signal is shown at the bottom of Fig. 22. As can 
be seen, it is a sawtooth signal With a peak ampli- 
tude of plus and minus half the quantization level. 
The rms error for a sawtooth is the peak amplitude 
divided by so the rms quantization error is 
Q/(2 n/ 3), Thus when a signal if. digitized, 
’’quantization noise” is added to the signal. The 
relative amplitude of this noise depends on the 
resolution of the digitizer; its component fre- 
quencies depend on the signal form. 


If the original signal is later to be recovered, 
quantization noise must be taken into account. 
Sometimes one needs to digitize signals that have 
an extremely wide dynamic range. For example, 
in a high-fidelity music passage, the dynamic 
range, from a quiet tinkling sound to a loud boom 
of percussion inst»*ument8, might be 80 dB — or a 
range of 10, 000 to 1. For this example, if the 
digitizing device did not have a resolution of at 
least 10, 000 to 1, the quiet passages would be 
missed: or alternatively, the very loud passages 
would be clipped or distorted. 
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Figure 23(a) showa the amplitude waveform of a 
continuously varying signal. If the signal is derived 
from a real physical system, there are limitations 
imposed on its maximum rates of change: its first, 
second, and third derivatives, etc. Because real 
physical systems have some mass (or equivalent 
parameter), they are not capable of having such 
large -amplitude excursions at high rates as at 
lower rates. 

The characteristics of a signal over a period of 
time can also be described by its spectral distri- 
bution, as shown in Fig. 23(b). The Fourier 
relationships tell us that, if we know the amplitude 
function of time, we also know the spectral distri- 
bution for the signal during that time. Conversely, 
if we know the spectral distribution and phase 
relationship between these components for that 
period of time we also know, uniquely, the ampli- 
tude function versus time. A knowledge of one 
yields a knowledge of the other. 

The spectral distribution will normally have a 
relatively high amplitude at low frequency and 
lower amplitudes at higher frequency, eventually 
tailing off to negligible amplitudes. The maximum 
possible rate of change for a particular parameter 
may be approximated by summing the maximum 
rates of change for all possible signal components 
of different frequencies (assuming that these signal 
components are in phase). 

For random noise, shown as a function of time in 
Fig. 23(c), the signal may fluctuate continuously. 
For wideband white noise, the spectral distribution 
extends over an "infinite” bandwidth, and thus the 
rate of change of signal level may approach infinity. 
In real systems, of course, there are physical 
limitations on the amplitudes of the noise (deter- 
mined by the noise pov;er) and on the possible rates 
of change (dependent on actual noise bandwidth). 

In Fig, 23(d), the spectral distribution of wideband 
white noise is shown as being flat. That is, there 
is an equal probability of energy distribution over 
each frequency interval. In most unfiltered sys- 
tems the noise bandwidth, regardless of the level 
of noise power, will be wider than the signal band- 
width. This means that the noise spectrum will 
extend out further in frequency, beyond that 
expected from the actual signal spectrum. 

Figure 24 shows a varying signal being sampled in 
two different ways — at closely spaced points and 
at not so closely spaced points. If we were to 
reproduce just the closely spaced points on another 
sheet of paper, and the limitations on rate of 
change of the signal were known, then (without 
reference to Fig. 24) the points could be inter- 
connected to yield a close approximation to the 
signal indicated in Fig, 24 by a solid line. If, 
however, the only points taken were those inter- 
connected by the dotted line, and if just these 
points were reproduced on another sheet of paper, 
then an attempt to interconnect the points using the 
same imposed limitations of rate of change wotUd 
result in the interconnection shown by the dotted 
line. 


The function indicated by the dotted line appears to 
be changing at a slower rate than the signal repre- 
sented by the solid line. The frequency of change 
implied by the dotted line is called an "alias" of 
the frequency implied by the solid line. 

Thus, if signals are not sampled frequently enough, 
false information consisting of aliases may result. 
The low- pass anti-aliasing filter and four times 
cutoff frequency sample rates reduce this false 
information (aliasing). 

When a signal is sampled (see Fig. 25), the time- 
varying function f^(t) can be considered to be 
multiplied by a sampling function fs(t). During the 
sampling time, the sampling function is considered 
to be of unit height. That is, at that time the 
varying function is multiplied by a function of zero 
height. The sampling fxinction fg(t) then is a series 
of equally spaced flat-topped pulses having a time 
duration of AT and a repetition rate of fs = 1/Tg. 

The resultant observed signal, then, will be that 
shown at the bottom of Fig. 25 at which only 
regularly spaced portions of the original time- 
varying function have been observed. 

For any particular sampling function, as in Fig. 26, 
there exists a corresponding spectrum that is pre- 
cisely calculable. As is well known, flat-topped 
pulses of constant repetition rate are represented 
by a series of lines in a spectrum. The line ampli- 
tudes are determined by the (sin x)/x function and 
their spacing is fg. Thus, there is a dc component 
in the spectrum cf the sampling signal equal to the 
average level of signal, which is clearly 1 x AT/Tg. 
Lines of lower amplitude occur at fg, 2fg, 3fg, 

4fg, etc. 

If it is understood that the original signal is multi- 
plied by the sampling signal having an amplitude of 
0 to 1, then a knowledge of transform relationships 
makes it clear that the original spectrum is being 
multiplied by the spectrxim of the sampling 
spectrum. 

In a digitizing-type sampling system, the resultant 
multiplied product is not quite as shown in Fig. 25. 
Only a single number or level is obtained for each 
sample, not a varying section. This single number 
obviously contains less information than an analog 
sample that indicates the variation during the time 
AT. This is because, unlike the digitized sample, 
the finite -time analog sample is capable of yielding 
trend or derivative information. Also, the single 
value obtained by digitizing falsely implies that a 
true single-point sample in zero time was made. 

But this is not what takes place physically. It has 
been shown that the action is approximately equiv- 
alent to having passed a single-point sample 
through a network that spreads the value over the 
time AT. As we know, the network that converts a 
single -point sample, or impulse function, to a 
pulse of length AT has the transfer function 
(sin x)/x, where the first null is at frequency 1/AT. 

The effect of not making point samples in zero time 
thus yields a result similar to having passed the 
original signal through a low-pass filter with a 
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(•in x)/x chArRcterifltic. Therefore, when aignali 
having high bandwidth are to be sampled, one must 
ensure that sampling interval is not too long to 
provide the desired accuracy. It is for this reason 
that sample-and-hold devices with short aperture 
times are frequently used prior to digitiaing, as 
in Fig. 16. 

We've alrep ly seen that if a signal (as shown in 
Fig« 26) i I rnultipUed by a sampling signal (as 
show > in Fi;. 25), each component of the sampling 
spec iun ' ill be multiplied by the original 
spectrum The original spectrum will be mxilti- 
plievl by a dc component, thus yielding itself; and, 
in addition, it will be multiplied by fs to produce an 
additi')ial spectrum symmetrically around fa. This 
occurs because, when two signals having different 
frequencies are multiplied, the product contains 
all the sums and differences of the two frequencies, 
Similari/, here will be symmetrical spectra 
around the larmonic frequencies Zfg, 3fa« 4fa, etc. 

If the or'g* al spectrum is characterir^ ^ by that of 
an "ideal" ilter, then, as shown in i 27, the 
spectra fre.n sampling multiplication ..ill be rec* 
tanguXar, The first one, starting at zero frequency, 
has .1 width of f^o, while the other spectra sur- 
rounding 2fa, 3fs, etc, , will each have a width of 
2 fco’ 

After sampling, then, the spectrum of the resulting 
signal is as shown in Fig, 27. If, however, we 
wish to recover the original signal, we must elim- 
inate the spectra centered around fs, 2fg, etc. 

This is achieved by filtering out frequencies other 
than the low-band spectrum which represents the 
original signal. 

Thus, i an ideal low-pass filter were employed 
wiuch p ssed the low-band spectrum but did not 
pass the spectra surrounding fs, 2fs, etc. , the 
output would be the origins' signal, unmodified. 

Suppose that our ideal low-pass filter had a band- 
width of just fco and thus was perfectly flat over 
the range from zero to fco* attenuating all other 
frequencies at an infinite rate. It would not modify 
the original low-band spectrum, and it would 
elin'iinate all components of the sampling spectra 
around f^, 2f^,, etc. 

For this tinreali itic supposition, we can immedi- 
ately^ see the lowest mathc Tiatlcally possible 
adequate frequercy of sam^’\ig. We can move fg 
downward in fre<mency until f^ - fco just equals 
fco* Then, solving fo^ co* It is quite clear that 
2fco = fs* if fs we I made lower than 2fco* the 
low-frequency end of the spectrum surrounding fg, 
which contains .e high-frequency end of the 
original spect im, would "fold over" into the 
original spe*:trum and the folded signals could not 
be separateu from or distinguished from the origi- 
nal sign , Thus, fg/2 is considered to be the 
lowest .equency at which a signal having frequency 
con onents up io f^^ can b^ sampled and still 
theoretically be recoverr d without inaccuracies. 
This is the theorem y'elding Nyquist's result; i, e, , 
equlspaced data, with two or more points per cycle 
of highest freq u nc/, allows reconstruction of 
band-limited 1 metions. 


However, real signals and their spectra are not 
characterised by ideal low-pass functions; nor are 
ideal filters available for recovery separation as 
previously mentioned. 

Figure 28 represents more realistically what may 
happen. The original signal spectrum tends to 
have a decreasing amount of energy at upper 
frequencies. Superimposed on the signal spectrum 
will probably be a noise spectrum that tails off at 
a much slower rate. When the combined signal- 
plus-noise spectrum is multiplied by a c^impUng 
signal rate fg, there results a new spectrum which 
includes the original spectrum plus spectra sym- 
metrically located about fg, 2fg, etc, 

Thus, for a real signal characterized by a spec- 
trum decreasing at a finite rate and never reaching 
zero amplitude, there is no real value of fg equal 
to twice the highest irequency present. Further, 
there is no real filter available that can cut off at 
an infinite rate and separate the spectra. 

Therefore, to determine the minimum allowable 
sampling rate in a practical situation we must 
carefully analyze the original spectrum and its 
noise content. And we must use our judgment to 
determine what degree of error and distortion can 
be tolerated. To decide on a lower limit for fg, 
the test engineer must know the quality of the 
filter that can be used to separate the original 
spectrum from the higher-order spectra. He musl 
make a judgment not as to whether there will be 
any signal folded into the original spectrum but as 
to how much he can tolerate — and not whether he 
can allow any noise to be folded into the original 
spectrum, but how much he can tolerate. 

The recovered signal must be obtained with a 
minimum time delay so that the sampled data can 
be used in a real-time control loop. This imposes 
further restrictions on the nature of the recovery 
filter. As mentioned earlier, the more complex a 
physical filter the greater will be its time delay. 

(A similar argument applies to digital computation.) 

Thus, quite clearly, a sampling system must be 
considered to be just that: a system. The overall 
transfer function must be considered if any realis- 
tic estimate of the accuracies to be achieved are to 
be obtained. The test engineer must take into 
account the noise level, the noise bandwidth, and 
the shape of the original spectrum. He must 
determine what amount of frequency folding he can 
tolerate and, then, by considering all of the factors 
combined, he can determine what minimum rate of 
sampling to use. 

It has been determined experimentally from recent 
computer controlled test systems that the digitizing 
process is adequate when the low-pass antialiasing 
filter has at least 8 poles (-48 dB/octave, 
Butterworth response), the A/D converter is at 
least a 12-bit device, and the sampling frequency 
is four times the 0-dB cutoff frequency of the filter. 
As mentioned above, one truly never gets rid of the 
aliasing errors but only reouces them to a level not 
worth worrying about for the control application. 
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6. Digital-to^Analog Conversion , Digital -to - 
analog (D/ A) conversion is necessary in the 
computer controlled test system in order to 
implement the following important functions. 

(1) Produce an analog drive signal to the power 
amplifier. 

(2) Produce analog signals to the oscilloscope 
display devices. 

(3) Produce drive signals to the X-Y recorder. 

Item 1 above necessarily includes, besides a good 
quality D/A converter, a programmable amplifier/ 
attenuator and low- pass filter. The filter is 
normally identical to the input anti-aliasing filter 
and is necessary to reconstruct a true analog sig- 
nal from the D/A converter. Outputs from unfil- 
tered D/A converters exhibit signals that are step 
voltages and appear as synthesized waveforms. 

The filter removes the discrete steps or points by 
eliminating the high frequencies associated with 
discrete steps. Figure 29 is a blocV diagram of 
the analog output requirements of a computer- 
controlled induced environmental te^t system. The 
programmable amplifier /attenuator determines 
the amplitude voltage requirement to the power 
amplifier in order to meet the test specification 
across the desired frequency spectrum. The 
actual value of attenuation or amplification deter- 
mines the loop gain as does the input amplifier/ 
attenuator. These gains are determined by the 
system during the conversational language mode 
of the control program and are dependent upon 
transducer sensitivity, internal scale factors, 
and test level. 

7. Trigger Logic. The trigger logic circuit 
(FigT lb) is necessary for the analysis function of 
the computer-controlled test system and is used 
for shock pulse analysis. The purpose of the trig- 
ger logic Is to start the analog-to-digital encoding 
when the test system senses the presence of the 
leading edge of an analog shock pulse signal. This 
provides an automated means of starting the data 
conversion a split second after detection of the 
shock pulse. This guarantees that the shock pulse 
will be captured by the test system. 

In some systems it is possible to trigger on either 
a positive or negative slope and, in addition, the 
operator has a choice as to where (in amplitude) 
on the slope of the leading edge to start the 
analysis. 

8. Signal Conditioning Errors . Table 3 lists 
error factors and potential error factor problems 
associated with the data acquisition and condition- 
ing systt:m illustrated in Fig. 16. This acquisition 
and conditioning system is also necessary for the 
third function that the computer-controlled test 
system serves to the environmental laboratory, 
which IS described next. 

C. Signal Analysis 

The third function that the computer -controlled 
test system provides for the Induced Environmental 
Test Laboratory is the time series analysis 



Table 3. Error factors from input data channel 

Analog signal 

Transducer nonlinearities. 

Programmable amplifier /attenuator 

Input impedance. 

Input current. 

Offset. 

Transfer accuracy. 

Linearity. 

Interference susceptibility. 

Programmable low-pass (anti-aliasing) 
filter 

Amplitude response. 

Phase response. 

Transient response. 

L ine^rity. 

Tiansfer .''ccuracy. 

MuUiplexer (for muii* ole -channel 
control /analysis) 

Cross taik- 
Transfer accuracy. 

Leakage current . 

Sample -and-hold 

Transfer accuracy. 

Droop. 

Offset. 

Aperture time. 

Acquisition time. 

Analog-to-digital converter 

Speed. 

Rate. 

Mode of conversion. 

Resolution. 

Monotonicity. 

Stability. 

Interference susceptibility. 


function, and in particular, the following fre- 
quency functions 

(1) The Fourier transform. 

(2) The powei spectral density function G^(0. 
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(3) The cross spectral density function G (f)« 

xy 

(4) The frequency- response function H(f)» 

Some or all of these functions are required for 
computer -controlled random vibration tests and 
transient shock tests. 

1, Windowing s In digital time series analyzers, 

It Is necessary to window the incoming data from 
the A/D converter prior to calculating some of 
the above frequency functions. This Is due to the 
spectral content of the sampling fxinctlon fg, shown 
in Fig, 26, and, in fact, the ill effects of this 
spectrum past the frequency 1/AT. The sample 
function fg Is also referred to as the window func- 
tion (in the frequency domain) or the '‘boxcar" 
function In the time domaln(lo), xhe raw estimate 
of the above -listed frequency functions is the con- 
volution of the true function with the sample func- 
tion in the frequency domain or their products in 
the time domain. The effect of the sample function 
causes errors in the estimated spectral values of 
these frequency functions. To reduce these 
errors, the digital analyzer portion of the control 
system must employ data weighting. Weighting or 
windowing such as the Hanning function, described 
below and in Fig. 30, smooths the raw data esti- 
mates by modifying the normal "boxcar" time 
domain function and its associated (sin x)/x type 
frequency function. Windowed analyses result in 
better frequency function estimates than non- 
windowed analyses. Many types of weighting func- 
tions are defined, and the choice ultimately is a 
factor in the upper bound of the control systems 
real-time processing capability. The Hanning 
window or weighting can be accomplished in che 
time domain or in the frequency domain after the 
spectral values of the frequency functions have 
been generated. For the latter case, assume an 
unHanned data block of spectral values stored in 
BLKl of elements x(k), k = 1, 2, n. After 

the Hanning weighting, define another data block of 
unnormalized weighted spectral values to be stored 
in BLK2 of elements y(k), k = 1, 2 n. Then 

Yj = I (Xj + X2> (14) 

'fz = *2 +1*3 


A general algorithm for performing this process is 
illustrated in Fig. 30. The above Hanned values 
are "tmnormalized" and must be multiplied by a 
factor of (8/3) assuming the function being 
Hanned is a Fourier transform and the x(k)'s are 
Fourier coefficierts. If the x(k)'s are spectral 
values of the PSD or cross spectrum function, 
then the multiplying factor is 8/3. This will 
become apparent when the PSD function is 
described in this section of the paper. It turns 
out that it is most desirable to window the Fourier 
coefficients provided these coefficients are used 
to generate functions of frequency and not the 
transforms (corresponding time functions) of these 
frequency functions. 

The windowing or weighting described above is a 
process that must be done on the input data from 
the A/D converter, or immediately after the 
Fourier coefficients have been generated in order 
to control a random noise vibration or acoustic 
test. There is another type of so called "window- 
ing" which is performed on the output signal to the 
power amplifier. This windowing is used to con- 
catenate output frames together in the random 
vibration control process. 

2. The Fourier Transform. ^ The Fourier trans- 
form is the basic frequency function that a 
computer-controlled test system or a time series 
analyzer performs. It enables the system to com- 
pute all the other time and frequency functions 
such as correlations, density functions, transfer 
functions, and the coherence function. 

The Fourier transform in essence is a mathemat- 
ical operation performed on a time history func- 
tion to determine and analyze its frequency content. 
The results of the Fourier transform process are 
therefore the magnitudes and the phases, called 
complex Fourier coefficients, that the signal 
possesses at the various frequencies within the 
control (or analysis) bandwidth. In digital systems, 
the input time signal is digitized as mentioned 
previously so that it is represented by a finite 
number N, called the frame size (FS), of ampli- 
tudes. The Fourier transform operation per- 
formed is therefore necessarily a discrete rather 
than a continuous process. The general form of 
the complex Fourier coefficients is: 
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(16) 
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k=0 

(19) 

I = 0, 1. 2 N-l (W = = -1) 


Vl = l*n-2 +i Vl +5*n 

= I‘*n-1 +*n> 

^Subsection 2 w.s written by W. Boctor. 


(17) 

where N is the frame size and xj^ is the kth digi- 
tized time sample. This method of calculating 
(18'i Ihe Fourier coefficients involves N complex oper- 
ations (i.e., complex additions and complex multi- 
plications) for each coefficient, i.e., 
complex operations in all. 


16 



CoQfliderable savings in time can be achieved by 
reducing this number from to N log^N, This is 
realized by performing the Fourier transform in 
steps* following the Cooley* Tukey( ^7) algorithm. 

If m steps are involved, then the number of oper- 
ations become N (ri + r2 I- . . . + *vich that 

(rj • r 2 • . . . • rm) = N, and r^ (i = 1, 2 m) 

are the radixes used in the different steps. 

Although the minimum number of operations is 
achieved by choosing all r^ = 3, rj = 2 proves to 
be much more advantageous for use by digital 
computers because of their binary nature. The 
loss in speed over r^ = 3 is only about 6%, This 
yields what is now known as the radix-2 algorithm 
and the Fourier transform thus calciUated is called 
the fast Fourier transform (FFT). This method 
provides the relationship 

N = 2*” 


for the inverse Fourier transform (IFT) also, 
which is used to compute the correlation functions 
from the spectral density functions in analysis, 
and more important, it is used in control to com- 
pute the signal to be sent out to the shaker or other 
exciter from the drive spectrum. The factor ISIGN 
is - I for an FFT and +1 for an IFT, 

It is worth noting that the cosines and sines used 
in the algorithm are generated in non- real-time 
and saved as a lookup table, and depending upon 
the frame size N, only every value 

is used, where maximum system 

frame size. Usually sines and cosines of angles 
within the range 0 to 90 deg only are generated 
because of the large memory storage required (in 
software FFT) or the large ROM (read only mem- 
ory) size required (in hardware FFT). 
Mathematically, 

cos (ISIGN, M, L) = cos (ISIGN ^ Ztt ^ (M-l)/L) 


and number of operations = N log^N, 

This sets the imperative requirement that FS be 
an integral power of 2, A flowchart that tests for 
this relationship is shown in Fig. 31. 

In this algorithm, the res\ilts obtained In the (n-l)th 
stage are utilized in calculating 2"^”*' sums of the 
form 


s. 
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(W„ = 
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( 20 ) 


This goes on until all m stages are performed, the 
results of each stage depending upon and spacially 
replacing the results of the previous stage. The 
FFT thus calculated is called an in-place FFT 
because the output block containing the real and 
imaginary pirts of the Fourier coefficients is 
coincident with the input block which previously 
contained the input time signal. In this case of a 
real input, the output block is skew- symmetric 
about the dc value, and only the positive frequency 
components need be calculated and retained. Thus 
an N-size real-time input block yields N/2 real 
frequency components and N/2 imaginary frequency 
components, and the FFT can still be performed 
in place. However, this requires the input time 
block to \indergo first a shuffling process so that no 
compon<>nt8 within it get replaced in any stage by 
others (resulting from computations in that stage), 
except when they are not needed any further. 

Fig, 32 shows an example for a flowchart of the 
shuffling process. 

An example of the main FFT flowchart is given in 
Fig. 33, It is seen that the main algorithm is used 
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sindSIGN, M. L) = sin (ISIGN Zit * (M-l)/L) 

( 22 ) 


The sine-cosine table, whether implemented in 
software or hardware, is usually in integer form to 
speed up the multiplication operations and is such 
that a sine or cosine value of 0 is represented by J 
whereas a sine or cosine value of 1 is represented 
by the largest integer number the computer regis- 
ters can accommodate (32767 in a 16-bit machine). 
This offers two great advantages; (1) it achieves 
the highest resolution possible with the existing 
machine, and (2) when multiplying any number (in 
the multiplicand register) by a value from the table 
after moving it to the multiplier register (again 
whether software or hardware), the product 
obtained in both registers may be chopped off with 
the original register (the multiplicand register) 
containing the correct 16 -bit value of the product. 
No further shifting or normalization is required. 
The reason here is that with the sine-cosine table 
in this form, all its values are inherently multi- 
plied by 2^“l, where n is the register si? in bits 
(i. e. , multiplied by 32767 in this example). And 
chopping off the right half of the product is equiva- 
lent to a division by that same factor, so that the 
end result in the multiplicand register is the cor- 
rect one. 


The FFT/IFT algorithm described above is suitable 
for both software and hardware implementation. 

The former is leas expensive than the latter but it 
is considerably slower (about 80 to 100 times on the 
average, as shown in Table 4). Hardware FFT/IFT 
is therefore more suitable in digitally controlled 
testing. In this case the sine-cosine table is a 
read-only-memory containing all fhe sine and 
cosine values within the first quadrant for fast 
access. The complex operations are all performed 
in integer form in hardware registers. The fast 
Fourier processor, as it is called, transfers the 
input block from the computer main storage to a 
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Table 4. Comparlaon of execution timet (in mt) 
of forward and inverse Fourier transforms in 
both software and hardware implementation (hard- 
ware data are according to highest speeds com- 
mercially available in January 1974) . 


Frame 

size 


Execution time, ms 


Software 


Hardware 


FFT 


IFT 


FFT 


IFT 


64 

32 

38 

0.58 

0. 58 

128 

75 

90 

1.09 

1,09 

256 

170 

200 

2.07 

2.07 

512 

380 

450 

4.9 

4.9 

1024 

850 

1000 

9.8 

9.8 

2048 

1650 

2100 

21.7 

21.7 

4096 

3200 

4300 

37.8 

37.8 


buffer block within it, performs the FFT or IFT 
on it, and puts the result back in the original 
main storage block in the computer memory, 

3. The Power Spectral Density Function . The 
power spectral oensity function (PSD) G (f) is 
also known as the normalized power spectrum, 
the normalized auto spectrum, and the auto 
spectral density function. In this paper, the 
function G (f) will be called the power spectral 
density function or simply PSD function. 

The G^(f) is defined as the self-conjugate product 
of the real and imaginary coefficients of the 
Fourier transform of the time function x(t). That 
is, 


G (f.) = 
x' r 


X (f.) X 


(fj) 


where 


X (f,) = a. + jb. 


and means complex conjugate. So 


G (f.) = (a. + jb.) (a. - jb.) 

X I ' i •' i' ' i i' 


(23) 


(24) 


(25) 


giving 


Gx(fp = a.^ + b^^ at frequency f^ (26) 


The coefficients a and b are the real and imaginary 
coefficients at the frequency f^ of the Fourier 
transform of x(t). The j operator drops out in the 


multiplication process (j^ « -1). What this means 
physically is that there is no phase information in 
the PSD function. 

For vibration analysis and control, it should be 
understood that the Fourier coefficients, as they 
are generated within the system processor, have 
units of peak voltage, and since the coefficients 
are derived from sine and cosine tables, these 
coefficients should be multiplied by a scale factor 
of 1/^/T in order to obtain rms amplitudes. 

If the scale factor is applied after Gx(^i) gener- 
ated, then the scale factor becomes 0* 5 times 
Gx(it)« In addition, if Hanning weighting is 
accomplished after Gx(fj) is generated, the overall 
scale factor becomes 1,^. 

If the signal being analyzed is random in nature, 
the PSD function can be normalized by dividing 
the results by the effective filter bandwidth used 
in the analysis to provide a characterization In 
terms of 1-Hz bandwidth. This allows direct 
comparisons of analyses to be made, T PSD, 
as defined here, is th 'formalized function that 

^he PSD of a random 
me strictest sense, and 
\>L i always implies an 


has rms units of g^/. 
variable is an estima 
the term PSD in this 
estimate. 


Figu.i ' 34 illustrates the basic process required to 
implement the PSD function using a digital machine 
designed around fast Fourier transform concepts. 

A PSD analysis is dependent upon the following 
parameters associated with all digital analyzers. 

(1) Frantie size (FS) and spectral lines (SL); A 
Fourier analysis requires li real values in 
the time domain in order to generate n/2 
spectral values (spectral lines) in the fre- 
quency domain. The frame size is the 
ntimber of data words in the time domain to 
be processed at one time. The spectral 
lines resulting in a PSD analysis is one-half 
the frame size number selected for the 
analysis. 

(2) Sample rate (SR), Sample period (SP) : The 
sample rate is tne points per second selected 
by the operator for the analog -to -digital 
conversion associated with digital ana- 
lyzers. The units of sample rate are points 
per second, samples per second, or con- 
versions per second. The sample period 

is the reciprocal of the sample rate. 

(3) Mode (M) : The mode, continuous (real- 
time) or discontinuous (non-real-time), is 
a function of the sample period. If the 
sample period SP is too short (i.e. , the 
sample rate too fast) for the number of 
spectral lines selected, the analyzer will 
Indicate and automatically go into a dis- 
continuous mode of operation. 

(4) Analysis bandwidth (AB) : The analysis 
bandwidth is defined as some fraction of the 
sample rate; the units are hertz. The analy- 
sis bandwidth must satisfy the minimum 
requirements of the sampling theorem 
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AB < 0. 5 SR H* (27) 


or 

AB<5|sHr (28) 


(5) Effective filter bandwidth (FB) ; The effec- 
tive filter bandwidth is defined as the 
analysis bandwidth divided by the spectral 
lines 


FB = AB/SL = SR/FSHz (29) 


(6) Load time (LT); The load time is defined 
as the time required by the digital analyzer 
to obtain one frame of data commensurate 
with the selected frame size. The load 
time is the indicated frame size divided 
by the sample rate 


LT = FS/SR (30) 

or 

LT = (FS) (SP) (31) 


The load time, digitally speaking, is 
equivalent to averaging time in the analog 
domain. The effective filter bandwidth 
and load time product is always one. 


(FB) • (LT) = 1 (32) 


(7) Analysis time (AT); The analysis time is 
the actual time required to calculate a total 
estimate of K statistical degrees of free- 
dom from time frames of input data. 


AT = K • LT/2 (33) 


(8) Statistical degrees of freedom (K) ; The 

number of degrees of freedom of a spectral 
estimate is a convenient description of the 
statistical reliability. The number of 
degrees of freedom K per frame is defined 
as twice the effective filter bandwidth 
times the load time. 


K/frame = 2 (FB) (LT) = 2 (34) 


For all digital analyzers, a single spectral 
estimate of two degrees of freedom is pro- 
duced from each time frame of Gaussian 
distributed data since the effective filter 
bandwidth and the load time product is 
always 1, 


As more time frames are analyzed, the 
degrees of freedom go up, indicating better 
reliability of the spectral estimate. It is 
important, therefore, to express a quality 
factor of degrees of freedom per second of 
analysis time. 


K/sec = 2/LT for continuous mode 

and 

K/scc = 2/(LT + FS • SP) for 
discontinuous mode 


The above parameters are all in effect during a 
PSD analysis. Some of the parameters must be 
selected prior to the process implementation illus- 
trated in Fig. 34. These parameters are gen- 
erally the analysis bandwidth and effective filter 
bandwidth as well as the number of averages 
required in order to obtain acceptable statistical 
reliability. The selection of these parameters 
sets the remaining parameters such as frame size 
and sample rate. Load time is a physical result 
from the previous choices as is the statistical 
reliability once the number of averages have been 
selected. The analysis time is a function of the 
algorithm, hardware and software. The param- 
eters are selected or calculated during the con- 
versational mode of the control program. 

The area under the PSD locus is the mean- square 
value of x(t): the rms of x(t) is the square root of 
the mean- square value. The digital analyzer must 
perform a numerical integration. This can be 
accomplished by implementing the trapezoidal rule 
or Simpson's rule. For vibration and acoustic 
analysis the trapezoidal method of integration is 
generally adequate and easy to implement by divid- 
ing the first and last spectral densities by 2 and 
adding these values to the sum of all the other 
spectral densities. The final sum is then multi- 
plied by the effective filter bandwidth to obtain the 
mean-square value of x(t). A square root algor- 
ithm is then required to obtain the rms value 
of x(t). 

4. The Cross Spectral Density Function . The cross 
spectrum function Gxy(^ is also called the cross 
power spectrum or, when normalized to I Hz, cross 
spectral density function. The cross spectrum func- 
tion relates to two time signals, both of which could 
appear as either forcing functions or response func- 
tions but usually an input forcing function x(t) and a 
response function y(t). In any case, the Fourier 
transforms X(f) and Y(f) are required to obtain 

Gv„(f). where, at frequency f,, 
xy 1 

X(f.) = aj+jb. (35) 

and 

Y(f.) = c. +jdj (36) 
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A cross ^conjugate product is performed in the From £qs. (39)» (40) and (42), 

computer control system as follows:^ 


or 

°xy<V = Vi + Vl + j< Vi - Vl> 

where ai, bj, c^, and d: are the Fourier coeffi- 
cients at frequency f^ of the time functions x(t) and 
y(t), respectively. 

It can be seen from £q. (38) that the cross spec- 
trum is a complex frequency function; that is, it 
has a real and imaginary part, unlike the power 
spectral density function Gj^(f). Let 


mag lG^y(fj)| = [(..Cj+b^d^)2 

1/2 

+ Kdl-b^Cj)2] (44) 

r 

MAO |0„«,ll . (Vl)' 


Pi = ajCj+b.dj 


and 


(39) 


1/2 

* (Vi)"] 


Q: 


a.d. - b.c. 

II 11 


(40) 


rearranging and cancelling some terms. 


Then 

G^yUl) = Pi+jQi (41) 

where P is called the coincident term or simply co 
term and Q is called the quadrature term or quad 
term. For rms amplitude units, Eq. (7) must be 
multiplied by a scale factor of 0, 5, and a scale 
factor of 8/3 for Hanning weighting, for a total 
scale factor of l.?3. 


MAO |0^(f,l| = (i,^ + V) 

1/2 


or 


Equation (41) is the rectangular form of the cross 
spectrum function. The function should be normal- 
ized to 1 Hz by dividing by the effective filter 
bandwidth. The normalized polar form is 


v'l' ■ K" * °i") 

(42) 


1/2 

MAO |Q^(Ii)| = [(.,^b;2)(ci^di2)] 

(47) 


From the definition of the power spectral density 
function (PSD), 


where 





TAN"* Q./P. rad 


(43) 


It can be seen that the cross spectrum function 
Gxy(f) yields both amplitude and phase information. 


MAG iG^ydi)! 


1/2 


(48) 


Therefore, either Eq. (42) or Eq, (48) can be used 
to calculate the normalized magnitude of the cross 
spectrum function from the digital analyzer. 


Gxy(A) could be calculated by the cross conjugate product of (a^‘fjb|)(ci-jd^). However, this product will 
result in a phase angle 180 deg apart from the frequency response function H(f), defined next in this 
paper. 
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The cross spectrum function by itself is not widely 
used. It is used to calculate the frequency response 
function described next in this paper. 

The cross spectrum function and the cross corre- 
lation function are transform pairs in the frequency 
and time domain, respectively. The power spec- 
tral density functions G^Cf) (input forcing function), 
Gy(f) (output response) oi x(t) and y(t), respec- 
tively, and the cross spectrtim of x(t), y(t), i. e, , 
Gxy(f)» will yield the total frequency information 
required for an understanding of the physical pro- 
cess resulting from the excitation of a structure, 
provided that the cross spectrum is measured at 
both transducer points simultaneously and both 
PSD functions are derived from that same data. 

Figure 35 shows the processes required to obtain 
the cross spectrum function. All of the param- 
eters defined for the PSD function are in effect for 
the cross spectrum function. These parameters 
include: 

(1) Frame size or number of spectral lines 
desired for the analysis. 

(2) Sample rate or sample period. 

(3) Analysis bandwidth. 

(4) Effective filter bandwidth, 

(5) Load time. 

(6) Analysis time. 


function H(f) goes by several names. It is also 
called the transfer function and the system func- 
tion, but strictly speaking, H(f) is the frequency 
response function, and the term transfer function 
should be reserved for the Laplace transforn^ of 
the impulse function h(t), l.e. , H(t). 

The frequency response function at frequency f. is 
defined as 


Y(f.) 

H(f.) = ^ (49) 


where X(f[) and Y(f^) are the Fourier transforms at 
frequency f^ of x(t) and y(t), respectively, in 
Fig. 36B, and 


X(f.) = a. +jb. (50) 

and 

Y(f.) = c. +jd. (51) 


where a^, b^, c: and d^ are the real and imaginary 
coefficients of tne Fourier transforms at frequency 
f^, and j is the imaginary operator 

Therefore, 


(7) Statistical degrees of freedom. 

For strictly analysis purposes, a two-channel 
analyzer is required with dual analog-to>digital 
converters. These converters must share a com- 
mon clock such that the conversions associated 
with each converter are simultaneous. Sets of 
Fourier coefficients are generated from the 
Fourier transforms of x(t), the input forcing func- 
tion, and y(t), the response function. A cross - 
conjugate multiply process is performed on these 
Fourier coefficients at each center frequency i of 
the effective filter bandwidth across the desired 
analysis bandwidth. For control purposes, it may 
be desirable to take the cross spectrum of certain 
signals associated with the control process. This 
will be discussed in a later section of this paper. 

5. The Frequency Response Function. Consider 
Fig. i6A, The reaaon x(t) is not identical to y(t) 
is that the test item itself is characterized physi- 
cally by its impulse response function h(t).^ The 
impulse response of the item under test can be 
cot .ndered as a mathematical model for the test 
item in the time domain, since it can be used to 
relate the input x(t) and the output y(t). In the 
frequency domain, the impulse response function 
h(t) transform.s to the frequency response function 
H(f), as seen in Fig. 36B. That is, the Fourier 
transform of h(t) is H(f). The frequency response 


H(f.) 


*l " j**! 


In polar form, Eq. (52) become. 

C, 


|H(£.)| = 




where 


and 


C. = 

i 


A = 
i 


The phase angle 






1/2 


1/2 




(52) 


i53) 


(54) 


( 55 ) 


7 


h(t) is also known as the 


weighting function of a constant parameter linear system. 



T 
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can be calculated by multiplying £q, (52) by the 
complex conjugate of a^ ^ jb^, i,e« » a^ • jb^ to 
obtain 




= TAN 




rad (56) 


This phase angle is identical to the phase angle of 
the cross spectrum function at all frequencies 
across the analysis bandwidth of the frequency 
response function. 

Equations (54) and (55) are the square roots of the 
power spectral density functions Gy(f^) and Gx(fi)» 
respectively (refer to the power spectral density 
function description in this paper). That is, 

Cj = (Gy(£.)) (57) 

and 


Therefore, 



From £q, (48), it is noted that 

1/2 


Dividing both sides of £q. (60) by the PSD of the 
input function x(t), i,e., Gj^(f|), the results are 


MAGlG (fj)! /Gy<^i>\ 

“ I VV ) 


(61) 


The right-hand side of Eq. (61) is precisely the 
magnitude of the frequency response function 
(£q, 59), That is. 


spectrum phase angle. The magnitude of H(f) is 
dimensionless. No scale factors are necessary 
for H(f), but scale factors are necessary for the 
original Fourier coefficients. 

Since both the PSD function and the cross spectrum 
function are estimates in the strictest sense, the 
frequency response function is also an estimate 
as defined and implemented here, and all of the 
previous frequency function constraints pertain to 
this function when a random noise signal is 
employed as the input forcing function to a Dnear 
system whose frequency response is to be ana- 
lyzed, Figure 37 describes the frequency response 
functions, utilizing the PSD and cross spectrur.i 
functions. 

The magnitude of the frequency response signal is 
sometimes called the transmissibiilt :.tion 
Tyjj(f), That is. 



The frequency response function is a very power- 
ful function which provides tremendous insight 
into the dynamics of mechanical systems. This 
function can be utilized in the random excitation 
control algorithm^^^). 

D. Servo Control of Tests^ 

1, Random Control , The essential feature of the 
automatic servo control system is feedback, 'ihc 
feedback is that property of the system which per- 
mits the output quantity to be compared with the 
desired reference level so that, upon the existence 
of a difference, an error signal arises which acts 
to bring the two into correspondence within the tol- 
erance limit. In a digital random control system, 
a reference power spectral density in g^/Hz is 
given as a fiuiction of frequency. In the case of 
the acoustic tests, the reference sound pressure 
level in dB may be given as a function of frequency. 
Other than the engineering units and data presen- 
tation format, random and acoustic control share 
the same servo concept, and the following control 
algorithm U equally applicable to both types of 
environments. 


mag|g (f.)l 
MAG|H(fi)| = a (ff^ ~~ 

x' i' 

Therefore, either £q, (59) or Eq, (62) can be used 
to calculate the magnitude of the frequency response 
function, and as was previously noted, the phase 
angle of this function is identical to the cross 


^Subsections D and £ were written by B. K, Kim, 


All of the system parameters that are necessary 
for control of a random/acoustic vibration system 
are not known a priori and therefore must be esti- 
mated. The driving and response signals are the 
only directly observable phenomena from which all 
control information must be deduced. The basic 
model for the system, as shown in Fig, 38, 
assumes an approximately linear transfer function 
for the shaker H(w), input driving function, x(t) 
additive system noise n(t), and a resultant output 
y(t). Then, 
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y(t) :: h(t) x(t) + n(t) 


(64) 


where h(t) is the impulse response function and 
indicates convolution product. If we assume x(t) 
and n(t) are uncorrelated, then 


Sy(«) = |H(u))! ^ • s^(ui) + Sjj(ui) ‘ E(u.) (65) 


where 

S^(u>) is the drive auto power spectrum 
Sy(w) is the response auto power spectrum 
S^(w) is the system noise auto power spectrum 

E(w) is the err^'r spectrum in estimation 

Note the conventional PSD is the auto power spec- 
trum (APS) normalized to 1 Hz as previously 
defined in this paper. The control problem is then 
to manipulate the drive x(t) or its APS S^(w) fo 
that the desired output spectrum Sy(w) is 
obtained. 

There are several proven control algorithms 
which differ slightly in strategy. The following 
explains the two strategies. 

a. Auto power spectrum control method . This 
method attcmpis to control such that Sy(u ) 

approaches the reference spectrun R(w) as shown 
in Fig. 38. This was initially proposed by 
Heizman (Ref. 1) and later implemented by 
Sotomayor( 1*^). After computing the response sig- 
nal level via APS routines, the computer compares 
the averaged APS to the reference spectrum and 
modifies the drive signal according to 


C(u) = 




RM- 


1/2 




C^(u) = ^u>) (66) 


where R is the reference spectrum 

S^ j is the averaged APS after frame i-1 

C. is the drive signal Fourier amplitude 
^ spectrum for frame i 


and 



represents the latest estimate of 
the system transfer function. 


Note the above equation involves entirely real num- 
bers only, and the random phase is introduced just 
prior to the inverse Fourier transform stage. This 
method of control has shown m practice that it 
tends to require longer convergence time to a given 
reference spectrum, and also each individual 


control spectrum shows greater statistical scatter 
because of the r- >dom nature of the phase. 

(18) 

b. Complex transfer function control method^ ' . 
This method has the distinct advantage of being 
able to yield an estimate of the transfer function 
H(w), which is computed from the exponentially 
averaged cross spectrum and power spectrum 
as shown In Fig. 39. 


H(w) 


E X (w ) ■ ’ ( w) 
£|X(u)i^ 


(67) 


Then the drive signal Fourier coefficients at 
frame i is 


C (u) = JEM . ^ c („) . (68) 

Hi fu) 


Note in this equation that H and C am complex 
numbers, thus preserving the e. >act phase relation- 
ship. The overall performance of this method 
shows faster convergence and a'higher statistical 
confidence level than the APS method. This is due 
to the tact that the binary psuedorandom-noise- 
generated phase angles are reflected in the trans- 
fer function computation process. 

The digital random control presents a new concept 
in interpreting the control data. This stems from 
a unique digital processing method of time series 
data which is frame-oriented and a discrete 
r'ourier transform technique which has no parallel 
in an analog control system. Owi^g to the statisti- 
cal nature of a randr^.m signal, the confidence level 
of a PSD estimate is a direct function of the number 
of frames averaged and the width of a scatter inter- 
val. Assuming the random noise is Gaussian, the 
averaged PSD tends to a chi-square distrib’^ted 
function. Thus it requires a judicious choice of 
number of frame averages to recognize the true 
wideband PSD le/el. This will depend on the mini- 
mum confidence level d esired, the maximum toler- 
ance band acceptable, and the total duration of the 
test. Methods of averaging can also vary as the 
linear or exponential weighing may be applied. 

2. Sine-Wave Control . This is the second of the 
main categories of control of vibiation tests using 
a digital computer. The digital system uses a 
combination of hardware and software algorithms 
to go through all the steps of a full test. The o^'cr- 
all block diagram has already been described in 
Fig. 2, and several of the functional requirements 
of the system have been given ample consideration 
above. They are: 

(1) The interactive routine to input all the 
required test parameters and frequency 
profiles. 

(2) The generation of the sine wave itself. 
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(3) The signal conditioning. 

( 4 ) The computation of the frequency and ampli- 
tude servo. 

The aspect of the system in this section is con- 
cerned with the servo control concept which forms 
the heart of the control system for sine testing. 

This control is achieved by an algorithm which 
basically sends out a signal at a specified fre- 
quency with a specified amplitude. This signal 
goes through all the conditioning steps and reaches 
the exciter in a smooth analog form. The analog 
response coming back is conditioned and digitized 
and its amplitude characteristics are computed. 

The servo algorithm then compares these ampli- 
tude characteristics at the prescribed frequency 
with the preset reference characteristics as deter- 
mined by the input parameters. The program cal- 
culates the difference between the two levels » thus 
giving an error to be algeb ically added to the 
previous amplitude of the signal sent out to the 
exciter, after being multiplied by a convergence 
factor. 

The foregoing is a brief description of the ampli- 
tude servo in a sine wave test. The other aspect 
of the servo, the frequency servo, has already 
been described under the sine-wave generation 
section of this paper. 

Before proceeding with a detailed description of 
the amplitude servo algorithm, it may be appro- 
priate to consider a brief description of the differ- 
ent types of spectra used here. They are all 
amplitude-vs-frequency spectra, 

a. Preference spectrum . This is the frequency 
profiTe in g’s |rms, peak or peak-to-p^ak, depend- 
ing upon the algorithm) which the test operator 
enters in the tes ^etup parameters. It is the 
spectrum to which the rervo attempts to control 
the response from the external load at all fre- 
quencies within the control bandwidth; i. e. , it is 
the d csired spectrum at the control point. 

b. Control spectrum . This is the frequency 
response of the external load to the driving signal. 

It is the spectrum which the servo attempts to con- 
trol to the reference spectrum by modifying appro- 
priately the driving signal. 

c. Monitor spectrum. This is the frequency 
response of monitoring acceleration signals, as 
selected in the test setup parameters. 

d. Error spectrum . This is the difference 
between the control and reference spectra. It indi- 
cates by how much control is off the required spe- 
c’^'-ation level. It is used to correct the driving 
signal to the external load in an effort to bring the 
control spectrum as close as possible to the refer- 
ence spectrum dynamically in real time. Theo- 
retically, 100% full control would be obtained by 
driving the external load such that a zero error 
spectrum is obtained. But this results in no sub- 
sequent correction and simulates an instantaneous 
open-loop condition. The driving signal thus starts 


deviating from the value required for proper 
control, making the error spectrum nonzero 
again, and the feedback locp is closed again. So 
the error spectrum, as any error signal in a feed- 
back system, is only a virtual zero* very close to 
zero but never being exactly equal to it. 

e. Drive spectrum. This is the frequency spec- 
trum of the signal actually sent to the external 
exciter. It is equal to the previous drive spec- 
trum corrected by the error spectrum after adjust- 
ing it and properly conditioning it, 

f. Description of the amplitude servo algorithm . 
Figure 40 shows a flowchart the steps involved 
in the amplitude servo algorithm. The servo first 
starts at the initial frequency and outputs a sine- 
wave cycle with an amplitude lower than the preset 
reference level at that frequency. It then looks at 
the control signal coming back, computes the error, 
and multiplies it by the reference spectral value 
and further by a convergence factor. This is then 
subtracted (algebraically) from the original ampli- 
tude, and the new drive signal is sent out to the 
external load at the updated frequency. The servo 
then looks again at the incoming response, corrects 
for the error, and outputs a further cycle. This 
continues until the frequency reaches the high limit 
of the control bandwidth. 

Mathematicallv, the function of the amplitude servo 
in a digitally controlled sine sweep testing system 
can be summarized in the following equation. 


AMPL = AMPL - REF • (CONT - REF) • o 

(69) 

where 

AMPL = amplitude of driving signal 
REF = amplitude of reference level 
CONT = amplitude of control level 
= convergence factor 

A word may be added about the convergence factor. 
The amplitude of the control spectrum is required 
to converge fast towards that of the reference spec- 
trum at all frequencies. However, no overshoot 
(or undershoot if control was previously greater 
than reference) should occur. So the case is here 
similar to a critically damped oscillator. It is 
therefore at once apparent that the convergence 
factor depends upon the whole servo system and 
test article dynamics. In fact, one of the param- 
eters in the digital sine control is the compression 
speed of the servo. This m'ght be low, normal, 
or high, causing some small occasional overshoots 
(or undershoots). The value of these convergence 
factors is not the same (for a certain nominal com- 
pression speed) at all frequencies, i. e. , a is a 
function of frequency. Consequently, look-up tables 
are generally provided within the algorithm to 
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supply those values. It would otherwise take 
too much time to generate them in real-time. The 
dependence of a on frequency will usually range 
anywhere from a linear function to a logarithmic 
function. Figure 41 shows a typical servo speed 
function in dB/sec vs log frequency for the cases 
of a 1 to 10 gain step (marked +) and a 10 to 1 gain 
step (marked -). Two different compression 
speeds are marked High (H) and Low (L). It is 
seen that a faster convergence is achieved in step- 
down rather than stepup jumps, and at higher 
rather than lower compression speed. 


Solving for ^f, this yields 

Af = - I ] (71) 

Now if gj and arc respectively the high and low 
amplitudes of the two levels, then 


g. Control to a step change in reference spec« 
trum . Since the convergence speed of the servo is 
finite, it is clear that programming an infinite 
slope step will never be controlled since this 
requires a very high servo speed, resulting in a 
high overshoot (or undershoot) which will probably 
exceed the abort limits preset in the test setup 
parameters. Steps should therefore have a finite 
frequency interval, so that a frequency profile 
switches from a level Mj at frequency f to a level 
M 2 at a frequency f + Af, where Af is nonzero. 
Figure 42 shows such a change. 


Substituting back in £q. (71) yields 


Af = 



1 


where y is given by 


Y 


20 > N * £n 2 
2.303 - 60 • S 




(73) 


(74) 


The Af is a function of the magnitude of the step, 
the frequency at which this step occurs, the sweep 
rate, and the servo speed. The higher the magni- 
tude, frequency, or sweep rate, the larger Af must 
be to enable the servo to control the jump. Con- 
versely, the higher the servo speed, the lower the 
necessary Af required for satisfactory control. 

To calculate Af, we proceed as follows: 


Let: 

N = sweep rate in octaves /min, 

M = step change in voltage dB, 

S = servo speed in dB/sec, 
f = step change frequency in Hz 

The number of octaves covered in going from 
frequency f to frequency f Af is 



The dependence of S on frequency is typically that 
shown in Fig. 41, So if N, M, and f are known, 

Af can be calculated. This should be used in an 
optimum system to calculate the convergence factor 
that will achieve the fastest realizable convergence 
without overshooting or undershooting. To this end 
the computer can generate a look-up table for y 
once N is known from the input parameters. The 
reference spectrum may then be scanned and all 
the steps defined in terms of magnitude and fre- 
quency. A Af may thus be calculated for each 
step and then used to determine the convergence 
factor by linear interpolation. This in turn deter- 
mines what amplitude the driving signal should 
have. 

3. Transient Control . The term ’’transient con- 
trol test" is defined as a vibration test that induces 
an oscillatory pulse of short duration that satisfies 
physical continuity to the specimen. This is also 
known as "transient waveform control" (TWC). 
Typical durations of these tests are in the order of 
1 sec or less, thus implying open-loop control. 

The classical shock pulses such as sawtooth, half- 
sine, rectangular, and triangular waveform as well 
as more general and arbitrary transient time his- 
tories are synthesized at a shaker control point. 

Since the introduction of digital transient control 
techniques for electrodynamic shakers by Favour, 
LeBrtxn and Young in 1969^^®), numerous labora- 
tories have conducted transient tests using simi- 
le r digital control systems. This opened a truly 
new capability for the digital system where there 
is no close parallel in the analog counterpart. 
Reproduction of transient waveforms on electrody- 
namic or hydraulic shakers is straightforeward 
and accurate using digital computers, provided that 
the test system under control is approximately 
linear. Recent studies of Hunter(^^) and 
Barthmaier( ^3) explore the digital capability in 
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highly nonlinear situations. The control methodoU 
ogy may be divided into the two following groups: 

a. Specific transient time history . Given a spe- 
cific transient time history, a digital system can 
equalize and synthesize the desired waveform, by 
compensating for the test system transfer func- 
tions^). Software requirements are shown in a 
flow diagram in Fig. 43. 

The main theory behind the TWC process is that of 
computing an instantaneous transfer function exper- 
imentally, assuming linearity of the system. It is 
most crucial to obtain a meaningful transfer func- 
tion within the dynamic range of the digital equip- 
ment. In order to guarantee the full frequency 
response of wide bandwidth, a delta function can be 
chosen for the initial calibration pulse fi(t). In 
reality, an ideal delta function is simulated by a 
finite -amplitude square-wave pulse with extremely 
short duration. 

The direct Fourier transform of f^(t) yields F|(w); 


Equation (78) appears to be straightforward algebra, 
but it requires particular care in practice. First, 
one must send out f-Ct) such that the physical system 
will produce a healthy frequency response in the 
bandwidth of the desired transient waveform. This 
requires the understanding of the physical limita- 
tions of the shaker -amplifier system. It is clear 
that one must analyze the desired test pulse wave- 
form x(t) beforehand to obtain its frequency band- 
width, peak acceleration, velocity, and displace- 
ment levels such that all the requirements fit withi»i 
the shaker limitations. Electrodynamic shaker 4 
behave like high-pass filters up to 3000 Hz* Lne 
hydraulic shakers behave like a low-par d filter. 
Second, the complex division of Eq \ 18 ) must be 
done within the dynamic range of the digital 
machinery. 

Once H(u>) is accurately determined, the required 
input waveform y(t) is synthesized in order to 
achieve the desired waveform x(t) at a specified 
location on the test specimen. The Fourier trans- 
form of the required transient wave y(t) is 

’'"■■I ■ Ss{ 


F.(w) 




f.(t) 


-jwt 


dt 


where X(u>) = »^[x(t)]. 



= — [sin wT + j(cos o)T - 1)] (76) 

U) 

and 

|f.(w) I = 4 yz ( 1 - CO 8 IJT) (76) 


The transform pair, functions fj(t) and |F|^(u))|, is 
shown in Fig. 44. If t is selected as one sample 
period of the D/A converter, then tt/2t is the 
total band vidth of F^(w). The choice of a delta 
function as the calibration pulse is not unique, 
although it is found to be most convenient for 

digital controllability. 


The calibration pulse f.(t) is transmitted to the 
shaker through the D/a converter, and the 
test specimen response signal fQ(t) is detected 
by a transducer and stored in the computer 
memory through the A/D converter. A real- 
time Fourier transform is performed on f.(t) 
and fjj(t): ‘ 

F.(u) = .F[f.(t)] 

FqIw) = .#-[£o(t)] 


(77) 


where is the Fourier transform operator. 

Then the transfer function is obtained through 
complex division, 


H(u.) 


^o<“> 

F,(u) 


(78) 


Then 


y(t) = JS”* [Y(«)] 


(80) 


All the frequency functions are complex functions; 
thus all the Fourier transform pairs constitute one- 
to-one invertible mapping. This eliminates the 
arbitrary phase assumptions which are associated 
with test specifications given as frequency spectra. 

Computational roundoff errors, truncation errors, 
and physical nonlinearity errors will occur in 
practice, and it is important to obtain quantitative 
information on the relative magnitude of the error; 
that is, the required function y(t) will not exactly 
yield the ^^.esired response x(t) 

I.et x(t) be the specimen response signal as y(t) is 
transmitted to the shaker. The error signal e(t) 
can be defined at each discrete point i in time as 

e. = X. - X., i = 1,2, • • • . N (81) 

N = digital frame size 

The error can be given as statistical quantities, 
i.e., a mean and a variance. The mean of the speci- 
fied wavefoim is 


N 



i=l 


The mean of the resulting error is 


N 



i=l 


(82) 


(83) 
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The associated variances are 


N 




) (84) 


i=l 


and 


N 

2 1 / .2 
% = ‘ 
i=l 


(85) 


Thus a measure of the error in time domain can be 
defined as 


percent variance error = 100 



( 86 ) 


For typical transient pulses generated on electro- 
dynamic shakers, the percent variance errors are 
in the neighborhood of 10%, and this is considered 
acceptable. TWC allows for improved control 
tolerances, which means that acceptable errors 
and tolerances can be realistically specified for 
this type of shock testing. 

b. Shock spectrum synthesization. This requires 
an additional computation to determine the exact 
transient time history that satisfies the shock 
spectrum. Since a shock spectrum cannot uniquely 
define a corresponding time history, there are 
several proven methods of achieving the goal. One 
method(23) is to use a series of symmetric half- 
sine wavelets to synthesize the given shock spec** 
trum. In due process, it allows the user to specify 
each wavelet's duration and time lag, thus giving 
control over general waveform shape in the time 
domain. This feature is desirable in some 
instances where the test specification requires time 
domain constraints in addition to the shock spec- 
trum. For example, total duration and the envelope 
shape may be specified as well as the ntimber of 
level crossings for fatigue consideration. The 
half-sine wavelet method also meets the shaker 
compatibility conditions such as the Initial and final 
displacements and acceleration to be zero. 


Each of the wavelets is an odd number of half cycles 
which has been v ighted by a half sine wave. Each 
of the wavelets may also have a delay with relation 
to the start of the data buffer. Two individual wave- 
lets are shown in Fig. 45 and the summation of a 
group of these wavelets to form a composite 
waveform is shown in Fig. 46. Thi^ composite 
waveform was synthe.sized to match a specific 
shock spectrum. 


W(t) = 


M 

§ 




f = 0 
m 


0 < t < t 


dm 


= ,in 2n b^(t - 2TtN^b^(t - 


‘dm " ‘ " ‘dm " 


‘ * ‘dm ^ ^m 


(87) 


where the frequency of each wavelet, A f , is 
^m^m* frequency of its half sine envefope is 
bml and is the number of half cycles. Then, 
Tm» the duration of the wavelet, is Tm = l/2b, 
and t^^ is the delay of the wavelet. 


m 


A solution for the amplitudes of the wavelets, Aj^, 
must then be found such that the shock spectrum 
of the composite waveform, W(t), is arbitrarily 
close to the specified shock spectrum. An approxi- 
mation of this wavelet amplitude, A„^, is formed 
by dividing the desired shock spectrum value at the 
wavelet frequency, Nj^b^^, by the number of half 
cycles Njy^. This approximation is used as a 
starting point for an iteration procedure which 
calculates the amplitudes, Arn» to the desired 
accuracy. Accuracies of 3% or less may generally 
be obtained in 3 or 4 iterations. Yang and Saffell 
discuss some of the restriedons necessary in 
defining the delays, t^j^, the number of half 
cycles and the spacing of the wavelet 

frequencies N b 
^ mm 


In one implementation, the delay and the number of 
half cycles are variable. This allows the operator 
a great deal of flexibility in defining the shape of 
the time history used The waveform may be made 
to look very much like earthquake, pyrotechnic, or 
other shock waveforms. 


E. Test Documentauon 

One of the most significant advantages of a digital 
vibration control system over an analog system is 
the ability of the mainframe computer to ’’remem- 
ber" everything that happened during the test by 
storing in its core memory, or on a mass storage 
device, the information gathered during the test 
run. The system is therefo.”e able to display at the 
end of the test a complete documentation of the 
conditions under which the test was actually run. 
This may include important information as to the 
maximum deviations of the control spectrum from 
the reference spectrum in a random or a sine -wave 
test and the frequency at which it happened. It may 
also indicate what the test duration and actual con- 
trol bandwidth was, wliether any of the abort or 
alarm lines have sensed any abnormal conditions, 
and so on. 


Using Yang and Saffell' s notation, the composite 
waveform is described inathematically by the fol- 
lowing equation: 


This so-called post-test documentation is especially 
useful when a test is aborted in an abnormal way, 
such as an operator manual abort, an abort due to 
abnormal conditions sensed on one of the abort 


lines, or a test specification that was exceeded 
during the test run. This allows the test operator 
to pinpoint the cause of the failure of the test and 
correct for it by either modifying the test specifi- 
cation or by checking the device connected to Ihe 
activated abort line, or by taking any other appro- 
priate measures. 

Additional information pertaining to random and 
sine testing might be a mention of which channel 
was in control over what frequency range. This is 
especially important in multichannel sine control. 
Two examples showing post-test documentations of 
normally and abnormally completed tests are shown 
in Tables 5 and 6. Both are for typical sine -sweep 
tests. A record is maintained of the test name 
and/or heading, the completion status (normal or 
aborted and the reason for aborting the test if in 
fact it was aborted), and which line caused the test 
to abort if it was aborted, "the report also shows 
the sweep number and frequency at which the test 
was aborted and the total duration of the test. Other 
documentation ^nown in these tables is the maxi- 
mum and average control error (deviation of control 
spectrum ^rom reference spectrum in dB) and the 
frequency ct wh’vh the error occurred. Also shown 
is the control .listory over the whole test run, that 
is, which input channel was controlling, and when 
it was controlling, its frequency range. 

Table 5. Typical post-test documentation 

for a completed sine test 

(Test duration: 6 min; maximum absolute control 
error; 0.89 dB at 661 Hz; average absolute control 
error: 0.147 dB) 


Control 

Frequency range 

channel 

Hz 


Sweep 1 


1 

20 “ 186. 3 

2 

186.3- 212.1 

1 

212. 1 - 859.1 

2 

859. 1 - 886.5 

1 

886. 5 - 2000 


Sweep Z 


1 

2000 - 

882. 1 

2 

882. 1 - 

855. 6 

1 

855. 6 - 

210. 2 

2 

210.2 - 

184. 8 

1 

184. 8 - 

20 


Sweep 3 


1 

20 - 186.5 

2 

186.5- 212 

1 

212 - 859.9 

2 

859.9 - 886.5 

1 

886. 5 - 2000 


Table 6. Typical post-test documentation 
for an aborted sine test 

(Completion status: aborted during sweep 1 at 
247. 3 Hz, abort line 5; test duration: 6 min, 

12 sec; maximum absolute control error: 0.71 dB 
at 53. 66 Hz; average absolute control error: 
0.067 dB) 



Control 

channel 

Frequency range, 
Hz 

Sweep 1 


1 

20 - 186. 1 


2 

186. 1 - 211. 8 


1 

211, 8 - 247 . 4 


V. CRITERIA FOR SELECTION, INSTALLATION, 
AND MAINTENANCE OF COMPUTER - 
CONTROLLED TEST SYSTEMS*^ 

The advent of digitally controlled systems for 
environmental testing and analysis will substan- 
tially affect the scope, operation, and personnel 
structure of environmental test laboratories 
employing such devices. The degree of impact will 
depend upon the types of applications to be sup- 
ported. Although only a few of these systems exist 
at present, dozens will be purchased in the future 
for test operations and control, analyses, and other 
related applications. The advantages of the digital 
test system over the older analog systems have 
been described in the literature(^) and are listed in 
a prior section of this paper and will not be repeated 
here except to say that the digital systems will 
replace the analog systems because the digital sys- 
tem functions are computer -controlled and are 
adaptable to many applications. The analytical 
capability of these systems can often justify their 
procurement almost independent of the control 
function. 

The purpose of this section of the paper is to 
suggest applications, list hardware, software, and 
personnel requirements, and discuss pertinent 
questions that must be answered and facts that 
should be known prior to designing or specifying and 
procuring computer -controlled environmental test 
systems. 

A. System Characteristics and Applications 

g 

1. Induced Random Vibration Test Systems The 
first computer-controlled test system was designed 
for random environment simulation!^*’ ^ Computer- 
controlled sine -wave test systems, a later develop- 
ment, are discussed in the following section. 
Random-type computer test systems are apparently 


^This section was presented originally in Ref. 24. 

g 

The terms ’’induced" and "induced environment" as used in this section refer to a mechanically produced 
or simulated environmental effect, such as a mechanical vibration or shock. There is no implication 
that the cause of an induced environment necessarily originates in a mechanical process 
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creating the most interest and greatest demand at 
present. Unfortunately, this category of computer 
control has been the hardest to implement. The 
general concepts for random test systems are 
described in this paper, and the basic algorithms 
to implement these concepts were originally pre- 
sented in Ref. 19. 

All of the original random -type digital control sys- 
tems were designed around existing main subsys- 
tems that were available during the construction of 
these test systems. The basic control system 
hardware can be divided into three main subsys- 
tems: a Fourier processor, a computer, and an 
interface controller that may be part of or separate 
from the first two subsystems. 

An interesting and probably worthwhile control 
concept modification was suggested by Sloane^^^l 
In this method of random vibration control, cross - 
spectral density techniques are used to estimate 
the transfer function of the vibration system 
including the test item. Spectrum control is 
achieved at the points s) of observation so that the 
output spectrum consists of two parts: a control 
spectrum representing the environment to which 
the structure is tested and a second spectrum 
representing uncorrelated self-noise generated 
from test specimen nonlinearities. Only a so^.ware 
change is required to implement this concept. 

Computer -controlled random environment test sys- 
tems should include the following multi-' ransducer 
channel functions: 

(1) Spectrum or spatial averaging control!^^* 

(2) Alert or shutdown for signal loss. 

(3) Alert or shutdown for spectrum out of 
specification. 

2. Sine -Wave Test Systems . Computer -controlled 
sine test systems are less expensive than compute r- 
controlled random test systems because the Fourier 
processor is not required for the sine system. 
However, another subsystem is required: the 
digitally controlled oscillator(®» H) or some 
other method of producing and controlling (in 
frequency and amplitude) a sine wave under com- 
puter control. 

The role of the computer for sine -wave testing is 
nebulous. One manufacturer suggests letting the 
computer provide a supervisory role only(^) and 
has provided special-purpose hardware to generate 
and control the sine-wave forcing function. The 
author believes that the system computer should be 
made to provide and control as many functions as 
possible for whatever application is needed, and 
that this capability should be achieved through the 
use of software rather than by attaching peripheral 
hardware, which decreases the system reliability. 
Such a solution holds down hardware costs but 
increases software costa. Perhaps an optimum 
compromise to this dilemma is to let the computer 
determine and control the frequency and amplitude 
of a sine -wave function generator as implemented 
by Norin^^®). In any case, a computer-controlled 


sine -wave test system should be able to provide the 
following multi-accelerometer channel functions: 

(1) Peak control. 

(2) RMS (or other type of averaging) control. 

(3) Displacement control. 

(4) Velocity control. 

(5) Acceleration control. 

(6) Multivibration-level programming. 

(7) Signal loss alert or shutdown. 

The system should also be capable of acting as a 
peak limiter for each channel of control, independ- 
ent of the functions above. If any one, some, or all 
of the control accelerometers exceed a predeter- 
mined level, the system should shut itself down in 
a controlled manner and advise the operator which 
accelerometer(s) exceeded the predetermined level. 

Ai with the random test system, the test specifica- 
tion for conducting a sine -wave test is entered by 
means of a pretest conversational language 
program. 


At present, the additional cost factor for a sine 
system over an n-channel (n > 1) random system 
is 1.1. That is, if the random system costs 1 unit, 
1.1 units will buy the sine system. 


Transient waveform testing capability may add to 
the cost, in the form of additional software; tran- 
sient waveform testing also requires a Fourier 
processor. If the required transient waveform is 
produced by the test system computer, however, no 
new hardware will be required If the vendor pro- 
vides the transient waveform software "free of 
charge, " then the cost factor of 1.1 applies for an 
n-channel random system with sine, shock, and 
transient waveform capability. 


3. Large-Scale Data Analysis and Presentation. 
Large-scale time series analysis and presentation 
is an extremely worthwhile by-product of the 
induced environmental test system This type of 
application is made possible by the Fourier proces- 
sor associated with the random-type control system 
(for a detailed description of this type of analysis 
and of the basic frequency function algorithms, see 
Ref. 29). The following time and frequency domain 
analyses are easy and inexpensive to implement: 

(1) Time domain data analysis and presenta- 
tion: auto-correlation, cross -correlation, 
and probability density and distribution 

(2) Frequency domain data analysis and presen- 
tation: auto spectrum (power spectral 
density), cross spectrum, transfer function, 
coherence function, and shock spectrum 


The computer-controlled test systems can support 
the large-scale analysis application merely by 
additional software and a high-speed output device 
that has hard copy capability. 


29 
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4. Natural Environment Testing . The compute r- 
comrolled test system has a dennite two-part appli- 
cation in natural environment testing: for control 
and for data presentation. No new technology is 
required. The Fourier processor is not required 
for this function, but an interface controller is 
required, along with pressure, temperature, and 
humidity controllers. These controllers can con- 
trol chamber environments and provide analog (or 
digital) information about the chamber environ- 
ments. This information can be displayed alpha- 
numerically on a cathode ray tube (CRT) and 
updated periodically. A permanent copy of the 
chamber environment parameters should be made 
periodically as well. Interrogation of the environ- 
mental status of any chamber from the test system 
teletype or CRT terminal would be a functional 
requirement. 

Two principal facto.-s govern the cost cf a natural 
environment computer control system: the system 
as an add-on to a raniom and/or sine system, and 
the system as specific \lly designed for natural 
environment control ar.d presentation. (The latter 
type of system, which provided by some manu- 
facturers, will not be discussed here.) The natural 
environment test system could be an add-on to a 
random system, a sine system, or a random-sine 
system. 

The natural environment test system should be able 
to control and specify the following parameters: 

(1) Chamber pressure control. 

(2) Chamber, fixture (heat exchanger), and 
test item temperature control. 

(3) Chamber humidity control. 

(4) Specification of alarm limits for pressure, 
temperature, humidity, chamber line 
voltage and power interruption, test item 
power supply voltage, ambient room 
temperature, and water, air, and nitrogen 
pressures. 

The following parameters should be available for 
presentation: 

(1) Chamber pressure. 

(2) Foreline pressure. 

(3) Backing pressure. 

(4) Chamber, fixture, test item, and ambient 
temperatures. 

(5) Chamber humidity 

(6) Pump status: on/off and valve positions. 

(7) Firm power and test item power status. 

5. General-Engineering Use . The basic computer- 
controlled test system makes an excellent general- 
engineering tool for solving those types of 
engineering problems associated with the 


environmental sciences, including mechanical and 
electrical engineering analysis and synthesis. This 
capability results from th«i general-purpose com- 
puter associated with these systems. All of the 
computers used thus far will support the BASIC 
user -oriented language, and some computers used 
can also support FORTRAN IV or a form of 
FORTRAN, although not conveniently without some 
sort of mass memory storage. Therefore, it is 
recommended that mass memory be considered for 
this application and the application to be discus -ed 
below. 

An alternative to the mass storage hardware is 
interfacing the computer to a larger computer sys- 
tem. This is called multiprocessing and makes 
the control system computer an extremely powerful 
tool at very little additional cost. Computer-to- 
computer (also called processor-to-processor) 
data exchange provides two computers, each with 
local file space, to communicate in such a way that 
each computer may access data stored in the other. 
This allows the control system to perform such 
tasks as interactive fixture design and to perform 
system analysis. However, 12,000 words of 
memory will allow engineers to do most of the run- 
of-the-mill types of mechanical and electronic 
analysis and synthesis problems required in test 
and control. 

6. Laboratory Management Use . The basic 
control system used as a laboratory management 
tool requires mass storage capability or computer- 
to-computer interfacing. The term '’laboratory 
management tool" is used to denote a method of 
implementing the mathematics of system analysis, 
which includes (but is not limited to) the following 
management sciences: 

(1) Linear, nonlinear, and dynamic 
programming 

(2) Reliability studies. 

(3) Statistical analysis and probability theory. 

(4) Decision and game theory. 

(5) Queue theory, including test and design 
scheduling. 

(6) Information theory. 

(7) Forecasting. 

(8) Program evaluation and review technique 
(PERT) (critical path method) and networks. 

(9) Human engineering 

(10) Financial analysis 

References 30 - 32 discuss most of these topics in 
a general way. Computer programs have been 
written to implement the algorithms associated 
with these management tools 

The control system computer and ^nput /output (I/O) 
devices can be used for general administrative 
tasks such as accounting and inventory and 
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The terms "natural" and "natural environment" as used in this paper refer to environmental effects that 
occur in nature, including temperature, vacuum, humidity, sunlight, sand and dust, fungus, salt spray, 
or any combinations of these. 
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equipment control. The system can also generate 
test procedures and modify or update these proce- 
dures, as required, since the computer manufac- 
turer includes in his software a means of updating 
and modifying source programs on a line-by-line 
basis. This means that long lists and text can be 
modified more easily and quickly by the computer 
than by a typist. 

7. Other Applications . Only the creativity and 
imagination of the vendor and the potential user 
limit the applications of computer-controlled test 
systems. Besides the above -listed applications, 
these systems are being used as transducer cali- 
brators, earthquake simulators, road sim\ilators, 
mechanical impedance analyzers, and nuclear 
reactor monitors. 

The number and type of applications will determine 
the size, complexity, and cost of the computer - 
controlled operating system. (An operating system 
is defined here to mean a test system under control 
of an executive program such that the operator 
merely specifies the test function(s) desired. The 
system does the rest, obtaining information from 
the operator as required.) 

Some of the applications above can be (and may be 
required to be) supported simultaneously. For 
example, the natural environment test control 
functions, as described above, can be carried on 
in parallel with an induced environment test. How- 
ever, if several parallel operations are necessary, 
two (or more) computers may be required. 

8. Criteria for Determining Applications . How 
does one recognize the requirements for computer- 
controlled systems? The answer, of course, is in 
knowing one's own testing requirements and knowing 
the advantages and capabilities of computer - 
controlled test systems, as well as carefully 
considering the following factors: 

(1) The test function performed by the environ- 
mental test laboratory. If the function is to 
provide induced random and/or sine -wave 
simulation and/or shock pulse or shock 
spectrum environments to complex test 
items requiring complex fixturing, then the 
laboratory may be a likely candidate for a 
computer- controlled test system. If the 
test program requires testing m items 
using n test specifications, where m is a 
relatively large number, and there is truly 
a need to ensure test repeatability by 
implementing the n specifications several 
times during the test program, there is an 
additional argument for the computer - 
controlled test system. If, in addition, any 
appreciable amount of data reduction, 
analysis, and presentation is required for 
test documentation or test specimen design, 
development, and reliability or structural 
behavior studies, there is all the more 
reason for considering computer-controlled 
test systems. If the laboratory is also 
performing natural environment testing on 
the same m items, there may be a need for 
an operating test system to provide all the 
test functions. 


(2) The requirements that may make the sys- 
tem feasible. Two conditions under which 
a computer -controlled test system may be 
justified are: when the procurement or 
construction of an automated test system 
pro ’ides the laboratory with greater teat 
control features and when operating costs 
may be reduced, thereby achieving a less 
expensive program. The following are 
additional considerations in making such a 
choice: 

(a) Increase in test volume (that is, a 
greater number of tests in a given 
period of time) for the same test 
facilities and space. 

(b) Safer and more accurate testing with 
reduced setup time and fewer human 
errors. 

(c) Repeatable testing and testing that can be 
exactly duplicated at other laboratories. 

(d) Data analysis capability without 
depending on external (department, 
plant, outside agency) facilities. 

(e) Automatic post-test performance 
documentation. 

(f) Flexible capability of implementing 
new or different test philosophies and 
expanding the test capability merely by 
changing or modifying software. 

(3) The existing test system control capability 
and the system's limitations and problems. 

If the existing system is adequate for the 
present and near-future requirements, it 
may not be possible to justify a computer- 
controlled test system at present. On the 
other hand, arguments for justifying the 
computer may lie in situations where, (a), it 
is necessary to subcontract such functions 
as fixture design and data reduction, anal- 
ysis, and presentation; (b), the laboratory 
is una^ le to do all the testing that is 
required; (c), the tests are based on what 
can conveniently be done, rather than on 
what should be done; or (d), the addition of 
a new system would increase the test con- 
trol cap.»bility. 

(4) The possibility that a new control system 
must interface with any part of the existing 
control system(s). A computer may create 
new control problems if it must be inter- 
faced to existing hardware. On the other 
hand, the old system can possibly be used 
in its existing form as a parallel redundant 
(but limited) test control system, completely 
independent from the new system. 

B. System Hardware and Software 

1. Add-on Hardware and Software . Table 7 lists 
the hierarchy of hardware and software for the 
types of applications that can be supported by 
computer- controlled environmental test systems on 
an add-on basis. 
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Table 7. Computer teat systems hierarchy of hardware and software 


Application 


I. Basic system, 
n channels of test 
control 


II. Random system, 
n channels of test 
control (for 
induced vibration 
and acoustic 
testing) 


III Large-scale data 
analysis and 
presentation 


Requirements 


A. Computer system 

1. General-purpose 
computer 

2. High-speed paper 
tape reader 

3. High-speed paper 
tape punch 

4. Analog multiplexer 
and A/D converter 

5. Priority interrupt 
lines 

6. Direct multiplex 
control/direct 
memory access 
channels 

7. Group I software 
package (equipment 
manufacturer’s 
software library) 

8. 8000 core memory 

B. Conditioning and 
conversion equipment 

1. Input /output 
amplifiers / 
attenuators 

2. D/A converter 

3. Input anti-aliasing 
filters 

4. Output low-pass 
filter 

C. Input /output equipment 

1. CRT/keyboard 
terminal with hard 
copy capability 

2. Group I cable and 
connector package 

A. Group I hardware and 
software 

B. Fourier processor 
(dual channel) 

C. Interface controller, 
version 1 

D. Additional 4000 
memory (12,000 total) 

E. Group II software 
package 

F Group II cable and 
connector package 

A. Group II hardware 
and software 

B. Group III software 
package 


Application 

IV. Sine (n channels) 
shock and tran- 
sient waveform 
analysis and 
control 


V. Natural environ- 
ment test system, 
n channels (less 
chamber 
controllers) 


VI. General 

engineering tool 


VII. Laboratory 

management tool 


Requirements 


A. Group I hardware and 
software 

B. Additional conditioning 
and conversion 
equipment 

1. Computer - 
controlled function 
generator 

2. Peak control 
hardware 

3. Average control 
hardware 

C. Interface controller, 
version 2 

D. Group IV software 
package 

E. Group IV cable 
and connector 
package 

A. Basic system 

B. Multiplexer expansion 
for computer output 
control 

C. Real-time clock added 
to computer 

D. Group V software 
package 

E Interface controller, 
version 3 

F. Group V cable and 
connector 
package 

A. Basic system 

B Additional 
4000 memory 
( 16,000 total) 

C. User-generated 
software 

A. Group VI hardware 

B Mass storage capabil- 
ity or link-up capabil- 
ity to large-scale 
computer 

C, Additional equipment 
manufacturer’ 8 
software 

D. User-generated 
software 

E Optional high-speed 
input device 


A basic system is defined in terms of three 
elements: the computer, conditioning and conver- 

sion equipment, and input/output devices. These 
elements, consisting of both hardware and soft- 
ware, are required independent of the application 
and, whatever applications are required, additional 
hardw'are and software must be added to the basic 
system. 

The next system in Table 7 s defined for induced 
random-excitation environmental testing. This 


application was chosen second to the basic system, 
rather than a sine-wave test system, because of 
present interest and demand. The la^ge-scale data 
analysis test system is essentially a by-product of 
the random test system and requires only some 
addi' ■ inal core memory and software. The sine- 
wave test system also requires additional hardware 
and different software but does not require the 
Fourier processor used in the random test system 
However, it may be possible to utilize the Fourier 
processor for sine-wave testing at some future 
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time since the hardwired Fourier proceaaora con* 
tain sine (and cosine) tables within ''read only 
memory" in order to implement the fast Fourier 
transform algorithm. A different version of the 
interface controller is required for thr applications 
of random and sine test systems as well as a natu- 
ral environment test system. These different 
versions are really insignificant in terms of tech- 
nology and complexity* and could easily be incor- 
porated into a single controller at very little if any 
additional cost. 

The main additional hardware for the natural envi- 
ronment test system is an expansion of the multi- 
plexer capability in order to control n number of 
chambers and m fvmctions for each chamber. In 
addition* a real-time clock is necessary in this 
application in order to conveniently present the 
chamber parameters (environments) at regular 
intervals of time. 

The test systems used as general engineering and 
management tools require only additional computer 
memory storage capability and specialized soft- 
ware, mostly user-generated. 

2. Control System Software. *A key question after 
the applications are chosen concerns the status of the 
required software; that is* are the software logic 
programs available to satisfy the proposed appli- 
cations? The specific applications? Or will the 
vendor develop this software while he is putting the 
system together? Or should the user's program- 
mer develop it? 

The software logic programs have been developed 
for "standard" random sine -wave control* and 
transient waveform control and analysis. Shock 
spectrum testing philosophy is as confused in the 
digital world as it is in the analog world. 

To the oest of the authors' knowledge* there has 
been no software development for natural environ- 
ment testing by the vendors in the business of 
supplying these new computer -controlled environ- 
mental teat systems* but some software develop- 
ment has been done by vendors supplying computer 
control systems specifically for chamber control. 

As special applications appear* special programs 
will be developed. If every vendor used the same 
equipment* a program bank could be established* 
but this is not likely to happen. The closest 
approach to such cooperation thus far on these 
systems is the use oi the Digital Equipment Corp. 
PDP-11 computer* although this is not universal. 

One factor is common to all the present control 
systems: the logic programs are generally written 
ill assembly language and, when these systems are 
procured, the buyer receives system tapes to 
implement the control functions. These system 
tapes were generated by the vendor using the pro- 
cedure in Fig. 47. If the user's organization struc- 
ture includes a programmer, as the author recom- 
mends* then the control system specification 
must include the source tapes used in generating 
the system tapes* whether the tapes are proprietary 
or not. 


All of the post-procurement programs should be 
written in the assembly language supported by the 
system assembler. (Higher-level computer lan- 
guages such as FORTRAN are probably too ineffi- 
cient for real-time control.) The vendor should 
supply updated programs and listings free of 
charge. The computer manufacturer will send the 
user updated object tapes and source listing for 
the computer. There may be a charge for these 
programs. 

C. Constraints and Tradeoffs 

1. General Specifications . The system is speci- 
fied when the applications are defined. However* 
there is more than one vendor available to build an 
operating system, and each vendor has a seemingly 
endless list of options. This section discusses 
some important parameters of the system and 
some of the input /output and memory options* along 
with other constraints. 

The five most important requirements of the total 
system are* in order of importance: 

(1) The system mu,it be flexible enough to meet 
current and anticipated needs. 

(2) The system must be safe to operate and 
work around. 

(3) The system must be human-engineered for 
ease of operation* maintenance* and 
personnel training. 

(4) The system must be adaptive and easily and 
economically expandable, both from the 
standpoint of hardware and of software. 

(5) The system must have a practical capital 
cost and practical cost of operation. 

2. Input -Output Options The I/O options probably 
deserve the most concern and contribute signifi- 
cantly to the overall cost. Table 8 lists some of 
the more pop\ilar I/O devices. A choice may have 
to be made between a teleprinter such as a Teletype 
Corp. ASR 33 or 35 and a CRT terminal such as a 
Tektronix 4010-1 with the 4610 hard copy unit (The 
designation "ASR" means "automatic send- receive.") 
The 33 machine, which comes "standard" with most 
of the computer-controlled teat systems, is made 
mostly of plastic and is inexpensive. It is a good 
machine for periodic use. but it is not a continuous- 
duty machine. A user who needs a continuous « duty 
machine should specify the ASR 35, which is made 
of metal and is more than twice as expensive as the 
ASR 33. On the other hand, the 3 3 is more conven- 
ient to use than the 35 for correcting and duplicating 
short programs. The author recommends an 

ASR 33 if the user is specifying a high-speed reader 
and punch and if he plans to obtain some other high- 
speed output device. A user who is teletype -limited 
should specify the 35 over the 33 

These teletypes (TTYs) will input data to the sys- 
tem from either the keyboard or prepunched paper 
tape (they have slow-speed feeler wire-type 
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Table 8. Input /output devices suitable for computer -controlled test systems 


Device 

Cost 

Functions 

Significant limitations 

factor 

Input 

Output 

Input 

Output 

TTY M33 

1.0 

Keyboard, paper 
tape reader 

Serial printer paper 
tape punch 

Slow 

Not continuous duty 

TTY M35 

2.4 

Keyboard, paper 
tape reader 

Serial printer paper 
tape punch 

Slow 


CRT termi- 
nal with 
hard cop) 
capabilit'^ 

4.0 and 
up 

Keyboard 

Alphanumeric CRT 
display with hard 
copy capability 

No tape reader 

No tape punch 

K«gh-spoed 
paper tape 
reader 

2.0 

Paper tape 
reader 




High-speed 
paper tape 
punch 

1.8-2. 5 


Paper tape punch 


Cost function 
is propor- 
tional to 
speed 

Combina- 
tion high- 
speed 
reader and 
punch 

2.0 

Paper tape 
reader 

Paper tape punch 

200 characters /s 

50 

characters/s 
(not really 
high speed) 

Plotter 

0. 7-5.0 


Alphanumeric graphic 
presentation 


Interfacing 
and D/A 
conversion if 
plotter is 
analog type 

High-speed 

printer 

3.0-10.0 


Serial or parallel 
(line) printer 


Requires 
cache storage 
for full 
speed 
advantage 


readers) and will print out data from the system 
onto paper tape and punch the data onto tape. It is 
also possible to punch a tape from the keyboard 
off-line: that is, when the TTY is not electrically 
connected to the system. 

The TTYs are slow (about 10 characters/sec) and, 
except for very short source programs, are 
virtually useless for a paper tape input device. 

A low-cost CRT terminal costs slightly more than 
four times as much as an ASR 33 TTY. This 
graphic terminal with the hard copy option will do 
everything the TTYs do and more, except that it 
will not punch or read tape. The hard copy option 
will allow permanent copies of the CRT display. 

The first copy takes about 18 sec. Test specihca- 
tions can be placed into the operating system via 
the CRT terminal keyboard. An optional TTY port 
interface is required with the Tektronix unit men- 
tioned above. The hard copy option cannot be added 
to the standard terminal after the fact. If a user 
will need the hard copy capability at a later date, 
he must specify that option when specifying the 
initial CRT terminal. 

The graphic display capability of the CRT terminal 
is the plus factor over the TTY. Software logic 
programs supported by the control system 


assembler are available from the CRT terminal 
manufacturer for a small fee. These programs 
conveniently allow the systems operator to utilize 
the full capability of the graphic display. 

If the user's organization will include a program- 
mer, as recommended, then a high-speed paper 
tape ptinch will be required for tape program gener- 
ation. (A high-speed paper tape reader comes with 
the control system.) High-speed punching is from 
SO to 110 characters/sec. (High-speed paper tape 
reading is from ISO to 300 characters /sec.) 

It is very desirable to have some type of high-speed 
printer for obtaining source program listings when 
developing new software logic programs. However, 
the CRT terminal could be used for the source 
listings, producing one page every 18 sec, 

72 characters /line, and 3S lines /page. A line 
printer operating at 300 lines /min, 120 characters/ 
line, and SO lines/page is faster. The decision 
must be based on the amount of programming 
required and on the applications. Another consid- 
eration, however, is the generation and cost of 
multiple copies, in contrast to single copies. 

Impact types of printers can make as many as five 
copies simultaneously, and TTYs can make two 
copies simultaneously. The CRT devices can make 
only one copy of the display at a time, but, of 


course, it is possible to copy the display as many 
times as necessary. 

Plotters can also be considered as a means of 
graphically displaying the results of a test or the 
output from the system. There are many devices, 
both analog and digital, that could be used. The 
plotter must be controlled by the computer. It 
should also be able to plot and label the output from 
an analysis program within, say, ?0 sec or the 
system will usually be waiting for the plotter. As 
it is, the I/O devices are the slowest links in the 
control system. 

Input amplifiers and anti-aliasing filters, although 
not truly considered I/O devices in the computer 
field sense, are important analog input devices, as 
discussed previously. Transducer conditioning 
amplifiers should be of the differential-input type 
in order to cancel out common-mode noise riding 
on the shielded cables from the transducers. 

There should be at least 90 dB of common-mode 
rejection at the nominal gain settings of these 
amplifiers. They must not be slew-rate limited 
for shock testing. 

The input anti-aliasing filters set the control band- 
width commensurate with the samplin> theorem. 

The low pass filters should be capable of at least 
36 dB/octave rolloff from the cutoff frequency 
(6 -pole filters). Their characteristics should be of 
the Butterworth type (maximally flat response) for 
vibration control, but they must have linear phase 
characteristics for shock test work: that is, both 
shock spectrum and transient waveform types of 
testing. 

3. Control Loop Cycle Speed Constraints . How 
fast^ust the internal processes of the control sys- 
tem be performed in order to implement the appli- 
cation concept? This is difficult to Answer, but 
the question can be restate'*: "Are ♦^he internal 
processes fast enough to mejt the application?" 

The speed requirement is oi interest only when one 
is using the system for real-time control or real- 
time analysis and, even then, is critical only for 
random-type testing Some of the latest special- 
purpose Fourier processors are so fast that the 
memory cycle time of the computer is becoming 
the limiting speed factor. Speed is important 
because a minimum number of algorithms must be 
executed by the system in order to implement any 
control application. It is desirable to have some 
time left over so that additional algorithms o: 
additional coding within an algorithm can oe added 
to the program 

a. Sine-wave testing . If special hardware is used 
to generate a swept sine wave, and special hard- 
ware is used to implement the concepts of peak or 
average control of several control channels, then 
no particular speed concern is required. However, 
if some of these functions are performed by the 
computer, then it is the vendor's responsibility to 
assure the customer that the processes can in fact 
be carried out fast enough so that his system will 
meet the functional requirements and specifications. 
(Another constraint on the vendor is that his 
algorithms are valid processes, mathematically 


speaking, and generate physically valid results 
\inder the concepts of environmental testing.) 

b. Random testing . The present-day computer- 
controlled random test systems are marginally 
fast enough to implement the real-time control 
function. In fact, it has been necessary in some 
systems to repeat, sequentially, several output 
data frames (time domain output forcing functions) 
in order to provide enough time to perform all the 
algorithms necessary for control. This repeti- 
tion, if done properly, is legitimate from the 
mathematical standpoint and probably does not 
degrade the overall induced random test philosophy, 
although it results in a significantly different audible 
noise from the shaker. 

The random test loop cycle time is affected by 
analysis parameters described in a previous sec- 
tion of this paper. The random system used as an 
analyzer gives the customer a large combination of 
effective filter bandwidths and analysis bandwidths, 
but, used as a controller, the system allows only 
a limited number of choices because of the time 
required to execute the logic programs during the 
control loop cycle. 

Th- number of mtiltiplex control channels is limited 
by the computer and multiplexer and by the maxi- 
mum A/D converter sample rate. Eight-channel 
control with analysis bandwidths up to 2500 Hz is 
available. The number of control channels can 
always be increased by spending more money, and 
30- or 40-channel multiplexing is feasible but 
expensive. Another problem, however, has only a 
philosophical answer. The problem is defined as 
follows. Suppose 30-channel control is required 
for an induced environmental test. Suppose also 
that we sample every 200 psec and generate 
100 spectral lines (200 frame size) for an analysis 
bandwidth of 2500 Hz and an effective filter band- 
width of 25 Hz. How long do we sample each 
response transducer before moving to the next? If 
our control loop cycle time is 80 msec and we 
sample each transducer for 1 input frame (40 msec), 
it wotUd take 2.4 sec to cycle through all the trans- 
ducers. Or should we sample the first transducer 
the first 200 psec, the second transducer the 
second 200 jisec, etc.? In short, what kind of 
spectral average control is proper for the test? 

Me cc itrol loop cycle time is available for 
reverb ration acoust testing because of the 
response time of the reverberation chamber, ^ 
particularly chambers larger than 140 m^(-5000ft\ 
Repeating output frames is desirable in this case 
because the time constants of large chambers are 
in the order of seconds at frequencies below 500 Hz. 
E.iternal generation of energy above 1000 Hz is 
usually not required. Convergence to the desired 
constant bandwidth test specifica,.ion must be very 
slow compared to random vibration tests, or 
instability results. 

c. Natural environment testing . Operating systems 
as described in this^apel-^an easily control and 
present environmental data from at least 24 cham- 
bers. This is an easy task for such a system The 
computer would be busy for only a few milliseconds 
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each minute or eo, unlee a it wae interrogated by 
the eyeteme operator or unleee a TTY rather than 
a CRT terminal wae ueed ior data presentation. 

4. Storage Options . The amount of core memory 
spec ified will depend on the applications for the 
control system. In any case, 16,000-word storage 
should be the maximum upper limit. If moie 
memory is required, slower- speed mass storage 
should be considered (mass memory storage access 
time is slower than core memory accesi 4me). 
Core storage is required for real-time control of 
sine- shock, random, and natural environmental 
testing. Mass storage can be a convenient or 
necessary storage supplement for the general 
engineering and management tool applications. 

If applications justify a neeu for mass storage 
capability, the user shoula specify at least 20 times 
more mass storage than would be specified for core 
storage. Mass storage is available in the form of 
drums, discs, or tape. Table 9 lists and defines 
several types of mass storage devices suitable for 
the applications defined in this paper. The tape 
cassette and cartridge mass storage devices are 
becoming very popular for minicomputer systems 
UiC^«'se they are less expensive than any of the 
other devices. For the same mass storage capa- 
bility, a tape cartridge vinit will cost only about one 
quarter as much as a fixed-head disc, based upon 
present prices. The trend in the minicomputer 
business seems to point toward tape cassette and 
cartridge mass storage. However, a r'^view of 
Table 9 reveals the tradeoffs required for the 
lower-cost tape tinits. 

D, Reliability and Maintenance 

1. Control System Reliability . The general reli- 
ability ^unction^(t) of the control system is given 
by 

B(t) = exp [ -yh(t) dtj (88) 

where h(t) dt Is defined as the conditional proba- 
bility of a failure in the interval t, t + dt, given 
failure -free performance up to time t. The func- 
tion h(t) is called the integrant hazard function, or 
instantaneous failure rate, or, sometimes, the 


force of mortality. Unfortunately, h(t) is unknown 
at the present time. If it is invar ia «t with time, 
then h(t) s \ (a constant) and h(t) is called simply 
the failure rate, in which cs' % 


R(t) = exp(-Xt) (89) 


It should be understood that there is no indication 
that the control system reliability model has an 
exponential form. Two years of experience with a 
computer -controlled random system suggests, in 
fact, that the failure rate is variable. But in any 
case, it is impossible to give mean time between 
failure (MTBF) or mean time to repair (MTTR) 
figures for these b- stems. It is possible to obtain 
reliability figures irom the computer manufacturer 
and perhaps the Fourier processor manufacturer, 
in the case of a random-type control system. A 
simplified series bimodal (failure by open or short) 
reliability model could then be constructed. What 
really is needed is the failure density function, 
which is defined as the probability density of the 
occurrence of failure as a function of time. 

2. Maintenance of the System . The purchasers of 
the c ompute r - c orit r ^le d sy ste m s are obligated to 
maintain and repair them. Some or ah the respon- 
sibility can be transferred to the vendor in the form 
of a service contract, or the responsibility for dif- 
ferent parts of the system can be transferred to 
outside service agencies. The authors recommend 
in-house support capability with a maintenance 
contract on the TTY, and vendor and manufacturer 
support, when required, with no service contract. 
The rationale for these recommendations is 
developed below. 

Thus far, the available control systems as 
described In this paper are available from vendors 
who manufacture only part of their control sys- 
tems, although this will not always be the case. At 
present, if a vendor buys a computer as an original 
equipment manufacturer (OEM!, will the vendor be 
able to maintain and repair it as well as the com- 
puter manufacturer? Hardly. It would probably be 
better (less expensive) to call for the computer 
manufacturer's services when required. It is 
doubtful that the vendor, as an OEM, will ever be 
able to really compete with the computer manu- 
facturer in servicing the system computer. At 


Table 9. Memory storage tradeoff parameters 


Storage device 

Cost 

factor 

Words storage 
per block 

Average 
accc«>s time 

Transfer rate, 
time/word 

Core 

l.O 

4, 000 words ^ 

16, 000 words ma:c 

1. 2 psec 

1. 2 psec 

Fixed head disc 

4, 0 

64, 000 words ^ 

256, 000 words max* 

17 msec 

16 psec 

Tape cartridge 

0. 8 

64, 000 words/track* 
4 tracks max 
90 m (-300 ft) max 

60 sec 

23 -m (-75 -ft) tape at 
19 cm/sec (7 5 in/scc) 

3. 55 msec 

^Maximum as defined in this table means a practical upper 
the paper. 

limit for the applications described in 
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bett» a service agreement between the vendor and 
the manufacturer of the computer could be 
arranged. The same reasoning applies to the 
Fourier processor, if one is required. 

The vendor can service and maintain the smaller, 
less expensive peripheral components, even though 
he may not manufacture them. One or two excep- 
tions might be such items as uisc or drum memo- 
ries, line printers, and perhaps CRT terminals. 


Table 10 lists some of thr preventive maintenance 
functions required for an all-round preventive 
maintenance program for random-type control 
systems. 

Table 10. Some preventive maintenance 
functions for corr r uter-controlled 
random test systems 


Period Maintenance function 


Another problem may be the location of the vendor's 
service facility. How far is the facility from the 
user's location? Will someone from that facility be 
able to service the system or will the factory be 
required to send someone? 

The minicomputers used for test control can be 
serviced in much the same way as the old black- 
and-white TV sets. First, one reads in the diag- 
nostic programs to determine what part of the 
computer is malfunctioning. (In the case of the 
TV set, the symptoms could be seen on the TV 
screen.) Then one can start plugging in spare logic 
cards, analogous to the spare TV tubes. With luck, 
the problem can be found. Reliability studies show 
that the integrated circuits are the least likely 
devices to fail. But the point is, the system ana- 
lyst (if there is one) has a good chance of finding 
and correcting a system failure, depending upon 
his experience and training. Several minicomputer 
manufacturers offer maintenance classes either 
free (for their customers) or at a nominal charge. 


Daily 


3 months 


6 months 


1. Check cooling fans 

Z, R\in control system integrity 
diagnostics 

1. Check and adjust master clock 
frequency 

2. Check all power supply voltages 

3. Inspect, clean and oil TTY 

4. Clean all air filters 

5. Clean and oil high-speed punch 

6. Align optical paper tape 
reader 

7. Clean and adjust line printer 

1. Calibrate all rms and peak 
reading meters 

2. Calibrate oscilloscope and CRT 
displays 

3. Calibrate charge amplifiers 

4. Calibrate I/O amplifiers and 
low-pass fi’ters 

5. Calibrate A/D converter(s) 

6. Calibrate D/A converter(s) 


The computer is not the only subsystem that may 
cause problems. The computer, the Fourier 
processor, and perhaps the interface controller 
are the most complex subsystems, however, in the 
control system. 


If the user plans to do most of his own maintenance 
(as recommended), it is very important to consider 
spare parts and diagnostic programs. 


One precaution is important: if any type of service 
is to be attempted by anyone, manuals, schematics, 
procedures, and diagnostic programs will be 
required. The user should specify them and see to 
it that they are kept up to date (as the system ana- 
lyst's responsibility)^^. The vendor has a funda- 
mental obligation to make this documentation 
available, to make sure it is accurate, and to keep 
it accurate as field modifications to his system arc 
made. This does not mean, however, that he 
should supply this documentation free of charg'^. 

Preventive maintenance which involves more than 
simply repairing hardware should be an in-house 
support capability, except for the TTYs. Teletype 
n‘)'»chincs are mechanical monsters and require 
careful looking after. They need to be almost 
immersed in oil to operate. Parts wear out and 
clutches need adjusting, springs need pulling or 
cutting; above all, TTYs must be kept clean. 
Although great care has been taken not to mention 
actual sums of money in this paper, it is important 
to know that a TTY will cost from $120 to $300/year 
to maintain, depending upon how many machines a 
company is having maintained on a service contract. 
The cost will also depend upon the part of the 
country in which a user is located. 


a. Spare parts . What parts of the system will be 
availabiV irom^ the vendor at a moment's notice? 
Will it be less expensive to obtain parts from the 
manufacturer rather than the vendor*^ (The vendor 
may also be the manufactiu er. ) What about subsys- 
tem parts such as power supplies, converters, 
TTYs, etc ? 


It probably will be worthwhile to stock some lo^’ic 
cards for the computer, interface, and the Fourier 
processor, if the system has a Fourier processor 
Some computer manufacturers sell maintenance and 
spar^* parts kits. Of course, this will add to the 
total cost of the control system. Other spare part« 
and supplies should include: 

(1) Ink ribbons for TTY and printer. 

(2) Paper for TTY and printer or hard-copy 
machine. 

(3) Paper tape for TTY and high-speed 
punch. 

(4) Oil and grease for paper tape punch and 
TTY. 
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Personnel responsibilities in general are discussed in later sections. 
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(5) Blank format tapes for high-speed printer. 

(6) Logic probe (very handy for 
troubleshooting). 

(7) Spare panel lamps. 

(8) Tape winder (this device is an absolute 
necessity) for other than fan fold tape. 

b. Diagnostics . If the vendor spends x hours in 
the development of software logic programs to 
meet an application, approximately 3x hours should 
be spent in developing diagnostic programs to 
determine whether, in fact, the application is being 
met. and if not. why not. This is a very important 
point, perhaps the most important point of this 
paper. When the system is not working correctly, 
how does one know what is wrong? Only the effect 
is observable. The cause is hidden. A properly 
organized diagnostic program will allow the system 
to diagnose itself. The computer manufacturers 
have diagnostics for their computers and their 
computer options. They have expended much money 
and effort in this type of software development. It 
is necessary for the vendor to do likewise. If the 
vendor refuses to develop good diagnostic pro- 
grams. the user must write his own, in which case 
he would require considerable knowledge of the 
vendor's hardware and software. The responsi- 
bility of writing diagnostics should fall on the 
vendor. It is for this reason that a good control 
system specification must include diagnostic 
requirements. The vendor certainly has the right 
to charge for this very important software 
development. 

Table 11 lists the type of diagnostic programs 
required to maintain a digitally controlled random 
test system. 


Table 11. Recommended diagnostics, computer- 
controlled random test system 


Computer^ 

Core memory test 
Central processing unit (CPU) test 
Paper tape reader and punch 
integrity 

Real-time clock test and 
calibration 

Line printer integrity 

System 

Algorithm integrity (dual channel) 
Internal data transfer test 
Core memory test (Fourier 
processor) 

Frame size selection 
Data transfer across interface. 
I/O test 

Control functions test 
Interrupt handling functions test 


Timing'^ 


^Usually supplied by computer manufacturer. 
^Requires logic hardware and oscilloscope. 


It is also very desirable to write some sort of 
timing diagnostic. It seems there is no good way of 
doing this without the help of some logic hardware. 
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Logic hardware should be included within the inter- 
face controller in such a way that* with the aid of 
an oscilloscope and a software logic program, the 
system timing can be mapped and diagnosed. 

E. Facility Requirements 

Some facility changes may be required in order to 
accommodate a computer -controlled lest system in 
a testing laboratory. The test system should not 
be located close to any strong electromagnetic 
interference source such as power ampliUers, 
shakers, shock machines, or cycling vacuum 
pumps. In addition, the equipment should be 
located in an air-conditioned environment where 
there is a minimal change in temperature through- 
out the testing day. It should be a clean environ- 
ment; the user should inspect his own plant's 
computer facility and observe that environment. 

These control systems will require input power 
from 3000 to 6500 W. which means a minimal cur- 
rent requirement of 25 A. The power source should 
be dedicated to the control system. That is. a 
separate transformer should be installed, and the 
control system should be the only equipment con- 
nected to the secondary of that transformer. Mul- 
tiple ac outlets must be provided from the power 
source. The dedicated transformer need not be of 
the isolated type. 

The line voltage. 117 Vac. must be constant, or 
computer problems will result. If the line voltage 
drops too low, the system computer will not operate 
properly, although no one may be immediately 
aware that it is malfunctioning. 

If the control system is to be operated from a 
remote location, a cable raceway will be required 
from the test site to the control system. Only sig- 
nal lines should be laid on this raceway. If power 
switching will be performed from the control sys- 
tem. the power switching lines should not be routed 
parallel to the signal lines. 

Television monitoring equipment and some means 
of audio communications may be required for 
remote operation. 

It may be necessary to reconfigure some existing 
equipment in order to accommodate the control 
system. All these costs must be considered when 
planning for the control system. If the control 
system is to be coupled to an existing computer 
facility, additional hardware including data sets 
(modems) and cabling will be required, increasing 
the total cost. 

F. Personnel Requirements 

The use of a computer in a testing laboratory will 
require new or retrained personnel if the manager 
is to avoid a drastically increasing operating bud- 
get; it is not financially feasible to have the vendor 
in residence at the laboratory. 

Moreover, consider the differences between analog 
and digital test systems: the analog systems con- 
sist of hardware performing a finite number of 
functions, while the digital systems consist of both 
hardware and software and can best be described 


as concepts. If necessary* one can change, modify, 
or add to the software of a digital system in other 
words, change the concept of the system. This 
cannot be done with the analog systems currently 
available. To make the best use of this adaptive 
capability, the manager requires personnel skilled 
in the applications of computer control and knowl- 
edgeable in the specific kinds of testing to be done. 
At least part of the testing organization will be 
fashioned into a small computer center. 

Figure 48 illustrates the relationship between an 
operating system and the skills required for an 
induced environment test group. These are working 
group functions required for performing induced 
environmental testing and in no way reflect the 
structure of the entire test laboratory. (Of course, 
one person may perform two or three different 
functions.) Some of these necessary personnel 
functions are further discussed below. 

1. Data Manager . The data manager (who is not 
the test laboratory manager) is in charge of and 
responsible for the functional structure illustrated 
in Fig. 48. He must have a keen understanding of 
both induced environmental testing and the computer 
sciences as applied to his operating system. His 
background should be mechanical as well as elec- 
trical (electronic) and must include a comprehen- 
sive understanding of the stochastic processes. He 
must also have a good understanding of test 
operations. 

2. Programmer . A programmer is required to 
change and add software logic programs as 
required. The vendor cannot be expected to make 
minor programming changes or write and assemble 
programs without a handsome charge. 

All the programs required by the operating system 
are written in assembly language or other special 
language, not in a symbolic user-oriented language 
such as FORTRAN. The computer manufacturer 
supplies an assembler with his machine to the 
vendor. In addition, a "complete” library of pro- 
grams written by the computer manufacturer is also 
supplied, usually free of charge, to the vendor. The 
vendor, whose task is to put a system together, 
uses some of these programs and some of his own 
to make up the tapes required for the operating 
system. Figure 47 shows the procedure required 
to make a system tape. It is the programmer's 
job to make system tapes, as required. 

If the programmer is to implement the laboratory's 
special requirements, he will surely need records 
of all the venae "'s and computer manufacturer's 
source listings. However, the vendors are reluc- 
tant to supply listings with their systems because 
they believe their ideas should be proprietary. (The 
programmer could, with a good deal of trouble, 
decode the systems tape anyway.) The point is that 
when a computer-controlled test system is speci- 
fied, the user should demand listings of all the 
source programs used in making up the operating 
systems tapes. 

The programmer should know assembly language 
programming as well as those user-oriented lan- 
guages that can be supported by the system 


computer. He should be extremely knowledgeable 
in the physical processes he is trying to imple- 
ment: in this case, vibration, shock, and acoustic 
testing and analyses. He must have such a clear 
understanding of the operating system that he can 
change its concept if required. 

3. System Analyst . The function of the system 
anal^t is to analyze and solve problems ir the 
hardware and to help the programmer debug and 
analyze the results of his programs, /ill pro- 
gramming modifications and additions must be 
approved and recorded by him. The system ana- 
lyst must be kept aware of every hardware and soft- 
ware change. It is his responsibility to record and 
log all system anomalies and his decision on 
whether vendor or computer manufacturer help is 
required. He controls the required inventory on 
spare parts, including logic cards, lamps, TTY 
ribbons and paper, line printer ribbons and paper, 
paper tape, and whatever else is required. He is 
in charge of a preventive maintenance plan for the 
entire test system. He is responsible for keeping 
all software logic program documentation updated, 
including vendor and computer manufacturer soft- 
ware updates and engineering field modifications. 
His educational background should be in the com- 
puter sciences and he must have a good working 
knowledge of the processes implemented by his 
operating system. 

4. System Operator. The operator is responsible 
for operation of the computer system and for seeing 
that the system Implements the required concepts. 
He inputs the test specifications via the conversa- 
tional language progranr^s or prepunched tapes. He 
is responsible for implementing the test specifica- 
tions and obtaining test results for the test engi- 
neer. He notifies the sy&tem analyst in case of 
trouble and works with the system analyst and the 
programmer. He is responsible for running the 
routine system integrity diagnostics executed prior 
to each day's testing. His educational background 
should include some assembly level programming 
and induced environmental te^t operations. 

5. Test Engineer . The test engineer needs addi- 
tional knowledge in the computer sciences and 
considerable knowledge in the field of time series 
analysis. He must be able to determine require- 
ments, procedures, and processes to implement an 
environmental specification and accurately analyze 
the results. 

Training . Training personnel to understand and 
operate a computer-controlled test system should 
be begun before taking delivery of the system. The 
training can come from several sources, including 
the system manufacturer (the vendor), the subsys- 
tem manufacturer, an out-of-plant educational 
facility or organization, or in-plant training facil- 
ity, The results of the training (there should never 
be a final results) should provide answers to the 
following questions: 

(1) What is the system supposed to do? 

(2) Is the system doing what it is supposedto do? 

(3) How is the system doing what it is supposed 
to do? 
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(4) Why 18 the system doing what it is doing? 

(5) What can be done if the system is not doing 
what it is supposed to do? 

(6) How can the system be made to do some- 
thing it is not supposed to do? 

The best way to obtain the answers to the questions 
above is obviously to design and build one's own 
control system, using available subsystems as 
building blocks. The next best way to obtain the 
answers is to become knowledgeable in the follow- 
ing disciplines: 

(1) Environmental testing and current test 
philosophies. 

(2) Engineering mechanics (structural dynam- 
ics), “jhysics, and vacuum technology. 

(3) Electronic engineering: analog and digital 
design and analysis. 

(4) Computer system architecture and 
organization. 

(5) Control system analysis and synthesis. 

(6) Digital communications and information 
theory. 

(7) Assembly language programming. 

(8) Time series analysis. 

(9) Probability theory and statistics. 

(10) Data management. 

The amount of training required to meet the orga- 
nization's needs will depend upon whether existing 
test personnel are trained or new people are 
brought into the test organization. If existing test 
personnel are trained, they will probably be versed 
in the disciplines of items 1 and perhaps 2 only, 
requiring some training in the remaining disci- 
plines. On the other hand, bringing new people 
into the test organization may require less training, 
provided these people have the proper educational 
background and/or work experience. For example, 
a new computer science (or information science) 
graduate will probably require training in the disci- 
plines listed from items 1, 2, and 5. A recent 
electronics engineering graduate should have a good 
working knowledge of most of the abov^' -listed 
disciplines except perhaps those covered in items 1, 
7, and 10. If new personnel must be added, it is 
desirable to hire recent graduates from the appro- 
priate sciences and then train these people as 
required, because the technology in this country is 
probably still doubling evtry four years and, in 
digital electronics and computer field, every two 
years. Hecent graduates provide updated knowl- 
edge, which can be advantageous to any 
organization. 


It is the test laboratory manager's responsibility 
to see tha^ people are qualified and kept prop- 
erly train. perform their work functions. Edu- 
cation is a er-ending requirement in any techni- 
cal discipline, including environmental testing. 
Environmental testing is a philosophy, not a 
science but the means by which the philosophy is 
implemented is a science and requires scientihe 
procedure. 


The bulk of the training burden should be placed on 
the vendor, as it is in the computer business. 

Some training can be obtained from the computer 
manufacturer; this would include items 4 and 6. 

The Fourier processor manufacturer should pro- 
vide the educational requirement to satisfy item 8. 
The remaining educational training will have to be 
provided by either in-plant or out-of-plant educa- 
tional facilities. Out-of-plant facilities include 
local colleges and universities. In addition, there 
are some organizations whose sole function is to 
provide computer science education and training^ 


Environmental science journals and organizations 
have an equal responsibility to provide their 
readers and members with the opportunity to 
become knowledgeable in the items listed above. 
This can be accomplished by inviting papers on 
appropriate topics fundamental to digital control 
systems and by organizing and implementing 
tutorial and workshop lectures at local and national 
meetings. 


G. Costs 

1. Hardware Acquisition Costs . Figure 49 illus- 
trates the hardware cost factor range versus the 
hierarchy of applications defined in Table 7, based 
upon vendor- supplied information. The basic sys- 
tem (group I, Table 7) has a cost factor of I.O units. 

To add group II capability, it will cost from 2.0 to 
2.6 units. That is, if the basic system costs $1 00, 
the group II system, which is a random system, 
will cost from $2.00 to 2.60. Having group II 
capability and desiring group 111 capability, it will 
cost approximately an additional $0.03 to 0.05, for 
a total cost for groups I, II, and III of about $2. 63. 
The cost function, as defined here, is a cumulative 
function starting from the required basic control 
system. No other information is implied from 
Fig. 49. That is, it is impossible to obtain the cost 
factor for going from a sine system to a random 
system from Fig. 49. For exrmple, this figure 
shows a cost factor of about 0.5 going from a 
random system (having the basic system) to a sine 
system. But the author has determined that a cost 
factor of from 1.1 to 1.2 exists going from a com- 
plete n-channel sine system to an n-channel random 
system. 

The important point about the Fig. 49 cumulative cost 
function is that it starts to flatten out above the 
group IV application. The steepest part of the 
curve is between group I and group II. To encom- 
pass all the applications in Fig. 49, group VII will 
require less than twice the money required to get 


^^Two such organizations are the Institute for Advanced Technology (Control Data Corp.)^ 5272 River Hoad, 
Washington, D.C. 20016, and RCA Institutes, 132 W 3St 31st Street, New York, N. Y. 10001. 



40 


to group II. The curve flares up at the end 
(group VII application) because the vendors are 
presently quoting fixed-head disc prices rather 
than tape cassette or cartridge prices. Again, the 
continuous function shown crossing the cost factor 
ranges does not really exist but is drawn in to illus 
trate the cost trend as the number of consecutive 
applications increases. 

The cost factors above pertain to hardware onl> . 

The total costs must rehect operation, mainte- 
nance, and risk, which may outweigh the capital 
hardware cost over some reasonable time period. 

Depreciation (useful life) costs may depend upon 
the particular application and should be considered 
in the gross cost figure. 

2. System Operating Costs . System operating 
costs are denned in this paper as five parameters: 
manpower, maintenance, training, supplies, and 
power. The most significant factor is manpower; 
the least significant cost factor is power. If the 
test system is to be supported by the ideal organi- 
zation previously discussed, there will be at least 
two new possible jobs (therefore, two additional 
people) over and above the personnel requirements 
for an analog test svstem. Both these new people 
(system analyst and programmer) could be in the 
same pay category as a senior engineer. 

Cost for power to operate the control system is 
negligible when considering the overall total 
operating cost. If the system consumes 4 kW, its 
cost per hour will be about $0.08, and, based upon 
normal operating time, the system should cost less 
than $200/year for electricity, although it could 
cost up to $700 /year if the system is left on all the 
time, as it probably should. 

Maintenance costs (that is, labor to maintain) can 
be estimated initially by using 1% of the total hard- 
ware cost of the system as a monthly figure. This 
is a reasonable rule of thumb. Maintenance con- 
tracts for minicomputers can range from 1% of 
total hardware cost per month up to 3.5%. Better 
estimates for new contract maintenance can be 
made if the vendor is willing to reveal his service 
call statistics for existing control systems. A 
computer service call from the manufacturer can 
be as high as $30/hour plus traveling time. 

Training costs can range from $450 for a one- or 
iwo-week class, plus transportation and subsist- 
ence, to free classes with only transportation and 
subsistence costs. 

Supplies to be considered when estimating operating 
costs include such items as ink ribbons, TTY 
paper, printer paper, copy paper, and paper tape. 
Supplies should also include spare parts, including 
spare subsystems if required. 

H. Customer- Vendor Relatio s 

I. I» elivery and Checkout of System. One of the 
current proolems is that, up to now, each delivered 
computer- controlled test system is in some manner 
different from the preceding systems, thus making 
the term "off the shelf* not quite applicable. There 


is certainly nothing wrong with this. What it means 
is that the user is 8p*2Cifying a system for his 
requirements. The vendor is accommodating the 
user. The degree to which the specified system 
(both software and hardware) departs from what has 
previously been delivered will determine the 
delivery time after purchase. 

The first commercially built computer-controlled 
test system took well o\ ‘T a year to deliver. Test 
systems now are being delivered in 40-90 days. 
However, if many "specials" are specified, a 
longer delivery period should be expected. 

Once the system is delivered to the user, a certain 
amount of time will be required by the vendor to 
get the system set up and operating. The amount 
of time required will, again, depend upon the 
special features and upon the vendor's diagnostic 
and tester evolution. Assuming that all hardware 
and software have been proved at the factory, 
several weeks may be required to get the system 
set up and operating. The user should realize that 
large computer systems and data communications 
systems require months of installation and check- 
out. The computer-controlled test system is not 
-s complex as these systems. However, the test 
system is complex enough so that more than one or 
two days will be required for vendor checkout. The 
checkout performed by the user may take weeks and 
will depend upon the prior training and on-site 
instruction received by the system operator, the 
system analyst, and the programmer. 

2. Division of Responsibility . It is the user's 
obligation to determine th^ division of responsibil- 
ity between himself and the vendor, and it is the 
vendor's obligation to inform the user of exactly 
what he, the vendor, is responsible for, prior to 
consummating the procurement. The division of 
responsibility must be stated in the procurement 
specification, and, even then, questions will come 
up, possibly during the procurement and surely 
after delivery of the test system. Again, program 
algorithms and listings fall in this category, to 
name only one important potential problem. A 
good deal of thought and effort on the part of the 
user must go into specifying the systeni require- 
ments. On the other hand, the vendor is obliged to 
answer the procurement specification truthfully, 
listing all the system constraints so that there will 
be no misunderstanding between the user and the 
vendor as to what the system can and cannot do 
when delivered. 
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Fig. 1, Analog vibration test system 

Fig. Z, Computer controlled vibration test 
system 
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Fig. 3. Analog sine scr'^o (conceptual) 

















Fig. 4. Analog random noise servo (conceptual) 
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Fig. 5. Input/output flip-flop characteristics 
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Fig. b. Binary random noise generator (conceptual) 




(C) PROeABILlTY DENSITY FUNCTION 
OF FtJEREO BINARY PULSES 

Fig. 7. Spectrur'i and amplitude distribution 
characteristics of binary random noise 
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Fig. 8. Four-stage shift register (conceptual) 





















Fig, 12. Pseudorandom noise generator, schematic diagram 
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Fig. 16. Data channel of input signal conditioning subsystem 
(conceptual) 
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locations 
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Fig. 2^, Quantization effect and 
the resulting error function 
(Ref. 15) 



Fig, 2*1. Effect of sampling at 
inc-' equate rate and, thus, 
creating "aliases” (Ref. 15) 
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Fig. 23. Continuously varying signal (a) and (b) and 
random noise (c) and (d) displayed as amplitude 
function of time (left) and as frequency spectra 
(right) (Ref. 15) 



Fig. 25. Effect of multiplying continuously varying 
signal b> flat-topped sampling function (Ref. 15) 
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^ig. 26. Relationship between time function and 
frequency function for dat- topped repetitive 
pulses (Ref. 15) 
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Fig. 30. Hanning data weighting 
(normalized) 
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Fig. 31. Flowchart for power-of 
two frame size test 
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Fig. 32. Flowchart for input shuffling process 


Fig. 33. Flowchart FFT/IFT algorithm 
































NOTE ^ DENOTES MULTIPLICATtON 

Fig. 34. Power spectral density- 
process 



Fig. 35. Cross spectrum process 
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Fig, 38 . Random servo control (APS method) functional block diagram 
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Fig. 44. Transform pair of the TWC initial 
calibration pulse 


Fig, 45a. Wavelet with 7 half cycles 
and 60 Hertz fundamental frequency 

Fig. 4 5b. Wavelet with 9 half cycles 
and 120 Hertz fundamental frequency, 
and a delay of 10 millisecondb 
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Fig, 46. Composite waveform from 
typical shock spectrum specification 
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Fig, 47. System program 
generator procedure 













Fig. 48, Operating system and personnel requirements 



Fig. 49. Cost factor range versus hierarchy of 
application 
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APPENDIX 

COMPUTERS FOR ENVIRONMENTAL TEST SYSTEMS^ 2 


I. COMPUTER TECHNOLOGY 

Modern -r'ay digital computation started in the 
mid- 1 94 J' 6 with the programmable vacuum tube 
computers developed for the military for purposes 
of computing ballistic firing tables. Figure A-1 
indicates the time chart of computer technology in 
the United States. Note that computers for envirou > 
mental test control were utilized starting in 1968. 

Computer technology today is doubling perhaps 
every four years. Computer hardware costs are 
being reduced by an order of magnitude every six 
years. These costs are linked to technology and 
technological advances associated with batch- 
fal ricated microcircuits. However, peripheral 
hardware, i. e. • that hardware hung onto the com- 
puters such as terminals, tape drives, printers, 
etc. . is not dropping in price as rapidly as the 
computer costs. This is because these peripheral 
devices are electromechanical in nature. Software 
costs, on the other hand, are increasmg in that 
software development is tied directly with people's 
salaries. 

It can be seen from Fig. A-1 that the advent of the 
minicomputer took place about 1964. With the 
growth of the minicomputer market and the heavy 
inroads the minicomputer has made into applica- 
tions areas, the mini is rapidly replacing many of 
the larger, more costly machines, despite the fact 
that the applications software is generally a more 
expensive development, and that the mini, by 
dehnition. cannot compete in terms of word and 
memory size with larger, more powerful comput- 


II. COMPUTER APPLICATIONS 

Figure A-Z shows the development of computer 
applications as computer technology advanced. To 
the application of computing ballistic trajectories 
was added the computation tasks of simple batch 
processing of business and scientific jobs. In turn, 
applications utilizing these machines were tried 
and developed at a rather rapid rate as the technol- 
ogy advanced even further. Operating systems 
were configured to implement sequential batch 
processing of data, asynchronous multiprogram- 
ming requirements, and input/output oriented 
operations. Additional applications stemmed from 
these applications in turn. The following applica- 
tion areas for present day computer systems are 
briefly defined: 

Process Control. This application includes 
machine tool control, production line control, 
oil and gas field monitoring, biomedical 
monitoring, and laboratory monitoring. 
Scientific Computing. This includes equation 
solving, engineering and system design, 


system analysis and simulation, and system 
modeling. 

Business Data Processing. Here payrolls are 
g ene rateci. account s paya ole are paid, accounts 
receivable are recorded, inventories are 
controlled, and in summary, the computer is 
utilized as a valuable management tool. This 
includes the very powerful techniques of linear 
programming, probability theory, reliability, 
nonlinear programming, statistics, decision 
theory, queue and game theory, and informa- 
tion theory. 

Military/Spacecraft Command and Control. 
These appUcVtions include guidance and track - 
ing. weapons control, navigation, communica- 
tions control, and cryptographic processing. 
Real Time Processing. This includes environ- 
mental test control, which is a subset of 
command and control applications, brokerage 
systems, and airline reservation systems, to 
name a few. 

Time Sharing. These applications include 
remote inquiry, scientific computing, and data 
base access. 

Data Communications Control. This is 
probably the newest application field, which 
includes message switching, telemetry, 
communications concentrators, and. most 
importantly, front ends for the "super com- 
puters". 

Prep r oc e 8 8 ing. This application includes 
buffering, peApheral control, display control, 
and communications con^-rol. 

The point here is that computer utilization for 
environmental test control is only a very limited 
application in the to^al hierarchy of applications 
to date. And still, the computer used for test 
control can and should certainly be utilized for 
other than real time processing within the test 
lab. 

III. THi: SPECTRUM OF POWER 

Computers can be categorized in terms of CPU 
(Central Processing Unit) capability, word length, 
memory capacity, input/output (I/O) capability, 
complexity of interfacing, and supplied software. 
Cost is another factor. Figure A- 3 defines the 
spectrum of power of computers as a function of 
word size and cost from the mini to the super 
number crunchers. The mini is. of course, the 
most popular type of computer found in the 
environmental labs these days, although not 
exclusively. 

Minis and other computers are characterized by 
I/O throughput speeds of up to I million words per 
second. (A computer word will be defined 
shortly. ) In addition, typical instruction speeds 
of most computers range from 100 to 600 thousand 


Some of the figures in this appendix were redrawn from Kef. A-1. 
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instructions per second (KIPS), The truly large 
multiprocessors and supercomputers range from 
1000 KIPS ov million-instructions per second 
(MIPS) up to 10 or 20 MIPS. Since the mini is of 
prime importance in this paper, it should be noted 
that the mini here is defined as a computer whose 
cost is less than $50, 000, the word sire is not 
greater than 18 bits, its core or dynamic memory 
is no greater than 32,000 words, and its physical 
size is small, and in fact, will generally fit in a 
19-inch equipment rack. In addition, it requires 
less than 30-A service at standard 117 Vac line 
voltage, and requires no special environment or 
air conditioning. 

IV. DEFINITIONS AND JARGON 

In general, minicomputers have the same basic 
elements found in their larger counterparts. The 
mini system elements fall into the following 
subsystems: 

Processor (CPU). Most minis employ a single 
address structure, utilize two’s complement 
negative numbers, and generally have no 
greater than 8 general purpose hardware 
registers. Usually, hardware multiply and 
divide is offered as an option. The minis 
usually have an interrupt handling capability 
(to be discussed) and a bus arrangement that 
allows transfer of data and control signals. 
Memory. Most minis in the past were avail- 
aible with, only core memory. Memory cycle 
times ranged from 600 nanoseconds to 1 micro- 
second. The new minis are now being made 
available with solid state dynamic memories 
whose cycle times range in the low hundreds of 
nanoseconds. The dynamic memories are 
faster than core, but more expensive at the 
present than core, and volatile in contrast to 
the nondestructive cores. In either case, 
however, memory blocks come in increments 
of 1024, 4096, or 8192 words. Most minis 
have a memory protect scheme that prevents 
the user from wiping out certain sacred areas 
of memory. In addition, a parity check is 
utilized in order to provide confidence of the 
integrity of the data retrieved from or 
deposited to memory. 

Input/Output. Data is transferred to or from 
a peripheral device by one or more of three 
popular techniques: Processor Program 
controlled I/O through hardware registers in 
the CPU; or Direct Memory Access (DMA); 
or Direct Multiplexed Memory Access (DMC). 
The DMA technique of data transfer is the 
most efficient and fastest scheme in that it 
allows the CPU to function while data transfers 
occur. 

Interrupt Structure. The interrupt structure 
concept is perhaps the most important parame- 
ter for minis used in real time control 
processes such as environmental test control. 

A computer interrupt is a means by which a 
signal occurs either internally to the computer 
or from a peripheral device that directs the 
computer to execute a special sequence of 
instructions as a function of the alerting 
signal. Priorities can be structured such that 


a signal from one device takes preference 
over a signal from a lower priority device. 
Figure A-4 is a flow diagram indicating the 
concepts of priority interrupts. 

Software. Software for minis falls into four 
categories. The first category is program 
development software that is required by the 
user to develop his programs for particular 
applications. This includes editors, assem- 
blers, debuggers, and utility routines and 
compilers such as BASIC or FORTRAN, 
Another category is Input/Output software 
for the system hardware and peripherals. 

This is the software that makes the termi- 
nals, line printers, reader, and punches 
work as well as the mass storage devices. 

A third category of software is the applica- 
tion programs. These programs are related 
to the task that the system is to perform and 
which is therefore unique to the particular 
system. The last category includes operat- 
ing system software, also called the execu- 
tive or system monitor, which tells the user 
what to do, when to do it, and what to do it 
with or to. 

Minicomputer Peripherals. Peripherals for 
minicomputers can be separated into two 
broad categories - input devices and output 
devices. Most peripherals are used for both 
input and o\itput functions. Input devices 
include keyboard terminals such as teletype 
(TTY) machines and Cathode Ray Tube 
(CRT) terminals, paper tape readers, fixed 
head disks, card readers and magnetic tape 
units including cassettes. Output devices 
include TTY's, JRT’s, card and paper tape 
punches, dis^' , line printers, magnetic tape 
units and an X-Y recorder (plotter). For 
environmental test control, the following are 
the most popular peripherals: 

(1) TTY or preferably a CRT terminal 
with hard copy capability. 

(2) High-speed paper tape reader and 
punch, 

(3) Low-cost disk (i. e. , the floppy disk) 
or cassette tape. It is highly desirable 
to have a line printer if your organiza- 
tion is going to assemble programs. 

V. MINICOMPUTER SIZE, EXPAN^JABILITY 
AND FLEXIBILITY 

Most minis fit in 19 -inch equipment racks. 

Optional expansion chassis are available for sys- 
tem expansion. These expansion chassis include 
power supplies and cooling fans and can be used 
to house additional core or dynamic memory and 
additional logic required for controlling line 
printers, disks, and other added peripherals. 

Many minis have a memory excnangeability 
feature that allows the faster dynamic semicon- 
ductor memories to replace the slower core 
memories by a board by board replacement (no 
mixing allowed). As faster memories become 
available, the old memory can simply be 
removed and replaced by a newer, faster one. 

With denser packaging technology, it is possible 
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to replace a 4000-vvord circuit board by a single 
8000-word board and later, no doubt, by a single 
16, 000-word board. Special interfacing can be 
housed in the expansion chassis. 

VI, DEFINITIONS 


either peripheral or remote terminal commun*i^a- 
tions control systems. Figures A-7, A-8, A-9, 
and A- 10 illustrate both systems. Note that the 
prc-. ess control application requiies a feedback 
configuration. This is a very important concept 
in environmental test control. 


The following definitions are listed here in an 

attempt to clarify some basic computer tei‘ms: 

Bit. A bit is a binary digit, either 0 or 1, 
Baud. This is a unit of rate ox infoimation. 

In general (but not always), a Baud rate can 
be divided by 1 0 to give an approximate 
character rate in units of characters per 
second (CPS). For example, 1600 Baud is 
about 160 CPS if the characters are ASCII 
characters. 

ASCII, This stands for American Standard 
Code for Information Interchange. This is an 
attempt to standardize certain keyboard 
characters in terms of octal numbers. 
Character. A character a symbol that is 
defined under the ASCII standard. These 
characters are generated from keyboards, 
CRT terminals and printers. They include 
letters, numbers, punctuation, and special 
symbols including space, carriage return, 
line feed, etc. 

Word, A computer word is a binary number 
which represents either data, an instruction, 
or an address. Minicomputer words are 
generally 8 bits, IZ bit'-, 16 bits, or 18 bits 
long. The bits can be grouped in three's and 
read as an octal number. 

Byte. A byte is an 8-bit code. The code 
represents an alpha-numeric character or 
word. Generally two or more bytes repre- 
sent a word, although by defiiiition, Z bytes 
can represent 2 words. Not all minicomput- 
ers are byte oriented. 

VU. NUMBER SYSTEMS FOR COMPUTERS 


IX. INDUSTRY STANDARDS 

Industry-wide standards for computers are 
either ill-defined or nonexistent. In most cases, 
the computer giant, IBM, h^ created its own 
standards, which to a great extent have, in turn, 
been adopted by the competition. These same 
IBM standards have also been adopted by the 
minicomputer industry. However, these trans- 
formed standards are generally for peripherals, 
not for cunjputers. For example, transmission 
rate agreements are; 


(U 

Low speed: 

1 1 0 to 300 bits per 



second. 

(2) 

Medium speed: 

300 to 2400 bits per 



second. 

(3) 

High speed: 

up to 50, 000 bits per 


second. 

Other standards (more properly, agreements) are 
132-column line printers where speed is measured 
in lines per minute (LPM). Medium speeds are 
around 300 to 600 EPM, and high speeds are 1100 
to 2000 LPM. The 80-column IBM Hollerith 
Punched Card is still one of the most widely used 
l/O media in the computer industry, although paper 
tape devices, being at the low end in terms of cost 
of l/O devices, are very popular in the minicom- 
puter industry. Medium- speed paper tape punching 
is in the range of 50 to 75 characters per second, 
while the medium-speed paper tape reader range 
is from 100 to 300 characters per second, 

X. MINICOMPUTER INDUSTRY 


Figure A-5 lists equivalent numbers in four 
number systems. Most mini programmers prefer 
to work in the octal numbering system since the 
binary bits can be grouped in three's and the octal 
number read directly fro*, i the grouping. Negative 
numbers are generated within the mini as two's 
complement. This is done by the machine by 
first one's complementing the bits (i. e, , 0 
becomes 1 and I becomes 0), then adding I to the 
results. An example is shown in Fig, A-6. The 
advpntage of two's cornplement is simply that 
binary subtraction results in a sum process and, 
after all, adding is the thing a computer can do 
most efficiently. 

Via the MINICOMPUTER AS A SYSTEMS 
ELEMENT 

The mini should be viewed by Ihe systems 
designer as a small but nevertheless important 
element in the cotal system. The most common 
minicomputer systems are cither man-machine 
or machine-ma ;hine oriented. The former con- 
sist of data acquisition, process control including 
environmental test control, and time sharing and 
problem solving systems, while the latter are 


There are four main segments of the minicomputer 
industry - the main-frame manufacturers, peri- 
pheral manufacturers, soxtware suppliers, and 
turn-key system suppliers, 

A. Main- Frame Manufacturers 


The minicomputer manufacturers were an out- 
growth of companies supplying circuit boards to 
other manufacturers. In 1970, there were 40 to 
50 companies manufacturing main frames Now, 
in early 1 975, there are less than a dozen minicom- 
puter main-frame manufacturers. A few of these 
companies are listed below: 


Manufacturer 

Honeywell 

Hewlett- 

Packard 

Digital Equip- 
ment Corp. 

Interdata 


Computer 


H 316, 516, 
716 

2100 A, 2114, 
2115, 2116 

DEC POP 8, 
POP 1 I 

Model 70 



Manufacturer 


Computer 


A, Memory 


Varian Data 620, 73 

Machin 

Data General NOVA 

B. Peripheral Manufacturers 


The second largest segment of the mini market 
is the peripheral manufacturers. Some of the 
mainframe manufacturers also make their own 
peripherals. The "peripheral" manufacturer buys 
or builds the end product and develops his own 
hardware/ software controller. He then sells the 
entire system to the end user. His system may be 
an improvement on the comparable system sold by 
the main-frame manufacturer in that his peripheral 
control software may be more efficient, thus 
requiring less main- frame memory, or perhaps 
his system sells for substantially less. Periph- 
erals include line printers, card readers, paper 
tape readers and punches, disk drives, plotters, 
CRT terminals and other keyboard terminals. 

C. Software Suppliers 

The third segment of the market consists of the 
software supplier .’'ho may be providing a more 
efficient, higher-level language oa* a general 
purpose modularized applications pa^-kage oan 
be run on a widely used mini. Many such packages 
exist for the DEC PDP 8 and PDF 1 1 as well as for 
the Data General NOVA. Some of the larger 
companies supplying software to end users of 
minis are Computer Sciences, Planning Research, 
Syster.' Development Corporation, and Informatics. 

D. Turn- Key System Suppliers 

The fourth segment of the marke. comprises the 
independent turn-key system suppliers who 
compete head-on with the main- frame manufac- 
turers, This supplier is generally at somewhat of 
a disadvantage tince he must purchase a large 
part, if not all, cf the hardware from the hardware 
manufacturer. The "value added" in thir case is 
the applications software and th.j intimate knowl- 
edge of the applications area. Herein lie, at least 
to some extent, the suppliers of environmental test 
systems. 


The memory unit provides the main storage for 
actual data as well as instructions that tell the 
control unit what to do with the data. Figure A- 12 
defines the memory organization concept. Perhaps 
the most difficult concept for a newcomer to com- 
puters to grasp is the difference between memory 
addresses (location and the contents of those 
addresses. 

3. Arithmetic and Logic Processor 

This subunit temporarily stores data received 
from memory and performs calculations and logic 
operations on this data. The arithmetic/logic 
processor contains one or more registers which, 
in turn, contain the data being operated on as 
defined by an instruction. 

C. Control Unit 

This logic controls the flow of data in the system, 
fetches instructions from memory, and decodes 
the instructions in one or more instruction regis- 
ters. This unit executes instructions by enabling 
the appropriate electronic signal paths and control- 
ling the proper sequence of operations performed 
by the arithmetic/logic and I/O units. Finally, 
it changes the state of the computer to that 
required by the next operation. 

D. Input/Output Unit 

This unit provides the interface and in some 
systems, buffering to peripherals connected to 
the computer, transferring data to and receiving 
data from the outside world. 

It should be noted that some of the more recently 
introduced minicomputers treat memory as an 
external device and communicate with it over the 
same I/O bus, just as they would with an external 
peripheral device. 

The central processing unit (CPU) is defined as 
the control unit, the arithmetic/logic processor, 
and all interrelated buses and registers. The 
basic elements of the CPU are: 

(1 ^ The Instruction Register(s). 


XL MINICOMPUTER ARCHITECTURE 


^2) ‘The Decoder. 


The architecture of a mini is the key to its 
performance and an indication of its relative utility 
and power. Machine architecture has a direct 
bearing on its operating characteristics, such as 
speed and efficiency, how instructions are carried 
out, and bow easy or difficult it is to program the 
machine. The basic building blocks of all minis 
are its memory, the arithmetic and logic proces- 
sor, the control logic, and the I/O logic. T1 ese 
subunits are interconnected by several lines that 
provide paths for data, control signals, and com- 
puter instructions. These lines are commonly 
called buses. Figure A-11 shows the bus line 
interconnections of the main building blocks of 
most minicomputers. 


(3) The Control Unit. 

(4) The Program Counter (Pointer). 

(5) The Adder and Comparer. 

(6) The AccuiTiula' r. 

(7) The Status Register, 

A typical minicomputer processing cycle can be 
divided into two subcycles - an instruction or 
fetch cycle and an execution cycle. A detailed 
description of each of these elements is beyond 
the scope of this paper (see Ref. A-l). However, 


the basic ideas of the central processing unit can be 
obtained by describing the chain of events that takes 
place during the fetch and execution cycles. During 
the fetch cycle, 

(1 ) The control unit logic fetches next memory 
address n f-:om the program counter. 

(2) The decoder decodes this address. 

(3) The control unit logic fetches ^he contents 
A (an instruction) of memory address n, 

(4) The control unit logic loads instruction A 
into the instruction register. 

(5) The control unit logic increments the pro- 
gram counter by 1 for the next cycle, i.e. , 
the counter acts as a memory pointer. 

(6) The control unit fetches the contents of 
instruction register. A, and decodes the 
instruction. 

The control unit is now ready to implement the 
second subcycle, the execution cycle. The steps fol- 
lowed in the execution cycle depend on the type of 
instruction to be executed. One of the most com- 
mon instructions is addition or ADD. A typical 
ADD instruction is based on the following sequence 
of steps: 


(7) A piece of data designated here as "a", 
stored in memory location m + n, is added 
to the contents of the accumulator. 

(8) From the previous execution cycle, 
assume that the previously stored quantity 
in the accun^ulator is ”b"; the new amount 
after additio. will be “a” + ''b”, 

(9) The new contents of the accumulator are 
stored in location m + n, which will now 
contain “a" + “b". 

The status register varies in size and complexity, 
but as a minimum, two single bit registers are 
used where the first stores the overflow from the 
accumulator, and the second register, called the 
test or link register, tests the condition of the 
accumulator. 

Figure A-13 shows the architecture of a micro- 
processor. This is, for all practical purposes, a 
CPU on a chip, A 40 -pin chip slightly over 2 
inches in length and about 5/8 of an inch wide* 

This should give the reader insight into things to 
come in the computer industry in the near future. 
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Fig. A-l. Evolution of computer technology 
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COMPUTER COST IN THOUSANDS OF $ 


^'ig. A- 3, Classification of computers 
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MINICOMPUTERS 


Stephen A, Kallis, Jr, 

Digifol Equipment Coq^oratton 
Maynard, Massachusetts 

A tutorial paper on minicomputers. Minicomputers share the characteristics 
of other general-purpose digital electronic computers in thot they are programmable, 
can be used in batch, real-time, or timesharing mode. In addition, their degree 
of rugged ness, when compared to larger electronic computers, permits them to be 
employed in more severe environments than the larger machines; their cost is low 
enough to permit them to be employed as instruments in laboratory and test 
si tuations. 


INTRODUCTION 

Since its introduction in 1963 with the PDP-5, 
the minicomputer has had a significont impact in a 
variety of onalyticol fields as well as the more con- 
ventional “computer applications" such as statistical 
analyses, financial calculations, and the like. The 
minicomputers available today have onalyticol power 
comporable to the medium-scale computers of a 
decode ago, yet their costs ore so low that a few 
people have purchased minicomputer systems for pri- 
vate use. Small businesses, and as important, small 
laboratories and research facilities are acquiring mini- 
computers to permit them to extend their activities 
into areas that it would have been impractical for 
them to have attempted only a few years ago. 

A minicomputer is defined as a general-pur- 
pose electronic computer with a digital word length 
of between 8 end 18 bits, a minimum of 4,096 words 
in its memory, and a cost of less than $20,000 in its 
minimum configuration. For clarity, this definition 
will be fined further. 

A “general -purpose" computer is one that can 
be employed for a variety of different applications. 
This is achieved through olterotion in the programs 
stored within the computer's memory. Thus, the same 
computer could be used in connection with laboratory 
instruments at one time, for inventory control at 
another time, and for mathematical analyses at yet a 
third time, merely by reading new progroms into the 
computer's memory for each new application. 

This should be contrasted with a “special- 
purpose" computer, which can only be used for one 
sort of application. The flight computer oboard o 


spacecraft, for instance, is an exomple of a special- 
purpose computer. 

The fundamental unit of information in a 
digital computer is a "bit, “ from Binary digIT. For 
operation in computers, bits ore grouped into 
aggregations known as "words"; common minicom- 
puter words are 12 and 16 bits long. 

The capacity of a computer's memory is 
generally considered the number of words it can 
store. Common values for minicomputer capacity 
include 4,096, 8,192, and 16, 384 words. All 
these numbers are powers of 2. In computer parlance, 
the term "K, " when used for memory copacity makes 
use of the fact that 2 is approximately equal to 
1,000; therefore, the computer industry uses "K" t© 
stand for 1,024 — and by this definition, a computer 
with a memory capacity of 4,096 words would be 
said to have a “4K word" memory, and so forth. 

Sometimes computer professionals work with 
another grouping of bits colled “bytes," This term 
is a bit flexible, but usuolly refers to a group of 
bits less than one computer word in length, and fre- 
quently o submul tiple of that word (e.g., a 12-bit 
word might be divided into two 6- bit bytes). To a 
number of specialists, though, "byte" seems to have 
a more fixed length, usually 8 bits. For discussion 
of the elements and applications of minicomputer 
systems, the term "byte" is not essential, and will 
be clearly defined in any areas where it is employed. 

The cost factor of the definition of the mini- 
computer was strictly arbitrary and stemmed from the 
price of the original commercially available mini- 
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comput*er, the PDP-5. The line of computers de- 
riving from this machine, the PDP-8 famiiy, has 
averaged a central processor cost reduction of opproxi- 
mately 10 percent per year, so the cost of the average 
minimum-configuration minicomputer system is 
opproximately an order of magnitude below that of 
the defined minimum. 

MODES OF OPERATION 

There are two ways a computer can be used: 

A$ Q "reel -time" device or os o “non real-time*' 
device. In computer terminology, "real-time" opera- 
tions are those in which the computer is able to 
assimilate information concerning a process, process 
the data, and present the solution or solutions with 
sufficient speed so that the solution can effect the 
process under analysis. Although this definition 
implies a direct control feedback, intervention 
could take place through the agency of a human 
operotor. The important concept for real-time 
Operations is that »he time spent in analysis should be 
extremely short when contrasted to the process being 
analyzed. 

"Non real-time" refers to an operation that 
does not require Immediate results. Traditional com- 
puter operations such as bookkeeping, financial 
arxilyses, generation of tobies, and the like, ore 
non real-time functions. 

COMPUTER SYSTEM ELEMENTS 

A computer system is composed of various 
physical devices, loosely cotegorized as "hardware, " 
and computer programs, colled "software." Within 
the category of hardware, there is the computer, 
peripheral devices, and interfaces. Within the 
category of software, there are application programs, 
general-purpose high- and low-level programs, and 
operating systems. 

Treating the hardware first, there is the com- 
puter itself, generally called the central processing 
unit, or "CPU." Within this unit, all the electronics 
necessary to make the computer operate will be found. 
In many cases, the computer's memory is considered a 
part of the central processing unit, although in some 
schemes, the memory is treated as a separate element 
of a computer system (see Fig. 1). 

In gross, the orilhmetic ond logical operations 
of the computer ore performed by the central proces- 
sing unit. 

To extend the capabilities of the computer 
system, a number of devices, known as "peripheral 
devices, " or "peripherals, " are attached to the CPU. 
Peripherols fall into two main categories; "Input- 
output" (or I/O) and "storoge" devices. The first 
category includes all units that establish a communi- 
cation chonnel between the computer system and the 


user. Included in Input-output devices are terminals, 
printers, plotters, card readers, paper ^ape readers, 
optical scanners, and paper tape punches. Because 
of the variety of oval table terminals, It is necessory 
to add that termiixrls may take the form of 
writers, supplying the user with o printed record of 
his communication with the computer system; or video 
terminals, where communications are dlsployed on o 
CRT screen. Those I/O devices producing tangible 
copies of computer communications ore known os 
"hard copy" devices; those displaying communications 
are called "soft copy" devices. 

Storage peripherals act to augment the com- 
puter's memory. Dota or programs are stored on a 
recording medium and ore ovailable for callup by 
the computer in its operations. One of the most 
common storage medio is magnetic tope. Magnetic 
tape is ovoiloblc on a number of different reel sizes, 
and is, relative to other media. Inexpensive. Mag- 
netic cassettes are considered os being in the "tope" 
category. 

While magnetic topes ore inexpensive, they 
are also relatively slow in action. Another medium, 
the magnetic disk, operates at many times the speed 
of a tape: Where tape access to files may take 
seconds to achieve, access to those same files if 
stored on a disk can be achieved in 70 milliseconds 
or less (u rider 20 milliseconds for a fixed -head disk). 
This rapid rate of data transfer permits some varieties 
of disk to be used in timesharing opplications, where 
the computer works with a number of terminals 
"simultaneously" — actuolly by a form of time- 
division multiplexing. 

A third magnetic medium used for data or 
program storage is the magnetic drum. These operate 
on the same principle as o magnetic disk, but with 
slightly faster access. Drums are not usually used 
with minicomputers. 

Minicomputer peripherals, like minicomputers 
themselves, are both rugged and inexpensive. They 
tend to be usable In a less controlled environment 
than that required for large-scale computers. When 
tlie minicomputer wos first developed, the only 
available peripherals were those designed for larger- 
scale computer systems, ond the minicomputer system 
designer was frequently required to tic a peripheral 
to the control processor that cost four times as much 
as the CPU did, merely because thot was the only 
unit avoiloble. Todoy, the low cost and durability 
of minicomputer peripherals are such that they are 
being used with large-scale computer systems. 

A final type of hardwore is the interface. 
When a computer is used to gather information from a 
loborotory instrument, or when one is used to control 
a process of some sort, appropriate electronics are 
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Fiaure l.A - MmicomDuter scheme 





B - Mini computer scheme in which the memory is considered separate from the centrol process* un?^ 


required to effect a connection between the computer 
and the process or instrument. In some coses, the 
interface con be a direct digital connection between 
the computer and a switch assembly (or equivalent). 
More frequently, however, some sort of conversion 
process is required to permit the computer to gather 
information from the process or generote instructions 
to it. 

Mony rfeodings from instruments or sensors are 
in the form of electricol magnitudes (analog signals). 
To convert them into a form compatible with the 
circuitry of the minicomputer, an analog-to-digital 
converter is used. To convert information from the 
computer to analog signals required to control some 
process, o digital-to-analog converter is employed. 


For long-distance connection between a terminal 
and o computer system, a special communication 
interface known as a "modem” (for MOdulotor- 
DEModulator) is employed. A modem permits commu- 
nication between a computer and a terminal (or 
between computers — or termlnols) over telephone 
ci rcui ts. 

Programs, or software, fall into several 
categories. Some programs ore very specialized, 
and permit a computer to perform a specific job only. 
These programs are called "application progroms, " 
because of their rwture. There are many and varied 
application progroms, ronging from simple tosks such 
as the generation of mailing lists to very complex 
programs such os those used to perform Fourier 
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analyses on incoming waveforms* Application pro- 
grams encompass o voriety of disciplines^ but each 
generally restricts a minicomputer to o single {ob or a 
single type of job. 

A general-purpose high- or low-level pro- 
gramming longuoge Is used to develop application 
programs. A "high-level" programming language Is 
one like FORTRAN, BASIC, COBOL, or FOCAL®. 
These languages are similar to English, and are rela- 
tively easy to employ. Most such languages are 
"machine independent, " meoning that the program 
written In, say, FORTRAN IV for one machine will 
run on another machine that was built by a different 
manufacturer. A high-level language Is effectively 
Independent of a minicomputer's Internal orgonizatlon, 
as far as the user Is concerned, as Is not the case with 
a low-level language. 

A low-level language is written to reflect the 
actual orgonizotion of the minicomputer and its 
working parts. Such an "assembly language, " requires 
knowledge of the way a particular computer ope rotes 
and is totally dependent upon the computer thot it is 
being run on. Whereas a high-level language Is easy 
to write and requires only a few instructions to execute 
a program, a low-level language requires utter pre- 
cision in writing every step of action required in the 
program. 

The advantage of a low-level language over a 
high-level language is quite simple: Although the 
Instruction list is much longer for a low-level language 
the actual memory that contains the final progrom is 
much smaller than for a high-level language. The 
reason for this is that the computer only operates in 
binary Instructions, and both the low-level language 
and the high-level language hove to be transformed 
Into binary instructions if the computer is to use them. 
In the case of the low-level language, this is o simple 
process; in the case of the high-level language, a 
large and complex progrom to tronslate the English- 
like program into binory instructions must olso be pre- 
sent in memory. This large progrom (called either an 
"inteqjreter" or a "compiler, " depending upon its 
operation), because of its generol purpose nature, 
occupies a large amount of a minicomputer's memory. 
Where memory is ot a premium, the user may be forced 
to consider u*ing a low-level language in order to 
conserve momory space. 

Operoting systems (also known as "executives") 
ore programs designed to direct the operation of other 
program elements in a computer system. An operating 
system moy be used to initiate a data transfer between 
peripherals, call up new programs (or parts of pro- 
groms) for execution, or for library functions (keeping 
track of data files, for example). Operating systems 
have been avoilable for computers for several years, 
and currently ore quite sophlstlcoted. 


MINICOMPUTER USES 

As noted before, minicomputers con be em- 
ployed In a variety of dif rerent oppllcotlons. In 
general, these foil Into the two brood categories of 
dota onolysls and control. Administrative duties and 
data gathering without analysis con be associated 
with the former category; data communication and 
messoge switching, ^e latter. But there are two 
other fundamental ways o computer con be used: As 
a "stand-alone" unit, or as a part of some sort of 
computer hlerorchy. 

The term "stand-alone" Implies the computer 
use: The computer is used as an independent unit, 
tied into no other computer system. The relatively 
advonced software available today mokes such an 
approach practical. Today, minicomputer systems 
can run ANSI-compatible FORTRAN IV programs in 
as little as 8K words of memory, wi^h the assistance 
of the appropriate perlpherol(s) such as a moss storage 
device. This degree of power means that computer 
systems that con fit atop a desk con deliver power to 
the engineer that only a few years ago wos available 
only through use of a large-scale centrally located 
computer system. 

Stand-alone computer systems permit a degree 
of protection for a facility's coerotlons that is diffi- 
cult to achieve by other means. A number of stond- 
alone systems scottered thro'jghout the focility result 
in a distribution of the total computer power among a 
number of independent stations. The workload of one 
system will not necessarily interfere with the opera- 
tion of another system, and In the event of o system 
failure, only the one station will be affected, not 
every station. 

Where a stand-alone computer is used for 
only a froction of a workdoy on a single oppll cation, 
it is frequently possible to re-progrom it for other 
applicotions in Its Idle time. Some facilities hove 
mounted their smaller minicomputer systems on casters 
so that they can be wheeled between Ic oratories, 
permitting them to be shared during the course of a 
workday. 

Minicomputers can also be used in hierarchies. 
A number of minicomputers can be connected, via 
communication lines of some sort, to a lorger com- 
puter. In this sort of arrongement, the minicomputer 
acts in a semi -independent mode, transmitting data 
to the larger computer but handling ordinary prob- 
lems itself. If, however, it encounters on extro- 
ordinory problem, beyond its copocity to solve — or 
beyond its copacity to solve in an occeptoble period 
of time — It may format its problem in o "predigested" 
form and forward the problem to the larger computer 
to solve. The solution is returned to the minicom- 
puter, which then handles the solution a$ If It had 
performed the necessary operations. 
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Figure 2, - A compufer hierarchy. In this example, fhe computer hierarchy resembles the corporate 5 truc^ure 
of 0 manufacturing facility, with each level of computer activity mirroring an equivalent human activity. 


The lorger computer may in turn be tied to an 
even larger computer. The minicomputer can, in turn, 
Ey linked to even simpler autorrwitic devices, such as 
progro:nmable controllers. Such a hierarchy porallels 
the hierarchy found in large organizations. (See 
Figure 2 above) 

While computer hierorchies are inte linked, 
there is also a certain degree of independence be- 
tween the various levels. Minicomputers, for example, 
ordinarily will work independently of the medium- 
scale computer. The only times that communication 
normally will occur between these machines is either 
when the medium-scole computer interrogates the 
minicomputer for information, or the minicomputer 
forwards data to the larger machine to obtoin the 
solution of a problem. The larger machine con for- 
word information to the executive machine; it also 
can coordinate the operation of the minicomputers it 
is monitoring. The effect is the most optimum use of 
computers in o focllity, with each being used at its 
maximum efficiency, with minimum interference from 
each other. Even with a short breakdown in communi- 
cation, the elements of the hierarchy would be able 
to operate independently. 


CONCLUSION 

The minicomputer has evolved into a power- 
ful, rugged, and reliable unit thot con be used as a 
supercalculator, an analytical device that con be 
attoched directly to o process, a control device, an 
administrative device, or an element of o computer 
hierarchy, Historicolly, the cost of minicomputers 
hos diminished while their power has been maintained, 
or has actually irrcreased. Software for minicom- 
puters has become highly sophisticoted, and mini- 
computers con perform operations today that only a 
f years ago were the exclusive province of much 
i<*rcer computers. The result is a tool that can be 
• I'c: rporated into o wide variety of systems for o 
lc-»;e number of varied applications. 


Note: •*PDP" and ’'FOCAL” are registered trode- 
marks of Digitol Equipment Corporation. 
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DIGITAL ANALYSIS AND CONTROL IN 
THE VIBRATION LABORATORY 


Peter Moseley 
Hewl ett-Packard 
Digital Signal Analysis Group 
Santa Clara. California 


A discussion of the differences between the traditif'r.u'i 
analog Instrument and the modern mini-computer based 
digital Instrument Is given. Examples of the kinds of 
measurements and new capabilities provided by a digital 
instrument are presented. Included is a brief discussion 
of how vibration control is implemented In a digital 
Instrument/system. Finally, a discussion is presented 
on the Impact of digital instruments on the typical 
vibration test laboratory. 


INTRODUCTION 


The advent of the mini computer has allowed the 
basic concept of an instrument to change, and 
this is why digital-based instrument/ systems 
are having such a dramatic effect on the vibra- 
tion test laboratory. 

This paoer summarizes the Important difference 
between analog and digital instruments and 
discusses some of the advantages of the digital 
approach. The use of memory to save data and 
numerical computation to make the measurement 
allows a single instrument/ system to perform 
a wide variety of tasks in the vibration 
laboratory. These capabilities range from PSD 
analysis to transfer functions (and therefore 
modal analysis) to complete random, sine and 
transient control. 

Analog and Digital Instruments 

Figure lisa block diagram of a typical 
analog instrument. 



The instrument includes input circuitry, front 
panel controls and some form of display; for 
exanple, a meter movement, numerals or a 
cathode ray tube (CRT). The measurement Is 
actually performed using fixed purpose analog 
circuitry which may typically be made to operate 
over several ranges, but still performs only 
one fundamental type of measurement. More 
sophisticated analog instruments provide "plug- 
ins" to extend the capability, and offer 
various outputs to transfer the measured 
information to another device; e.g., a 
recorder, printer or a computer. 

Figure 2 is a block diagram o^ a typical digital 
Instrument. 



Superficial ly , it looks much like an analog 
instrument, but there are many important 
differences both for the manufacturer and the 
user. The input circuitry consists of anti- 
aliasing filters and analog-to-digital con- 
verters. The front panel controls usually are 
no longer knobs, but more generally are push- 
buttons. These can ranae from alphanumeric 
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keys on a CRT terminal to pushbuttons on a con- 
trol panel (see Figure 3). 



FIbURE 3 

A Digital Instrument/Syste^i Control 


The display can range from a CRT to a full 
graphics capability terminal. The auxiliary 
outputs can Include complete graphics capa- 
bility. 

The fundamental difference between analog and 
digital instr indents is In the way a measurefiient 
1* made. Provided the a nit-aliasing filters 
are properly set* the data that is gathered 
by the A/D convener and stored In the instru- 
ment nienory represents all the Information 
available about that signal during the time 
interval 1i was observed. 

Once the data has been sampled, the software 
In the instrument detem^ines the form in u.nrh 
the information is presenter to the user, i,e,* 
the "measurement.'* The power of the digital 
instrument lies in the fact that once the data 
has been properly sampled virtually an> measure- 
ment ma> be 

The Infonm-Lioi: tret is available is detennln^ 
;>y the sampli.^g a’(* *ers used to set the AdC 
whkn In turn Is go^ ^ned by the relationships 
showr In 4, Hit:,? relationships are 

laws njcurr do v.'t v„ry b?tv-een instru- 
p<?ut ■'lanuuct Ji nor .'e^ween vibration 
laboratories. 


The speed of the analog-to-digital ^.onverter 
deternHnes the bandwidth of the instrument, and 
the number of bits {and the computer word siie) 
determine the dy ramie range* For example* 
with a 10-bit ADC and a 16-bit word size a 
narrow band power spectrum may be computed 
with an 80 dB dynamic range. 

Table I lists a representative set of measure- 
ments which may be made from the keyboard of 
a digital instrument. All of these represent 
the same information transformed by one pro- 
cess or another to the most useful funn. 


Another important difference between analog 
and digital instruments Is the use of 
digital memory. This allows the original data 
to be saved and the 'measurement'* to be per- 
formed after the time waveform Itself is gone. 
Further* many "measurenents” may be made on the 
same data. For example* it is possible to 
compute, among other things* the histogram, 
power spectrum, auto correlation function, and 
rms power of a single set of random data. To 
make a similar number of measurements using 
analog instruments would require many ins+ru- 
mnnts making measurements either simuUanecusly 
or one after the other on recorded data. 

Speed, Bandwidth , Overlapp ed 
P rocessing and "Realtime 

Ocassionally, it is desirable to make the 
'measurement" (1,e., sample the data and 
perform the computations necessary to transform 
the data to the desired form and display it) 
in "^eal time"; that is, without loss of date 
over many contiguous observation times T. In 
this case* the performance of the machine- is 
limited by the amount of information (band- 
width) 4 nd processing speed of the instrument. 
Processing speed depends upon the computer cycle 



TYPICAL MEASUREMENT CAPABILllV 
OF A DIGITAL UiSTPiiHEMT/SySTEM 


TABLE I 


RMS Power 

Peak 

Average 

Linear Spectrum 
Auto Power Spectrum 
Auto Correlation 
Histogram 
Integration 

'entlation 

Bai.w lelectahle Fourier Analysis 
1/N Octave Analysis 
ShocF Spectrum Analysis 
Amplituue i^odulatlon 


Frequency liodulation 
Phase Modulation 
Hilbert Transforms 
Cross Power Spectrum 
Cross Correlation 
Transfer Function 
(Transmissibil Hy, frequency 
response, impedance) 

Impulse Response 
Convolution 

Coherence (signal-to-noise) 
Multiple Coherence 
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time and upon the organization of the software. 

In machines equipped with microcoding capability 
up to a 10 to 1 speed improvement may be obtained 
over conventional machine language code by 
implementing certain highly repetitive opera- 
tions in microcode(l). 

When making measurements o.. random vibration 
data it is desirable to get the maximum vari- 
ance reduction per unit of time. The greatest 
variance reduction comes when the data frames 
are overlapped prior to performing the Fourier 
Transform. Overlapping capability is particu- 
larly useful in analyzing low-frequency random 
data, i.e., less than 500 Hz, because the 
record time T becomes relatively long and to 
get the desired variance reduction a great 
many averages may be required. By over- 
lapping, the time to get the eouivelent 
variance reduction is significantly re- 
duced. For example, if the measurement 
bandwidth is 50 Hz and ICO averages are 
required, overlapping can provide a 30 to 1 
reduction in measurement time. 

As the bandwidth is increased, t'le processing 
time becoMes a larger and larger fraction of 
the data input time T. When they ar - t 

equal, this is the "real time" limit Mi 

machine for the given measurenent. Beyond 
this limit data is lost. For example, the 
real time limit for a typical Fourier Analyzer 
is 3000 Hz when performing a power spectrum 
average. 

Mass Memory 

If the process under observation can be con- 
sidered stationary, then he loss of data 
above the r^l time lim.t is of little conse- 
quence unless there was a limited amount of 
it to begin with and not enough averages can 
be made to achieve the desired variance reduc- 
tion. In this case, some alternative pro- 
cedure is desirable to capture all of the avail- 
able data. The data must be put into a mass 
mer'^^V large enough to hold the data. This 
can be done if the digital instrument has the 
capability to digitize and send the data to a 
magnetic disc or tape in "real time" that is 
without loss of data. This capability is 
referred to as "ADC throughput." Once the 
da^a is on the disc or tape, then the measure- 
ment can be made after the fact using all of 
the data. 

The bandwidth limitation is determined by the 
throughput rate to the disc or magnetic tape. 

A typical specification tor this is 39 KHz 
throughput rate ^>er channel, which means that 
single-channel measurements can he made over 
a bai.dwidth of 19 KHz, and dual-channel over a 
$ 'z bandwidth. 


A digital instrument equipped with a magnetic 
disc or tape not only can have throughput 
capability but also general data and program 
storage capability. This extends the useful- 
ness of the machine and provides, among other 
things, the means for establishing a digital 
data base for such applications as modal 
analysis and signature analysis. 

Most vibration laboratories record test data 
on analog magnetic tape and perform detailed 
analysis after the test is complete. Al- 
though digital magnetic discs and tapes are 
extremely useful tools in conjunction with 
a digital instrument, they do not replace the 
function of analog tape recoi dings of test 
data. These records are often just that— 
records— and they need to be kept in an in- 
expensive form. In addition, to digitally 
record, say, 200 channels f data from a 
feur-minute spacecraft vibiation test out to 
a bandwidth of 2.5 kHz is not, as yet, prac- 
tical, as it would require throughput rates 
In excess of 2 megawords/second to record on 
a single device. If parallel devices are 
used, it is feasible, but expensive. In 
addition, if phase information is to be main- 
tained between all the channels, then they 
must be simul caneously sampled. This would 
imply multiple multiplexers slaved together 
and, again, it becomes expensive. 

For modal surveys, nowever, where the test 
levels are lower and the test duration can be 
lengthened, then the digital throughput 
method works extremely well as a means of 
gathering and maintaining the data in an 
organized and readily accessible fc^ i. 

Some Considerations in Making 
Measurements with Digital -Sampled 
Data Instruments 

As noted before, almost any measurement can be 
performed on he sampled data. The measurements, 
to be properly made, must rest on a strong 
theoretical background. Yet, the user of the 
instrument should not need to be intimately 
familiar with th*^ theory, so long as he can 
properly make the measurements he needs. 

Fixed Resolut jn Bandwidth Measurements 

A practical problem in this regard comes up 
when measurements are made on random data. 

For example, it is an easy matter to measure 
the rms power of a signal using digital 
techniques and with an ordinary true rms meter. 
But what if they disagree? VJhich instrument 
is right? Most often they are both right. The 
difference is the bandwidth of the instrument. 

The true rms meter might have a bandwidth from 
10 Hz to 1 MHz, and the digital system might 


be set to have a bandwidth frocr DC to 2 KHi. 
Clearly, the digital system will measure less 
rms power. 

This same difficulty arises In ccMnparlng spectrum 
analyzer outputs. There are often differences 
between wave analyzers, "speed-up'* analyzers 
and digital FFT-based Instruments* The 
differences almost always come back to a failure 
on the part of the person making the measurement 
to make the measurements using comparable resolu- 
tion bandwldths and properly normalizing the 
results* 

Figures S and 6 show two measurements of a sine 
wave with harmonics superimposed on broadband 
random noise* The measurement In Figure 5 was 
made using a spectrum analyzer which uses analog 
sweep techniques. Figure 6 shows the same meas- 
urement made using a Fourier Analyzer which 
uses FFT techniques* The measurements were not 
made to quite the same overall bandwidth noi 
to the same resolution bandwidth* 

Figure 7 was made with the analog analyzer 
using a 3 Hz bandwidth and now we see how 
closely the two agree. 

The Important point is that the different 
ft Her bandwldths are performing different 
convolutions on the data and that it is the 
area under the curve that gives us the 
Information, not the peak . Notice that the 
peaks In Figures 5 and f differ considerably, 
but the areas are approximately the same. 

Variable flesolution Bandwidth 
Measurements 


Many specifications such is MIL 810C allow 
random data to be measured using two different 
bandwldths, say, 25 Hz from 20 to 2000 Hz and 
50 Hz from 200 to 1000 Hz, 

What do you do if you*' digital instrument can 
o^'ly give you constant bandwidth measurements? 
In general, tnere Is no difficulty getting 
enough resolution, but sometimes more than 
enough Is constdered too much, A very narrow 
line containing only a small amount of energy 
may peak up beyond the test limits when anal- 
yzed with more resolution than specified. 

There are several solutions to this problem* 
One is to measure the data twice using two 
different bandwidths and merge the data to- 
gether to form a composite spectrum. Another 
is to measure the data with a single measure- 
ment of narrower bandwidth than required and 
then apply a smoothing function to the lower 
frequency portion to effectively broaden the 
filter shape* Other techniques are possible 
using digiUl filtering methods* 


FIGURE 5 

Analog Spectrum Analyzer 
Resolution Bandwidth - 30 Hz 


FIGURE 6 

Digital Spec t nan Analyzer 
Resolution Bandwidth = 1*02 Hz 


FIGURE 7 

Analog Spectrum Analyzer 
Resolution Bandwidth » 3 llz 








This example emphasizes the unique capabilities 
of a digital instrument to "massage" data to 
suit the needs of the user. However, this 
capability must not be incorrectly or improperly 
used or erroneous results may be obtained. 

SOME NEW CAPABILITIES 

Transfer Functions 

The concept of a transfer function is old, but 
its use in the vibration test laboratory is 
still in its infancy. Partly, this sttms from 
the fact that the traditional method', for 
making transfer function measurements are 
quite slow. A transfer functi^* computed using 
baseband Fourier analysis teciir ques (that is, 
measurements which cover the range from OC to 
some maximum frequency) can be riu le in a matter 
of seconds and can provide a great deal of 
useful information. 

An example of this is the use of low-level 
random noise excitation to measure the trans- 
missibility between critical parts of a space- 
craft before and after sinusoidal qualification 
testing. A comparison of the data before and 
after quickly reveals if the structure has 
changed and at what frequency. This technique 
has pinpointed trouble that might possibly 
have escaped detection on several occasions. 

The reason the technique is fast is that a 
broadband stimulus is used to make the measure- 
ment and all the frequencies are measured at the 
same time. The stimulus can be random or 
pseudo-random or even a transient(2). Small 
hammers with load cells mounted on them make 
ideal "impulse generators" and are now avail- 
able commercially. 

The use of transients for making transfer func- 
tions has several advantages. First, it is 
easy to set up and, second, the measurements 
can be made very rapidly. Complex shaker fix- 
turing is not required, and it provides 
a very good "quick-look" capability. A 
disadvantage comes when the structure under 
test is very non-linear or requires a carefully 
controlled force input. Often, the problems 
with non-linearities may oe overcome by pre- 
loading the structure. More often than not, 
controlled force is not required because of 
the large dynamic range capabilities of the 
digital instrument. 

Another example of the utility of easily made 
transfer functions is f'xture verification. 

In general, fixtures are linear and pseudo- 
rardom noise Cc. i be efficiently used to excite 
the fixture and measure its t»'ansmissibility 
and verify 'he design. 

Modal Analysis 

The most general use of t\ ^n.>‘*er function 


measurements comes in modal analysis* which 
is discussed in detail in References (3), 

(4) and (5). There it is pointed out that 
the transfer function contains all of the 
information necessary to identify the modal 
parameters of a system; that is, the resonant 
frequency, damping, and complex residue. 

Further, a matrix of transfer functions can be 
used to completely model a mechanical system. 

By examining the complex residues of a series 
of transfer functions, the mode shape vectors 
can be obtained and displayed in animated 
form. 

Resolution 

One 0 ^ the objections raised about digital 
Fourier transform-based transfer functions is 
resolution. The traditional transform is a 
baseband measurement; that is, it provides 
uniformly spaced frequency samples from DC 
to some maximum frequency. What do 'ou do if 
you need resolution of a millihertz at a 
kilohertz? Baseband digital methods would re- 
quire an enormous data record and would take a 
long time .o process a lot of unwanted informa- 
tion. A new solution is Band Selectable Fourier 
Analysis. 

Band Selectable Fourier Analysis 

Band Selectable Fourier Analysis, or BSFA, is 
a measurement technique that makes it possible 
to perform Fourier analysis over a frequency 
band whose upper and lower frequencies are 
selectable. BSFA provides improvement in fre- 
quency resolution greater than 100 times 
compared to standard Fourier analysis without 
increasing the number of spectral lines stored 
in the computer. The dynamic range is main- 
tained at 80 dB when the major signal is with- 
in the selected band and is increased to 90 dB 
or more when the major signal is outside the band. 
The dynamic range is considerably greater than 
that obtained with analog hetrodyne techniques 
because of the sharp filters that can be imple- 
mented digital ly (6) . 

The way BSFA is achieved is to multiply the 
sampled data by a sine and cosine wave at the 
selected band's center frequency and then 
digitally filter the data, resample it and 
Fourier tran:iform the result (see Figure 8). 

The sine, cosine multiply produces a complex 
time waveform shifted so that the center fre- 
quency is now at DC; the digital filtering 
removes the unwanted infonnation outside the 
selected bands; and resampling, from the filtered 
data stream, causes the maximum frequency of the 
new data to be one-half the bandwidth of the 
selected band. The Fourier transform (now com- 
plex) then coverts the information to the fre- 
quency domain. Once in this domain, power 
spectral density and transfer functions can be 
computed. It should be pointed out that the 
sampling laws of nature described above still 




Consider whai is happening in the freguency domain after 
each of the operaltons shown above The region of interest 
IS highlighted by giving it a f1 shape. The femainder of the 
spectrum is assumed to be flat 
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This signal is then Fourier transformed into N equally spaced 
spectral lines between -Fj“2 arnS F*v2. 
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Important BSFA Relations 

M =Ft^'N - F, 'n N = I T 


Where 


Ft' ” Resampling Frequency 
N = Number of Frequency Lines 
Ft “ Original SampNr>g Frequency 
n 5= Int^er Determined by Ft and BSFA Bandwidth 
T' ^ Length ol Time Record 


FIGURE 8. BSFA signal flow. A frequency shift 
by ff. centers the analysis band on the frequency 
of interest and digital filtering removes in- 
formation outside the band of interest. 



85 






apply and that to achieve fine resolution still 
requires a long meajjureitiefit time. For example, 
to obtain 1 mil li hertz resolution requires 1000 
seconds of data. Figures 9 and 10 show the ad- 
vantages of BSFA in power spectrum measurments. 
Figure 11 shows a twc-channel transfer function 
BSFA measurement. 





Frvquancv (Mil 


FIGURE 11, In this mechanical impedance mea- 
surement on a squeaky automobile disc brake 
rotor, BSFA made possible an accurate 
determination of tht- r'vies of vibration and 
the dampl ’^9 MSSOciatU rfith each mode. 


With BSFA capability, it -s now an easy matter 
to make a quick baseband measurement, decide 
where to look more closely and thef^ zoom in 
on it. (Hence, the name Zoom t^ransfom is 
sometimes used to describe BSFA.) Ho longer 
is resolution a deterent to the use of the 
transfer function In the vibration test labora- 
tory. Indeed, the BSFA technique is essential 
to successful modal analysis using broadband 
transfer function techniques. 

The BSFA Tneasurement 'an be made on data saved 
on a mass memory using ADC throughput capa- 
bility. This means that one set of data can 
re-examined with any resolution up to that 
dictated by the record length itself. 


1/Nth Octave Analysis 

Just as band selecUble Fourier analysis makes 
use of digital filtering techniques, so does 
\/n octave analysis. 1 /H octave analysis is 
a measurement technique that utilizes convolu- 
tional digital filters to separate the spectral 
components of a sampled waveform into 1/N 
octave spaced bands of 1, 1/2, 1/3, 1/6 and 
1/12 octave spacing. Because of the digital 
nature of the Instrument, the information is 
stor 0 d in manory in 1 / 1 2 th octave and the 
final result may be viewed in any of the 
octave s pacings with any standard weighting; 
i.e. , A, 8 , C or D, 

Figure 12 shows a plot of a typical l/12th 
octave measurement performed on acoustic data. 
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Shock Response Spectrum 

Another measurarent that can be implemented 
using digital filters is the shock response 
spectrum. Unlike the Fourier transform, this 
Is a non-linear, irreversible process. 

Historically, It predates even most electronic 
analog instruments, having come from Biot‘s 
efforts to determine the frequency spectrum of 
earthquakes using tuned mechanical "filters." 
This "poor-man’s" spectrum analyzer has been 
enulated in both electronic analog and 
digital instruments, but seems archaic now 
that so many more powerful tools are avail- 
able (7). However, its use persists and the 
fact that It can be implemented in a digital 
instrument demonstrates the versatility of 
such devices. 

Special-Purpose Keasurements 

As should be clear by now, the measurement 
capability of a digital instrument system 
lies in the software that implements the 
Measurement, If a system has mass storage 
capability such as a disc then a great deal of 
software can readily be loaded into the computer 
to change the system capabilities. Special- 
purpose measurement packages such as signatu**e 
analysis and modal analysis can be implefnented 
this way. These make use of the basic measure- 
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8SFA SPECTRUM 30 Hi/DlV 


B/tSEBAND SPECTRUM 2.5 kHz/DIV 


5010 Hi 


FIGURE 9. The resolution and dynamic range improvenents of 
Band Selectable Fourier Analysis (DSFA) are demonstrated 
by these power spectra of the output of a sine wave 
generator. The conventional baseband spectrum shows the 
5-kHz fundamental and the third harmonic at -69 dB. The 
noise f’ 'Or is 72 dB down. The BSFA spectrum shows that 
the funujmental is actually at 5010 Hz and that there are 
60-Hz and 120-Hz sidebands. The noise floor is greater than 
80 dS down. 
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merit capabilities described above, but present 
the capabilities to the user at a much higher 
level. Thus, a single button push can produce 
an RPM spectral map or an animated mode shape. 

Vibration Control 


As if all of the above measurement capabilities 
aren't enough, a digital instrument/system can 
provide vibration control capability including 
sine, random and transient (shock) control. 
Figure 13 shows how a digital analysis instru- 
ment is turned into a digital vibration control 



FIGURE 13 

Block Diagram of 
Digital Vibration Control 


In the Hewlett-Packard Vibration Control Sys- 
tem, the generation of the random signal (8) 
and the sinusoidal signal are both implemented 
making use of the microcode capability of the 
machine. 

The use of microcode to keep the signal genera- 
tion overhead to a low level, as well as per- 
form several other time critical operations, 
makes it possible to implement 4-channel 
random control and up to 12-channel sinusoidal 
control. Four-channel random control is per- 
formed by averaging four power spectra computed 
from simultaneously sampled time waveforms using 
the system ADC. Twelve-channel sine control 
requires a multiplexer and includes peak 
select in notching capability. 

Figures 14, 15 and 16 show post-test plots of 
random, sine and transient tests run on a 
structure having the acceleration-to-amplifier 
voltage, transfer function shown in Figure 17. 



An nalog-to-digital converter and a filter, 
along with special software, are added to give 
signal generation capability to the system. 

Although the basic elements are the same as 
those available in the general purpose 
instrument/system, the uses and requirements 
of a vibration control system are different 
and more emphasis needs to be given to simple 
operator controls to avoid mistakes and clarify 
the operation. The system control keyboard is 
used for this purpose (see Figure 3). 

If the vibration control system is disc-based, 
changing from Fourier analysis to random to 
sine or to trans"»ent control can be almost 
instantaneous. 

The successful implementation of digital 
signal generation and control over 5 KHz band- 
width requires a high-performance computer. 

The overhead required to generate the signals 
takes away from the processing time available 
for measurement and control and, if it becomes 
too great, makes for a poor control situation. 


Random Control Spectrum with 
Test Limits Superimposed 



Sine Control Envelope with 
Test Limits Superimposed 
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Transient Waveform 
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Transfer Function of Device Under Test 


Perhaps the most remarkable contribution of 
digital vibration control is the time saved 
conducting tests. All or the "knobs" are set 
by operator responses to questions, and these 
test parameters can be recalled and set up in 
a matter of seconds. The high accuracy and 
repeatability of digital vibration control 
greatly simplifies test operation by, in most 
cases, eliminating "equalization runs" on 
bare fixtures. 

Some of the other advantages of digital vibra- 
tion control over conventional analog tech- 
niques are increased resolution and dynamic 
range in random tests; extremely flexible 
test strategies including multiple-channel 
notch control in sinusoidal testing; and 
shock spectrum analysis and synthesis in 
transient waveform control. 


A modern digital instrument/ system is a truly 
powerful tool in the laboratory, but if it is 
to be used to its full capabilities it must be 
managed by people who understand the capabili- 
ties and limitations and who have the skills 
to implement the solutions to their unique 
measurement problems. This represents a cost 
after initial purchase which is sometimes 
overlooked ana may result in underutilization 
of a valuable resource. The written documenta- 
tion and training provided by the manufacturer 
are an important part of what is purchased. 
Therefore, the maintenance of the software 
and manuals in an organized and systematic 
manner is necessary to protect the investment. 

Since most of the measurements listed in Table 
I can be performed using the keyboard pro- 
gramming language, there is ofU^n little need 
for the user to write special software. How- 
ever, the capability Is there. If it is used, 
careful documentation is required tc make the 
results useful to others. 

CONCLUSION 

The broad range of measurements that can be 
made, coupled witn the vibration control 
capabilities, make a digital Fourier trans- 
former-based instrument/ system a necessity in 
a modern vibration laboratory. When properly 
managed and utilized, production can increase 
many times on routine work. Valuable measure- 
ments that used to be impractical, and were 
overlooked, can now be made. 
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